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Abstract 

In this review some recent multi-loop results obtained in the framework of perturba- 
tive Quantum Chromodynamics (QCD) and Quantum Electrodynamics (QED) are 
discussed. After reviewing the most advanced techniques used for the computation 
of renormalization group functions, we consider the decoupling of heavy quarks. 
In particular, an effective method for the evaluation of the decoupling constants is 
presented and explicit results are given. Furthermore the connection to observables 
involving a scalar Higgs boson is worked out in detail. An all-order low energy 
theorem is derived which establishes a relation between the coefficient functions in 
the hadronic Higgs decay and the decoupling constants. We review the radiative 
corrections of a Higgs boson into gluons and quarks and present explicit results 
up to order and af, respectively. In this review special emphasis is put on the 
applications of asymptotic expansions. A method is described which combines ex- 
pansion terms of different kinematical regions with the help of conformal mapping 
and Fade approximation. This method allows us to proceed beyond the present 
scope of exact multi-loop calculations. As far as physical processes are concerned, 
we review the computation of three-loop current correlators in QCD taking into 
account the full mass-dependence. In particular, we concentrate on the evaluation 
of the total cross section for the production of hadrons in e~^e~ annihilation. The 
knowledge of the complete mass dependence at order has triggered a bunch of 
theory-driven analyses of the hadronic contribution to the electromagnetic coupling 
evaluated at high energy scales. The status is summarized in this review. In a 
further application four-loop diagrams are considered which contribute to the order 
QED corrections to the ^ decay. Its relevance for the determination of the Fermi 
constant Gp is discussed. Finally the calculation of the three-loop relation between 
the MS and on-shell quark mass definitions is presented and physical applications 
are given. To complete the presentation, some technical details are presented in the 
Appendix, where also explicit analytical results are listed. 
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1 Introduction 



Nowadays the Standard Model (SM) of elementary particle physics is well established. 
Some parts of it (e.g. the properties of the Z boson) have been tested to an accuracy far 
below the percent level — mostly at the CERN Large-Electron-Positron collider (LEP, 
Geneva), at the SLAG Linear Gollider (SLG, Stanford) and at the Fermilab TEVATRON 
(Ghicago). Up to now no significant deviation between theory and experiment has been 
found. For other parts of the SM, related to GP violation and quark mixing, the B 
factories like BaBar at SLAG, Belle at KEK (Tsukuba) or HERA-B at DESY (Hamburg) 
will provide deeper insight, and significant improvements in the determination of the 
corresponding parameters will be obtained. Gurrently mainly the scalar sector of the 
SM eludes from direct experimental observation. This affects both the generation of 
the particle masses and the existence of the Higgs boson itself. Probably Run II of the 
TEVATRON and certainly the Large Hadron Gollider (LHG) at GERN will provide more 
sureness — not only in connection to the Higgs sector but also to possible extensions of 
the SM. Once the Higgs boson would be discovered it immediately would become subject 
to precision measurements. In particular at a future e^e~ linear collider such as DESY 
TESLA, a precise study of its properties would be possible. 

In the recent years there has been an enormous development in the evaluation of 
radiative corrections. It is fair to say that the major part of it was initiated by the 
fundamental works of 't Hooft and Veltman in 1972 |^, Q where dimensional regular- 
ization (see also ||^) was established as a powerful tool in the evaluation of multi-loop 
diagrams^. Since that time a whole industry has been formed to develop techniques for the 
computation of complicated Feynman integrals. At one-loop order the procedure of the 
computation has been systematically studied quite some time ago p, |^, ||, Q . Neverthe- 
less also nowadays it is not completely straightforward to evaluate an arbitrary one-loop 
diagram — in particular if many legs and complicated momentum configurations are in- 
volved. One can easily imagine that at two and more loops one arrives quite soon at the 
limit where the occuring mathematical expressions can not be solved. At two-loop order 
certain classes of diagrams can still be treated by either using a combination of analytical 
simplifications and fast numerical routines, like in the case of two-point function with 



several non-zero masses |jTO|, or applying purely analytical methods, like in the case of 
massless digrams with four external legs fill]. However, at three-loop order it is essentially 
only possible to solve one-scale integrals. A systematic study at four or more loops is still 
missing. 

QGD, the field theoretical realization of the strong interaction, constitutes an impor- 
tant part of the SM and also of most of its extensions. At low energies the coupling 
constant of QGD, a^, is large and perturbative calculations are not possible. However, 
due to the phenomenon of asymptotic freedom the value of gets smaller with raising 
energy and perturbation theory is an appropriate tool to evaluate radiative corrections. 

Up to now the vast majority of the multi-loop calculations have been performed in the 



^Dimensional regularization applied to infra-red divergences and mass singularities has first been 
considered in p. 
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framework of QED and QCD. One reason is certainly that the calculations are simpler as 
compared to the full SM since there are less parameters. Furthermore, there is a strong 
hierarchy both in the quark and lepton masses which also simplifies the calculations. On 
the other hand, the higher-order corrections are indeed necessary. In QED there exist 
precise experiments which require high theoretical precision and although the coupling 
constant is quite small sometimes high loop orders are necessary. For example, in the case 
of the anomalous magnetic moment of the electron, four-loop corrections are needed to 
match the experimental precision. In QCD the coupling is roughly a factor of ten bigger. 
Nevertheless it is often still small enough to perform a perturbative expansion. However, 
the higher order terms are significant and can not be neglected in the cases where high 
precision is required. 

In this work some recent developments in the calculation of multi-loop diagrams are 
reviewed. Thereby we will mention the most important methods which have been used in 
the computation of higher order quantum corrections and explain a few selected ones in 
greater detail. At the same time, we discuss the present theoretical status of important 
physical quantities. In particular, the renormalization group functions in the modified 
minimal subtraction scheme |1T2| , |13[ (MS) are provided up to the four-loop order. As 
is well-known, four loop running must be accompanied by three-loop matching at quark 
thresholds. The corresponding decoupling relations are presented in Section 0. 

The hadronic Higgs decay is closely connected to the decoupling relations as we will 
show is Section ^ Parts of the quantum corrections can be computed in the framework 
of an effective Lagrangian where the coefficients can be determined from the decoupling 
relations. The origin of this miraculous connection lies in the use of the dimension-four 
operators, which constitute an important ingredient of the effective Lagrangian. Another 
application of the dimension-four operators are the quartic mass corrections to the cross 
section a{e~^e~ hadrons), which is also discussed in Section |[ 

The last issue is again picked up in Section H, where also QCD corrections to the 
production of hadrons in e+e~ annihilation are computed. Putting together all terms one 
arrives at a complete picture up to the quartic mass corrections of order a^. The main 
purposes of Section ^ are practical applications of asymptotic expansions^. Besides the 
diagonal correlators also the non-diagonal ones are considered. 

The Fermi constant, Gp, the mass of the Z boson, Mz, and the electromagnetic 
coupling constant, a, are the best known parameters of the SM. Mz has been measured 
at LEP with an accuracy of a few per mille to be Mz = 91.1876 ± 0.0021 GeV. In this 
review we want to discuss quantum corrections to the other two parameters. An essential 
ingredient to the running of the electromagnetic coupling from = to = M| is 
provided by the cross section a{e~^e'' hadrons). The correction terms discussed earlier 
have been used to obtain so-called theory-driven results. The different approaches are 
discussed. As further applications we present the status of the QED corrections to the 
muon decay and the relation between the MS and the on-shell quark mass. Let us in the 
following discuss the individual issues in more detail. 



■See Appendix A.l for details 
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As far as radiative corrections are concerned, a crucial role is played by the renor- 
malization group functions. In particular the functions (3{as) and 7m(«s) governing the 
running of the coupling and the quark masses comprise a significant part of the higher 
quantum corrections. They are in particular very important to re-sum large logarithms 
to all orders in perturbation theory. Only a few years ago the four-loop terms of order 
have been evaluated for and 7rrt(tts) US], |1^]. The latter has been computed by 



two groups using completely independent methods. Both methods are based on the fact 
that the pole part of a logarithmically divergent diagram is independent of the masses or 
momenta. In [l^ this is exploited together with the technique of infra-red re-arrangement 
(IRR) in order to obtain a factorization of the four-loop integrals into massless three-loop 
and massive one-loop ones. In Refs. 0, |16[ all lines were assigned to the same mass, and 
all external momenta were set to zero, which leads to a special class of bubble diagrams. 
From them only the pole parts have to be computed. In Section |^ we want to review both 
methods and explicitly demonstrate the way they work. 

In this review special emphasis is put on the construction of effective theories in the 
framework of QCD. In Section ^ an effective QCD Lagrangian is constructed for the case 
where one of the quarks is much heavier than the others. The construction is made explicit 
by specifying the relations between the parameters in the full and effective theories. These 
relations provide at the same time the well-known decoupling constants which have to be 
applied in QCD every time a particle threshold is crossed. The most prominent example 
for their necessity is probably the computation of as{Mr) from as{Mz) or vice versa. In 
the latter case five quarks are active whereas in the former one only three quarks are 
present in the effective QCD Lagrangian. 

In Section ^ a slightly different point of view is adopted. Here the scalar operators 
of dimension four are considered in QCD. In a first step they are used to construct 
an effective Lagrangian describing the coupling of an intermediate-mass Higgs boson to 
quarks and gluons. The top quark is considered as heavy and manifests itself in the 
coefficient functions of the effective Lagrangian. Once the latter has been found, the 
imaginary part of the Higgs boson correlator in the effective theory leads to the total decay 
rate. As a central result of Section ^ we derive a low-energy theorem which considerably 
simplifies the computation of the coefficient functions as they are related to the decoupling 
constants of QCD evaluated in Section ^. In Section [4.1.3| the background field method 
is introduced as a convenient tool for the computation of higher-order corrections. As an 
example, the coefficient functions describing the decay of the Higgs boson into gluons is 
also computed in this framework. 

In the second part of Section ^ another important application of the scalar dimen- 
sion four operators is discussed, namely the quartic mass corrections to the cross section 
cr(e+e~ hadrons). Mass corrections of order [m? / s)'^ and [rn? / sY are obtained rel- 
atively easy as in QCD there are no non-trivial operators of dimension less than four. 
However, the quartic corrections require the inclusion of the dimension-four operators 
with all their renormalization and mixing properties. We will explain the techniques and 
present results obtained recently at order a^. 

The last part of this review. Section ^ is devoted to the discussion of some results 



4 



obtained with the help of asymptotic expansion accompanied with conformal mapping 
and Pade approximation. This method has been developed in the recent years and has 
been applied successfully to a number of important processes. The underlying idea is the 
following: only in rare cases it is possible to compute three-loop diagrams involving more 
than one scale. However, if a certain hierarchy exists between the scales it is promising 
to apply an asymptotic expansion. This effectively reduces the number of scales present 
in the integrals which are subsequently significantly simplified. 

In particular we will discuss the corrections of order to the photon polarization func- 
tion. Its imaginary part is directly connected to the physical quantity R{s) = a{e^e~ — > 
hadrons)/(T(e+e~ — > /i"*"//"). The apphcation of our method leads to the full mass de- 
pendence. Combining the results with the quartic corrections given in Section [4.2| , one 
obtains a prediction for R{s) up to and including 0{a^gm^ / s^). Only recently also the 
non-diagonal current correlator formed by a massive and a massless quark has been com- 
puted. In this application special emphasis lies on the extraction of information about the 
threshold behaviour which has some relevance in the framework of heavy-quark effective 
QCD. 

As an application of the knowledge of R{s) to high perturbative order, we discuss the 
evaluation of q;(M|). The electromagnetic coupling is defined at vanishing momentum 
transfer. However, its evolution to high energies constitutes the dominant part of the 
radiative corrections to electroweak observables. The accurate determination of a(M|) 
is thus essential for any precise test of the theory. At the same time the indirect de- 
termination of the masses of heavy, hitherto unobserved particles, e.g. the Higgs boson 
or supersymmertic particles, depends critically on this parameter. Of particular impor- 
tance in this context is the hadronic vacuum polarization. It is nearly as large as the 
leptonic contribution, but cannot yet be computed perturbatively. However, it may be 
related through dispersion relations to the cross section for hadron production in electron- 
positron annihilation. The integrand can thus be obtained from data, phenomenological 
models and/or perturbative QCD, whenever applicable. In Section [5.2.2| we will discuss 
the developments in the evaluation of a(M|) which took place in the recent two to three 
years due to the knowledge of the complete mass dependence of R{s) at order (cf. 
Section |5.2.1|) . 

Gp is defined through the muon lifetime, and the decay of the muon, as a purely 
leptonic process, is rather clean — both experimentally and theoretically. The one-loop 
corrections of order a were computed more than 40 years ago whereas only recently 
the two-loop corrections of order have been evaluated [18, 19 1. The large gap in time 
shows that this calculation is highly non-trivial. The inclusion of the two-loop terms 
greatly reduced the relative theoretical error of 1.5 x 10^^ which was an estimate of the 
size of the missing corrections. The remaining error on Gp now reads 0.9 x 10~^ and is of 
pure experimental nature. Upcoming experiments will further improve the accuracy of the 
muon lifetime measurement and therefore the O(a^) corrections to the muon decay are 
very important and constitute a crucial ingredient from the theoretical side. In Section p.3| 
we discuss the results obtained with the help of asymptotic expansion. 

In the SM the quark masses have still relatively big uncertainties. This is mainly due 
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to the confinement property of QCD wliich prevents tlie production of free quarks. It is 
also important to have a convenient definition of the quark mass in order to perform a 
comparison between theory and experiment. Recently there has been quite some activity 
connected to the precise determination of the bottom- and top-quark masses. The bottom- 
quark mass is determined with the help of QCD sum rules where a proper mass definition 
helps to reduce the error. In the case of the top quark, studies have been performed for 
an e^e~ collider with a center-of-mass energy in the threshold region of top-quark-pair 
production. An energy scan which provides the measurement of the total production 
cross section would provide an error of about 100 MeV in the top-quark mass. Also here 
a special mass definition has to be employed. In both cases the three-loop on-shell-MS 
conversion formula is needed in order to obtain the corresponding MS quark mass. The 
latter is important for processes not connected to the threshold. In Section |5]^ these 
issues are discussed in detail. 



2 Renormalization group functions in QCD 

In perturbative QCD the renormalization group functions play a very important role. In 
particular the (3 and 7m functions governing the running of the strong coupling and the 
quark masses, respectively, are indispensable when evaluating physical observables. For 
this reason we decided to discuss the techniques used for the computation of the four-loop 
contributions to [3 |]14| and 7^ |]15|, |l^] in more detail. 



In the MS scheme the knowledge of the renormalization group functions is equivalent to 
the knowledge of the corresponding renormalization constant. Thus, in order to compute 
a renormalization constant at n-loop order it is sufficient to evaluate the ultra-violet (UV) 
poles of n-loop diagrams. Nevertheless it is often also quite useful to have also a handle on 
the infra-red (IR) poles. More details on the UV and IR structure of Feynman diagrams 



are given in Section |2.1| . Afterwards, in Sections p.2| and p.3| , two practical methods are 
described which have been applied at the four-loop level. 

2.1 Ultra-violet and infra-red counterterms 

Before we want to consider explicit examples in the next two sections the theoretical 
background needed for the higher-loop calculation of renormalization group functions is 
introduced in this Subsection. In particular we want to demonstrate the interplay be- 
tween UV and IR divergences in dimensional regularization and show their connection to 
asymptotic expansions. We refrain from presenting the material in a mathematical rigor- 
ous framework, for which we refer to the original literature, but exemplify the important 
points at explicit diagrams. 

As the properties discussed in this section are independent of the particle type we 
consider only scalar propagators of the form 
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where also M = is allowed. We furthermore assume that the integration momenta are 
denoted by ki, k2, ■ ■ ■ and introduce the abbreviation J- = J hi / {2t:)^ . 

In order to remove the UV divergences the so-called R operation has been intro- 
duced 



20, 21 



It is a recursive subtraction scheme where the UV divergences are re- 
moved from the Feynman integrals in a way compatible to adding local counterterms to 
the Lagrangian. Formally, the R operation can be written as a sum where each term is a 
product of operators acting on subsets of disjoint one-particle-irreducible subgraphs. In 
this way counterterms are generated. They have to be inserted into the vertices of the di- 
agrams which remain after shrinking the corresponding subgraphs to a point. Non-trivial 
diagrammatic examples can, e.g., be found in the text books p2|, |23|. The R operation 
applied to a one-loop propagator-type integral with external momentum q leads to the 
equation 



R 



[M2-A;f][M2-(A;i + g)2] 



1 



1 [M2 
finite , 



A;2][M2-(fci + g)2] 



(2.2) 



from which the renormalization constant Z^^^ is determined. In the MS scheme it reads 

1 / r 



167r2 



(2.3) 



Once Z^^^ is known one can turn to the two-loop order. The application of the R operation 
to the diagram shown in Fig. |2.1| leads to 



R 



1 J2 [M2 - {ki - A;2)2][M2 - {ki + g)2][m2 - A;2]2 



n J2 [M2 - {ki 
finite , 



A;2)2][M2-(A;i + g)2][m2 



k^V 



2 



which fixes Z'^2) to 



^(2) 





2 


■ 1 


1 ■ 


(l67r2. 


) 


.2^2 


~ 2e. 



zm 

(2.4) 
(2.5) 



A formal definition of the R operation can, e.g., be found in . 

There are two special features of dimensional regularization which prove to be a pow- 
erful tool, especially in the evaluation of renormalization group functions. The first one is 
that dimensional regularization is able to regularize simultaneously both UV and IR diver- 
gences. Furthermore, as was realized in p5|, in renormalization schemes based on minimal 
subtraction all UV counterterms are polynomial in the momenta^ and masses. This means 
that the divergence of logarithmically divergent Feynman integrals are polynomial in 1/e 
and there is no dependence on a dimensionful scale. As a consequence, in the computation 



^This has to be the case for each meaningful renormahzation prescription. 
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Figure 2.1: Scalar two-loop diagram. 

of the corresponding coefficients one is free to make arbitrary re-arrangements of masses 
and external momenta — provided no IR divergences are introduced. This was for the 
first time observed in Ref. (see also |2^, |24 



In the example of Eqs. ( p.2|) and ( p.4|) this means that the masses M could be set to 
zero without changing the results for Z^^^^ and Z^'^\ 

This procedure is often referred to as infra-red re-arrangement (IRR) [^. One has in 
mind a transformation of the IR structure of a diagram in such a way that its U V divergent 
part can be easily computed. Next to nullifying masses or momenta, also differentiations 
with respect to masses are allowed in order to achieve the simplifications. 

The requirement that no IR divergences may be introduced is quite restrictive and 
makes the application of the IRR very tedious. Often one has to remain with Feynman 
integrals which are not as simple as one would like to have them. To overcome this 
disadvantage the R operation has been generalized and the so-called operations has 
been developed ^ in order to deal not only with the UV but also with the IR 
divergences. Thus arbitrary re- arrangements of masses and momenta are allowed as the 
R* operation takes care of all occuring IR divergences. 

The i?* operation can be written as R* = RR = RR, where R is responsible for the UV 
divergences and R subtracts the IR ones. For R there exists a recursive definition which 
is in close analogy to the one of the R operation. The precise definition and the criterion 
for the subgraphs, which have to be considered while applying R, is quite involved and 
requires the introduction of a lot of mathematical terminology The way Z — the 



renormalization constant generated by i? — is computed in practice can best be seen by 
looking at the examples discussed below. At this point we only want to mention that a 
renormalization constant Z for the whole diagram only appears for scaleless integrals. As 
a consequence, the R operation does not commute with the limit of masses or external 
momenta going to zero (cf. Eqs. ( p.8D and ( |2.9| ) below). 

The R* operation is essential to prove the following theorem |]30| : 



Any UV counterterm of a {I + l)-loop Feynman integral can be written in 
terms of the poles and finite parts of appropriately constructed l-loop massless 
propagator-type integrals. 

This theorem is very powerful. It states that — at least in principle — all renormalization 
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group functions at four-loop order can be computed from the knowledge of massless three- 
loop two-point functions. The latter are, e.g., provided by the package MINCER [^. In 



practice there are problems connected to the (in general) large number of contributing 
diagrams as the prescriptions for the R* operation given in ||2^, ^ have to be applied 
individually to each of them. Thus, for practical applications further improvements are 
necessary. We will come back to this point later. 

Let us consider the one- loop integral of Eq. ( p. 21) with M = and q = 0. By definition 
the resulting massless tadpole integral is set to zero in dimensional regularization. On the 
other hand, the application of the R* operation gives 







R 
R 



-m-ki 
1 

1 



1 [-kl][-kl] 



(2.6) 



In the first step the R operation is applied resulting in the diagram itself and the countert- 
erm Z^^^ rendering the expression in the square brackets of the second line UV finite. The 
subsequent application of R generates the counterterm Z^^^ which corresponds to the IR 
divergence of one-loop integral. Once the application of R* is resolved scaleless integrals 
are set to zero which leads to the last line of Eq. ( p. 61) . Using Eq. ( p.6|) in combination 
with Eq. ( p^.3| ), one obtains for the IR renormalization constant of the one-loop two-point 
function 



1 1 



(2.7) 



In analogy, applying the R* operation to the diagram of Fig. |2]l| with all masses and 
external momenta set to zero leads to 







R* 



R 



1 J2 



1J2 [-k!][-k|n-{k^-k2)'] 



[-kl? 



+ ^(2) 



(2.8) 



This equation defines the IR renormalization constant Z^'^\ 

Let us come back to the two-loop diagram of Fig. with m = 0. In this limit 
the diagram contains both an IR and UV divergent subdiagram. The application of R* 
generates the following terms 



'1 J2 



[M2 - {ki - A;2)2][M2 - {ki + g)2][-A;i]2 
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II J2 [M2 - (k, - A:2)2][M2 - {ki + g)2][-A;2]2 



1 [M^ - kj][M^ - {ki + qy 




(2.9) 



and 



= finite , 

where R after the second equal sign acts on the IR divergent integral J2 1/ 
generates the factors Z^^^ . Note that no term Z(2) appears as the original integral involves 
the scales M and q. After the third equal sign the results of Eqs. (|2.3|), ( |2.5| ), ( p.7| ) and of 
the Appendix ( [A .41) have been used in order to explicitly show the finiteness. In general 
the logic is the other way around: one chooses the masses and external momenta in such 
a way that the very diagram, i.e. the first term in the second line of Eq. ( |2.9| ), can easily 
be evaluated and uses Eq. ( p.9|) in order to determine Z^'^\ 

At the end of this Subsection we want to work out the connection of the R* operation 
to the asymptotic expansion with respect to large masses. The hard-mass procedure 
provides a prescription on how to evaluate Feynman integrals where one of the internal 
masses is much larger than the others (cf. Appendix [A.l|) . Thus one can adopt the 



point-of-view to introduce light masses which regularize all IR divergences and apply the 
hard-mass procedure. For illustration let us consider the diagram in Fig. with one of 
the propagators carrying mass M doubled in order to avoid UV divergences which keeps 
the formula more transparent. For m 7^ the diagram is also IR finite. Applying the 
hard-mass procedure in the limit 
m2/M2 and g2/M2 leads to 

1 



M2 ^ m? , g2 and keeping only the leading terms in 



1 J2 [M2 - {ki - A;2)2][M2 - {ki + g)2]2[m2 



.212 



+ 



2 [M2 - {ki - k2f][M^ - klf[-klf Ji [M2 - kff J2 K - fci]2 



+ . . . 



(2.10) 



where the ellipses represent terms of order q^ /M"^ and jM"^ . The general rules and 
explicit examples for the hard-mass procedure are discussed in Appendix |A.1|. The ap- 



plication of Eq. (^]2) to the two-loop diagram at hand leads to the two contributions 
which are listed in the second line of Eq. ( |2.10| ) . The first term corresponds to the whole 
diagram which according to the rules of the hard-mass procedure has to be expanded in 
the small quantities q and m. To our approximation this means to simply nullify both q 
and m. In the second contribution the hard subgraph consists of the one-loop subdiagram 
where all lines carry the heavy mass M. The expansion in the external momenta leads to 
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Ji 1/[M2 — kf]^ which finally leads to the second term on the right-hand side of Eq. ( p.lO| ). 
Note that the first integral is IR divergent whereas the second one is UV divergent. Their 
connection to the R and R operation becomes clear after adding and subtracting the term 



Eq. (2.10) = 



1 



+ 
R 

+ 



1 J2 [M2 - {ki - fc2)2][M2 - fc2]2[-A;2]2 
1 / f 1 



1 



[M2 - A;2]3 



1 [M2 - A;2]3 \J2 [m2 - A;2]2 



[M2 - {ki - A;2)2][M2 - A;f]2[-A;i]2 



R 



[M2 - A;2]3 



2 



A;?12 



(2.11) 



where in the last step Eq. ( |2.6|) has been used. In Eq. ( p.ll| ) the correspondence between 
the IR and UV divergences, which are introduced through the hard-mass procedure, can 
nicely be observed. Furthermore, it can be seen how they have to be combined in order 
to arrive at the final form which contains the application of the R and R operations. 



2.2 Global infra-red re-arrangement and the quark anomalous 
dimension 

In the previous Subsection it has been demonstrated that the IRR in connection with the 
R* operation provides a very powerful tool for the computation of renormalization group 
functions at higher orders. In this Subsection we want to discuss its practical application 
in the case of the quark anomalous dimension 7^- 

In the MS scheme the running of the quark masses is governed by the function 7m(as) 
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[2.12) 



where the coefficients 7^,4 are known up to the four-loop order ^ |T5|, 
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(2.13) 



with Cs ~ 1.202 057, (4 = vr^/90 and (5 ~ 1.036 928. The superscript nj indicates the 
dependence on the number of quarks. 

For the computation of 7^ one has to know the quark mass renormahzation constant 
in the MS scheme, Zm, which relates the bare mass, m^, to the renormahzed one through 



m 



(2.14) 



Zm can be obtained form the vector and scalar parts of the quark propagator. Its inverse 
reads in bare form 



m 



m 



(2.15) 



where the functions and depend on the external momentum g, the bare mass m° 
and on the bare strong coupling constant From the requirement that the renormahzed 
quark propagator is finite one gets in the MS scheme the following two equations 



Z2 



Z2^V 



Z2Zm^% 



(2.16) 



where Z2 is the wave function renormahzation of the quark. The operator extracts the 
poles in 1/e. The Eqs. (|2.16| ) are solved recursively for Z2 and Z^, and 7^ is computed 
from Eq. ( p.l2| ) using ( p.l4| ) and the fact that /i^dmo/d/i^ = 0. This leads to 



In 



d 



-/3^(a,)7r7-— InZ^ , 
da. 



(2.17) 



where 13^ {as) = —e + /5(a<j) is the D-dimensional /3 function and [3{as) is defined below 

in Eq. (|2^ ). 

For the computation in the MS scheme one has to evaluate the pole part of the 
fermion propagator. Some sample diagrams contributing at one, two, three and four loops 
are shown in Fig. p.2| . As the diagrams contributing to and Sy are logarithmically 
divergent the computation can be performed by setting all internal masses to zero and 
keeping the external momentum finite. This leads to massless /-loop propagator-type 
integrals for the evaluation of 7^,^-1 • AH occuring integrals are free from IR divergences. 
At one-, two- and three-loop order the package MINCER can be used in order to perform 
the computation. However, this method fails to be practical for the computation of 7m,3 
as currently massless four-loop integrals are not available. Alternatively one could set 
the external momentum to zero but keep the quark mass finite. Again the technology is 
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Figure 2.2: Sample diagrams contributing to the fermion propagator. 



available to perform the three-loop calculation |^|, however, the corresponding four-loop 



vacuum diagrams are currently again out of range. 



A tempting approach for the computation of 7m.,3, which actually was considered in W5 



see also fS^), is the following: due to the properties of the IRR p6| one can set to zero 



the external momentum and choose an arbitrary subset of the internal lines to have a 
non-zero mass. A clever choice is to allow only the quark propagator which is attached to 
the left vertex a non-zero mass M. This has the advantage that the /-loop integrals can 
be solved by iterating a massive one-loop vacuum and a (/ — l)-loop massless propagator- 
type integral where the external momentum of the latter exactly corresponds to the loop 
momentum of the former. Thus, even at four-loop order at most massless three-loop 
two-point integrals have to be evaluated. It is sufficient to compute their finite parts, as 
only the 1/e poles of the /-loop integral are needed. However, there are two subtleties 
connected to this choice of IR structure. The first one is connected to the asymmetry 
one introduces due to the choice of the massive line. As a consequence the "left" vertex 
has to be renormalized differently from the "right" one. This we will explicitly see in the 
formulae we derive below. 

The second disadvantage is the occurrence of IR divergences. At this point the idea 
is to use the R* operation in order to subtract them. However, an effective evaluation of 
the roughly 6000 diagrams is only possible if the computation of the IR renormalization 



constant can be performed in a global way. This was achieved in Ref. |15] and will be 



13 



described in the following. 

As the explicit formulae for Z2 and are not yet available in the literature and 
as they will be published elsewhere [37] we present the derivation of a formula for the 
renormalization constant of the vector current correlator, Z^™ ^ which is quite similar to 
the one of Z^ and Zm- Z^™ evaluated at four- loop order immediately leads to corrections 
of order a\ to the cross section o[e^e^ —>■ hadrons) [p^ . 

Our starting point for the computation of Z^'^ is the renormalized vector current 
polarization function^ which can be cast in the form 



(2.18) 



Il^''^'^{q) is the bare correlator as indicated by the index "0". Z2 corresponds to the wave 
function renormalization constant and Zy is the renormalization constant for the vector 
current j"" = ip'y'^ip defined through 



(/ 



^2 



[2.19) 



bare 



Note that Z'^™ as defined in Eq. 



T^ ) receives contributions starting from one-loop order. 

we use for 11°'^'^ (g) a Dyson- 



In order to derive a formula for the computation of Z"^ 
Schwinger-type representation containing the full fermion propagator, G", and the proper 
photon-quark vertex function, F'^''^, leaving out one integration over the final loop mo- 
mentum which we call p. The resulting diagrammatic representation is visualized in 
Fig. p^(a). After contraction with g^^, one obtains 



^1 Tr 

Z2 



p 



i2n) 
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z^,G%p,a',)T'''^{p,q,a',)G'ip + q 



1-D) 



(2.20) 



where the minus sign in front of the first term on the right hand side accounts for the 
closed fermion loop corresponding to the p integration. In Eq. (|2.20|) all colour indices 
have been suppressed and the dependence of G° and F^''^ on the bare coupling a° is made 
explicit . 

Next we exploit the finiteness of 11^ in Eq. ( p.l8| ) and the fact that Z'^™ only contains 
poles in 1/e. This leads to an explicit formula for Z'^^ which reads 




d^p 



D 



«7m 



G°(p,a°)F°'^(p,g,a°)G°(p + g, 



a. 



(2.21) 



In this equation we consider the limit g — on the right-hand side. This is possible as 
in the MS scheme Z^™ does not depend on any dimensional scale. In order to evaluate 

For a precise definition see Eq. (Islsl) . 
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(a) 



(b) 



Figure 2.3: (a) Graphical representation of Eq. ( p.20| ). In the case of the fermion propa- 
gator the analogous diagram looks like the one in (b). In both cases an artificial mass is 
introduced in the propagator of the fermion line attached to the left vertex. 



the four-loop contribution to the polarization function and T^''^ have to be inserted at 
tree-level and at one-, two- and three-loop order. 

Let us next have a closer look to the renormalization constant Zy. For the "right" 
vertex we have Zy = Z^- However, due to the artificial mass which we want to introduce 
in the fermion line connected to the "left" vertex the equation Zy = Z^ can not be applied 
immediately. Instead we write for the "left" vertex Zy = 1 + bZy and thus obtain from 

Eq. [ wm 
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(2.22) 



In the first line the presence of the artificial mass is made explicit in the factor p^/ (p^— M^) 
which effectively replaces the massless propagator by a massive one. On the other hand, 
in the second equation M can be set to zero, as SZy already contains the effect of M. 
In a next step we want to perform the limit g — 0. In particular, we use ^ and 
apply the hard- mass procedure to the first term in Eq. ( p. 22 ) which transforms it to 
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(2.23) 



where --^^SV 0) represents the sum of all one-particle-irreducible (IPI) hard subgraphs 
containing the mass M. It gets contributions starting from one-loop order. In Eq. ( |2.23|) 
we introduced the R operation in order to treat the IR divergences of those terms which 
lead to massless tadpoles for g — 0. Note that in general also the first term in ( |2.23| ) 
contains IR divergences. However, they originate from the hard-mass procedure and are 
correlated to the UV poles of the second term in Eq. ( |2.23| ) (cf. the discussion around 
Eq. (|TlD). 

The IR divergences which occur in the last two terms of Eq. ( p. 231 ) for g = are taken 
care by introducing appropriate renormalization constants. They can be determined in 
a global manner by setting g = in Eq. (|2.21|) and applying the R operation. As the 
left-hand side is unaffected we have 
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The application of R generates only one term, namely the IR counterterm, as we set 
massless tadpoles to zero. This counterterm is a global one as it treats the IR divergences 
of the sum of all diagrams. We want to mention that this convenient aspect is new as 
compared to older calculations where the R operator has been used. In Ref. [^, e.g., 
the R operator has been applied to each diagram individually which in practice is quite 
tedious. 

Inserting Eq. (|]2|) into (^^231) finally leads to 
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(2.25) 
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Note that at this point the relation SZy = Z2 — I can be used. 

At one- loop order only the first line of Eq. (|2.25| ) contributes where all renormalization 
constants can be set to one. Furthermore the functions and T^'^ take their tree-level 
values and only the one-loop diagram shown in Fig. |2.4| has to be evaluated for external 
momentum zero and finite quark mass in one of the fermion lines. In this case no IR 
divergences occur. 




Figure 2.4: One- and two-loop diagrams contributing the vector current correlator. 



A non-trivial contribution from the terms in the second line of Eq. ( p.25| ) occurs for the 
first time at two- loop order. In particular, for STf^^ the one-loop photon-fermion vertex 
corrections have to be evaluated where one of the quarks carries mass M and the other 
one is massless. All external momenta equal to zero. In the three two-loop diagrams (see 
Fig. |2.4| ) contributing to the first line we consequently assign a mass M to the fermion 
(not the anti-fermion) line attached to the "left" vertex and set the masses in all other 
propagators to zero. For vanishing external momentum in each diagram a massless one- 
loop two-point function can be identified which can be solved with the help of Eq. ([A.4|) 
in Appendix |A.2| . The external momentum coincides with the loop momentum of the 
remaining one-loop vacuum integral which is also expressible in terms of F functions (cf. 
Eq. O). 

In general, for a /-loop calculation the renormalization constants occuring in Eq. (|2.25| ) 
are only needed at loop order (/ — 1). The loop integrals which are necessary to compute 
the first expression in Eq. ( p^.25| ) are massless (Z — l)-loop two-point functions and one-loop 
vacuum integrals. For the evaluation of STjj^, {I — l)-loop vacuum integral are needed. 



In particular, for I = 4 the occuring integrals are very well studied (cf. Refs. 



In the above equations the bare coupling constant is used as a parameter. It has 
to be expressed through the renormalized version using the (/ — l)-loop formula for the 
renormalization constant Zg (cf. Eq. ( |3.3|) ). We want to note that in our case the mass 
M needs not to be renormalized as in the diagrams no subdivergence is present which 
could induce a corresponding counterterm. 

Besides the application to the vector current correlator, which leads to corrections 
of order to R{s) |3^, the method described in this Subsection has successfully been 
applied to the four-loop fermion propagator and the correlator of scalar currents in order 
to evaluate the the four-loop contribution to 7^ [0] and corrections of order to the 
decay of a scalar Higgs boson respectively. 

Due to the large number of genuine four-loop diagrams an automation of the compu- 
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tation is mandatory. The four-loop calculations of [T^ have been performed with 

the help of the package GEFICOM [^. Within this framework the computation of STjj^ 
is straightforward as one- two- and three-loop vacuum integrals are directly accessible. 
For the four- loop contributions, like in the first line of Eq. ( p.25|) , some tricks are neces- 



sary. In particular completely massless four-loop diagrams are generated. In a next step 
a Mathematica [^] program is used to identify the "left" vertex and to introduce the 
mass M. At this point the topology of the massless three-loop diagram is fixed and the 



mapping of the momenta according to the notation of MINCER can be performed. One 
has to exploit that the whole diagram is logarithmically divergent and the mass dimension 
is given by the artificial mass M. Then the massless integration and at the end also the 
one-loop massive one can be performed. 

2.3 Massive vacuum integrals and the /3 function 



A different approach has been employed in |T6[ in order to compute the four-loop 
contribution to the (3 and 7^ functions. Also here the basic property of the IRR has 
been exploited and the integrals have been simplified by modifying the IR behaviour of 
the Feynman diagrams. This time, however, a different attitude has been adopted than 
in the previous section: in all denominators, also the ones of the gluons and ghosts, a 
common mass parameter, M, has been introduced (see also This avoids completely 

the IR divergences. On the other side, however, both gauge invariance and useful Ward 
identities are broken by this method. Moreover, multiplicative renormalization, which is 
very convenient in practical renormalization is lost. Nevertheless the R operation applied 
to the individual diagrams still works and in principle it can be used to compute the over- 
all renormalization constant. However, this is not at all practical, especially as in the case 



of the (3 function roughly 50 000 diagrams have to be considered |T^ and an automatic 
treatment is absolutely mandatory. This is achieved by introducing effective vertices and 
propagators which incorporate all lower-order renormalization constants |T^. Further- 



more counterterms have to be introduced which correspond to a renormalization of the 
gluon and ghost mass and an overall renormalization constant for the gluon propagator. 

In the following we want to illustrate this method by considering two-loop QCD cor- 
rections to the photon propagator and evaluate the corrections to the wave function 
renormalization defined through 

A^'f" = ^/Z^A^", (2.26) 

where is the photon field. For convenience we also introduce the renormalization 
constants Z7 and Z2 via 

eV°7^^°^J = Zje^Yi^Af,, (2.27) 
where e is the electromagnetic charge and ip is the fermion field. The corrections up to 



order aas to Z3 are obtained from the diagrams of Fig. |2.4| which contribute to the photon 
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polarization function. The latter can be written in the form 



1^. 



(2.28) 



The canonical way to compute would be to consider all particles as massless and 
evaluate the occuring massless two-point function in dimensional regularization. In this 
case St('?^) is proportional to q^ and ^^.(g^) = which is due to gauge invariance. One 
observes that the sum of the bare two-loop diagrams only contains a simple pole in e and 
that it is independent of the QCD gauge parameter, ^. Furthermore, due to the Ward 
identity connecting the photon-quark vertex to the quark self energy no renormalization 
is necessary and the bare diagrams of Fig. ^]4| already lead to the final answer where the 
formula 



^3 



1 



(2.29) 



can be used. 

In the method proposed in [TB, H3| a common mass M is introduced in all lines 
for the computation of the renormalization group functions. As a consequence gauge 
invariance is broken and Si(g^) is not zero any more. Furthermore, both Er(g^) and 
Ei(g^) contain terms proportional to M^. One also observes that the sum of the bare 
diagrams contains a pole of the form ^/e. Clearly the R operation applied to the indi- 
vidual diagrams would still lead to the correct answer. However, we want to "immitate" 
multiplicative renormalization as close as possible and introduce effective vertices and 
propagators. In the case of the photon-quark vertex we write 



r 



Yz^ = Yi^ + szl 



(2.30) 



which actually corresponds to multiplicative renormalization. The quark propagator is 
modified to 
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(2.31) 



Note that for M = the terms proportional to 8Zl and bZl are proportional to the 
Born diagram and a cancellation takes place due to the Ward identity Zl = Zl- For 
M 7^ 0, Z7 and Z^ have to be determined form the vertex correction and quark self- 



energy, respectively, 
namely 



To our approximation they take the same values as for M = 0, 



Zl 



Zl 



a,CV/l 1 
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(2.32) 



where ^ is the gauge parameter defined in Eq. 

The use of Eqs. (|2.3CI| ) and (p.31|) for the computation of the one-loop diagram of 
Fig. 2A induces terms of order aUg which render the pole part independent of ^. However, 
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U'^" still contains terms proportional to g^'^M'^ . They are removed by introducing a local 
mass counterterm for the photon of the form M'^A^A'^. Also for the gluon field a similar 
counterterm has to be introduced. It is contained in the effective gluon propagator as it 
is needed for the cancellation of sub divergences. 
Finally one arrives at 



TT \3e TT 8e 



[2.33) 



In Ref. |T^ this method has been applied to obtain the four-loop contribution of the 
QCD P function. It has been applied to the ghost-gluon vertex, the gluon propagator and 
the ghost propagator in order to obtain the corresponding renormalization constants Zi, 
Z3 and Z3 (cf. Eq. (p. 3D) and finally the one for via 
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The (3 function is defined through 
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where is the number of active flavours. For completeness we want to list the results 
for the coefficients which are given by ||^, |n| 
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(2.36) 



C is Riemann's zeta function, with values C2 = tt'^/Q and ~ 1.202 057. 

One of the main new achivements of Ref. [1^ is the treatment of the four-loop vacuum 
diagrams. The task is simplified due to the fact that only the divergent parts in e are 



needed. The method of integration-by-parts has been used to derive recurrence rela- 
tions which reduce a general four-loop integral to a linear combination of simple integrals 
and two (difficult) master integrals. 
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3 Decoupling of heavy particles 



Quantum corrections to processes involving only light degrees of freedom contain in gen- 
eral the whole particle spectrum. In particular also heavy particles with masses much 
larger than the energy scale of the considered process contribute. It is highly desirable 
that in the limit where the heavy mass, M, goes to infinity its contribution to the light- 
particle Green function must tend to zero like /i/M where /i is a typical scale of the 
process. This is exactly the content of the so-called decoupling theorem which is proven 
in 



To be more precise let us consider an example, namely the production of heavy quarks 
in e+e~ annihilation. Due to the hierarchy in the quark masses there is a clear separa- 
tion into light and heavy. For center-of-mass energies, ^/s, of about 40 GeV we are well 
below the production threshold of top quarks and thus we expect their influence to be 
suppressed by y^/Mj. Analogously the contribution of bottom quarks to the production 
of charm quarks close to the threshold must be proportional to y/s/Mi^. In our example 
the heavy quarks enter the first time at order a1 via the diagram pictured in Fig. |3.1| . 
If the contribution of this diagram is computed in a momentum-subtraction scheme one 
indeed observes this behaviour^. However, in this example, and also in most other QCD 
processes it is much more convenient to use a mass-independent renormalization prescrip- 
tion, like the MS |]T2[ or its popular modification, the MS ||13[ scheme. These schemes 
are characterized through the fact that their renormalization group functions are mass- 
independent which makes renormalization group improvements much more transparent. 
It also has the advantage that the computation of the renormalization group functions 
themselves is significantly simplified (cf. Section ^). On the other hand there is the big 
drawback that the decouphng theorem of Appelquist and Carazzone |^ does not hold in 
mass-independent schemes. This is due to the mass-independence of the renormalization 
group functions which implies that, e.g. in case of QCD, the top quark and the down 
quark have identical contributions. 

Note that in the broad classes of momentum subtraction schemes the decoupling the- 
orem is valid. However, the calculations are much more complicated and in general a 
coupled system of differential equations involving also the quark masses and the gauge 
parameter has to be solved in order to obtain the running of the couplings. 

Coming back to our example this means that the cross section a{e^e~ ^ bb + gluons) 
evaluated in the MS scheme with y/s = 40 GeV does not behave like y/s/Mt but still 
contains logarithms of the form ln(s/M^^), which arise at order from the diagram in 
Fig. 

Clearly, both from the theoretical and practical point of view this is not acceptable. 
The way out is the explicit construction of an effective Lagrangian where the heavy 
particles are integrated out. This means that the dynamical degree of freedom of the 
heavy quark is removed, which manifestly leads to power-suppressed contributions of the 
latter. This will be performed in Section p]T] for the case of QCD. The effective Lagrangian 



^Actually the proof of the decouphng theorem in |4q] is performed for a momentum-subtraction scheme. 
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Figure 3.1: Feynman diagram contribution to the process e'^e bb which gives rise to 
top-quark-dependent terms. 



only depends on hght degrees of freedom where the couphngs0 are multiphcatively re- 
scaled by the so-called decoupling constants. These constants are universal and have 
been computed in the case of QCD up to order |]49 |. 



Concerning the practical consequences we again want to consider the example of e~^e~ 
annihilation. For a calculation at order one computes all relevant diagrams including 
the one of Fig. |^T^. As mentioned above the result diverges proportional to the logarithm 
of the heavy quark mass. Now one has to remember that the coupling]] has to be changed 
according to the known rules which describe the transition to the effective Lagrangian. 
Thus one arrives at a physical observable expressed in terms of parameters of the effective 
Lagrangian and it can be explicitly checked that the dependence on the heavy quark mass, 
Mt, is power-suppressed — in the case at hand it goes like s/M^. 

Pioneering work in the computation of the decoupling constants has been done in . 



In Ref. 1 51] the decoupling constant for has been computed at the two- loop order. The 



crucial idea of the method is based on the fact that the decoupling theorem |^ works in 



momentum subtraction schemes. Thus after relating the corresponding coupling constant 
to as defined in the MS scheme both in the full and the effective theory, it is possible to 
derive differential equations for the decoupling constant. They can easily be solved. For 
the corresponding integration constant a two-loop calculation is needed. 

Thirteen years later the authors of evaluated the corrections of order to the 
total decay rate of the Z boson induced by a heavy top quark. In order to make the 
decoupling explicit the two-loop result of is needed. However, it turned out that after 
expressing the decay rate in terms of effective parameters the top quark did not decouple. 



^Here we mean "normal" coupling constants, but also masses, gauge parameters, etc.. 

^Note, that, in the case where the light quark masses are neglected, only as remains as a parameter. 
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Thus in a second evaluation of the decouphng constant for as has been performed 

The result of was confirmed in . In the meantime, 



with a different result as in 151 



the authors of Ref . [ 51 1 have revised |^ | their original analysis and have found agreement 
with Ref. 

The method of Ref. for the computation of the decoupling relations is based on the 
evaluation of a top quark contribution to a physical quantity. In particular the authors 
of |52| considered corrections to the massless quark propagator with an additional zero 



momentum operator insertion. In order to obtain the decoupling relations to order a 
three-loop calculation corresponding to the order corrections is necessary. Correspond- 
ingly, the decoupling relation to order would require a four-loop calculation. In our 
method, which is described below, only the computation of three-loop vacuum diagrams 
are necessary. 

3.1 Construction of an effective Lagrangian in QCD 

The main idea of effective theories is that the dynamics at low energies does not depend 
on the details of the dynamics of high energies. Thus the low-energy physics can be 
described using an effective Lagrangian which does not depend on the additional degrees 
of freedom present at high energies. The only effect of the high-energy parameters are 
modified couplings of the effective Lagrangian with respect to the full one. 

In the following we want to describe the construction of the effective Lagrangian in 
QCD with one heavy quark of mass rrih and rii light quarks. 

Our starting point is the full QCD Lagrangian which reads 



£QCD 
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2(1-0 



/=i 
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9 - mf) ipf 



abc J} 



(3.1) 



where the field strength tensor is defined through G"''^'^ = d'^G'"''' -d''G^'^+gsr^''G^'''G''''' . 
/"'"^ are the structure constants of the QCD gauge group, gs = y/4:7ias is the QCD gauge 
coupling, ipf is a quark field with mass ruf, G"''^ is the gluon field, c° is the Faddeev-Popov- 
ghost field, and n/ = rii + l is the total number of quark flavours. -D^ = d^ — igsiX"" /2)G1^ 
is the covariant derivative in the fundamental representation and A° are the Gell-Mann 
matrices. For convenience we list the gluon propagator resulting from (p.lD 



DM 



-9 



(3.2) 



In this convention ^ = corresponds to Feynman gauge and = 1 to Landau gauge. 

For later use we define the renormalization constants connecting the bare and renor- 
mahzed quantities in Eq. (|3.1|): 



9s 



fJ'^Zggs , 



(3.3) 
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In addition one introduces the renormalization constants of the quark-gluon, three-gluon, 
four-gluon and gluon-ghost vertex which are denoted by Zip, Zi, Z^ and Z^. The Slavnov- 
Taylor identities connecting the different renormahzation constants can, e.g., be found 



m 



It is clear that in the effective Lagrangian C^^^ all explicit trace to the heavy quark 
must have disappeared. However, the mathematical structure must be identical to the 
one of C^'^^ in Eq. (|3.1| ). This is because represents the most general Lagrangian 

describing the interaction of quarks and gluons and respecting the symmetry properties 



imposed by the Becchi-Rouet-Stora-Tyutin invariance ||5J]. Of course, this also has to 
be respected by a nr flavour theory described by C!^^^ ■ Nevertheless the parameters and 
fields of L^^^ are different from the ones of the full theory. It is convenient to define 
the corresponding relations in analogy to the renormalization constants of Eq. ( |3.3| ) and 
introduce multiplicative factors — the so-called decoupling constants Q. Thereby it is 
advantageous to consider the decoupling relations in the bare theory which is indicated 
by the superscript zero. The renormalization is performed afterwards. Thus we define 

9V = O's , rn'/ = , - 1 = cue - 1) , 

(3.4) 

where the primes mark the quantities of the effective n^-flavour theory. 

Taking into account these considerations we can write down a defining equation for 
the bare effective Lagrangian in terms of the full Lagrangian with re-scaled parameters: 

^Tial K GJ- C°) = C^^'^ig'/, m'/, f''; Gj'-, c°''-) , (3.5) 

where /^Q*^^ is given in Eq. ( p.lD , q represents the ni light-quark flavours and (f collectively 
denotes all bare decoupling constants of Eq. ( |3.4|) . Once they are explicitly computed the 
effective Lagrangian is completely determined. Green functions of light fields obtained 
from agree with the ones of C^q^ up to terms suppressed by inverse powers of the 

heavy quark mass. 

In the language of effective theories the computation of the decoupling constants is 
referred to as matching calculation. It can be performed in a more or less complicated 
way. As it is even nowadays highly non-trivial to apply the methods of |^ and 



at order a^, we developed a procedure which relates the n-loop decoupling constants of 



Eq. ( |3.4j) to n-loop massive one-scale integrals. It will be described in the next Subsection. 



3.2 Computation of the decoupling constants 



In order to compute the decoupling constants of Eq. (|3.4]) we have to find convenient 
Green functions which have to be considered both in the effective and full theory. For 
this reason we define the bare two-point functions for quarks, gluons and ghosts as follows: 



1 



m — ^ — S°(p) 



d^a;e'^^(T<(x)V5°(0)) 
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-9i^u + 



-p2(i + n°G(p2)) 



+ ... = i / d^a;e'^'"(TG'°'"(x)G°'^(0)) 

) 



-pMl + no(p2)) 



I I d^xe'P''{Tc'''%x)c'''\0)) . (3.6) 



The ellipses in the case of the gluon propagator indicate the longitudinal part which we 
are not interested in. 

From the two-point functions of Eqs. ( |3.6| ) we will be able to obtain Cm' C3 ^^^^ C3 • 
In order to get a relation involving one has at least to consider three-point functions 
where the coupling Qs already appears at Born level. As the vertex between the gluon and 
ghost is the least complex one we will take it for the computation. In amputated form it 
is defined through 

P'g's {-^r'' [1 + ^G-cciP, k)]} + ... = z'J dxdye^^^-^+'-y'^ {Tc''%x)(^'\0)G''^'^{y)y''' , 

(3.7) 

where p and k are the outgoing four-momenta of c and G, respectively. The ellipses 
indicate other colour structures we are not interested in. Note that we pull out a factor 

on the left-hand side of Eq. { \i.7\) as it is already present at Born level. 

Let us start with the decoupling constant for the gluon field. It is obvious that the 
gluon propagator constitutes a good candidate to compute Cs- Up to terms of order l/ruh 
we have the following chain of equations 



5"' {-9^.. + 



d^a;e'P"(TG'°''^(x)G'°'''(0)) 



I I d^a;e*f^(TG°''''^(x)G°'''''(0)) 



I 

1 



(3. 



where in the second step Eqs. (|3l^ ) has been used. Note that 11^ (p^) only contains light 
degrees of freedom whereas IIq{p^) also receives virtual contributions from the heavy 
quark h. Eq. (|3.8|) provides a formulae for ^3 

_ 1 + uUp') 

~ i + u'^{p^)- ^^-^^ 

From the construction of the effective theory it is clear that the decoupling constants do 
not depend on the momentum transfer. Thus also the right-hand side of ( |3l^ ) has to be 
independent of the external momentum p. This means that we can choose any convenient 
momentum for the computation. In particular it is possible to choose p = 0. This has 
the advantage that only vacuum diagrams have to be considered. Since we work in the 
framework of dimensional regularization all scaleless integrals can be set to zero. As a 
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consequence we have n^(0) = and the contribution to 11^ is given by the diagrams 
containing at least one heavy quark hne. In the following we attach to the corresponding 
contributions an additional index "h" and refer to it as the "hard part" . Finally we arrive 
at the compact formula 

= i + n°/(o). (3.10) 

The n-loop contribution to is related to the n-loop vacuum diagrams where the scale 
is given by the mass of the heavy quark. At one-loop order only one diagram contributes 
to C3. At two loops there are already three diagrams and at three-loop order altogether 



189 diagrams have to be taken into account. A typical example is pictured in Fig. |373 . 
We want to mention that the bare decoupling constants may contain non-local terms like 
ln(/i^/m|)/e, which is in contrast to the renormalization constants. 

After the computation of the bare diagrams the parameters (in our case as, ^ and the 
heavy quark mass, m^) have to be expressed in terms of their renormalized counterparts. 
The finite decoupling constant is obtained from Eq. ( ^.4] ) after expressing the bare fields 
in terms of the renormalized ones via Eq. (|3] 



^'^'g, = JCsG,. (3.11) 



Note that Z'^ depends on a'^ and They have to be transformed to as and ^ with the 
help of (g and Cs, respectively, which are needed up to (/ — l)-loop accuracy if Eq. (|3.11|) 
is considered at /-loop order. 

The ghost propagator can be treated in complete analogy and one arrives at the 
following formula for the bare decoupling constant for the ghost field 

4° = l + n°''^(0). (3.12) 

There is no diagram which contributes at one-loop order and one at order a^. Also at 
three-loop order the number of diagrams is moderate and amounts to 25. 
The renormalized version of Eq. (|3.12|) reads 

Cs = (3.13) 

In order to obtain expressions for (2 and Cm one considers the light-quark propagator 
which leads to the following chain of relations: 



m- f- T.{j)) mil- S^(p2)) - ^{1 + S^(p2)) 

1 



(3.14) 
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Nullifying the external momentum p in the self energies leads to 

C° = l + S°/(0), (3.15) 



s°/(o) 

S°/(0) 

and the finite expressions are obtained from^ 



c = : . xz ^ (3-16) 



-0 



Z2C2 
I 
2 

Cm = (3.18 



C2 = (3.17) 



Similarly to 11°''^ (0) there are no one- loop diagrams contributing to Sy'*(0) and S^''*(0) 
and at two- and three-loop order again one and 25 diagrams, respectively, have to be 
considered. Typical specimen are depicted in Fig. Actually, through three loops, 
Eq. ( |3.18| ) simplifies to = 1 — Sy'*(0) — S°''^(0). It should be noted that the vector and 
scalar parts separately still depend on the QCD gauge parameter ^, but ^ drops out in 
their sum, which is a useful check for our calculation. 




Figure 3.2: Typical three-loop diagrams pertinent to Sy (0) and (0). Solid, bold- 
faced, and loopy lines represent massless quarks q, heavy quarks h, and gluons G, respec- 
tively. 



The derivation of a formula for (g is slightly more involved. As a starting point we 
choose the full, i.e. non-amputated gluon-ghost Green function and get (using Eqs. (|3.4|)) 
the following equations 



9sG'Gsc\P,k) = ^'y'dxdye*(^""+'=-J^)(rc°'"(x)c°''(0)G°''='^(l/) 

f f dxd?/e*(P-^+'=-J') (rc°'''(x)c°'^(0)G°'"'^(l/) 



C3 V C3 



AO /AO 
S3 V S3 



^Note that the same symbol is also used for Riemann's zeta function C2 = 7r^/6. However, as they 
appear in a completely different context confusion is not possible. 
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So 



S3' 



(p, k) 



(3.19) 



where p and /c are the outgoing four-momenta of c and G, respectively. At this point we 
should amputate the Green functions by multiplying with the inverse propagators of the 
external gluon and ghost fields in the full theory. From the decoupling relations derived 
above we get an additional factor (C3)^C3 the right-hand side of Eq. (|3.19| ) which leads 
to 



(3.20) 



Here, Tq^^{p, k) is defined through the IPI part of the amputated gluon-ghost Green 
function as introduced in Eq. (|3.7| ). Setting to zero the external momenta we obtain 



/■o //•o 

S3 V ^3 



where 



i + r°GL(o,o) 



(3.21) 



(3.22) 



is the decoupling constant belonging to the gluon-ghost vertex. Again, the renormal- 
ized version of the decoupling constant is obtained with the help of the renormalization 
constants in the full and effective theory as 



(3.23) 



Thus, in order to compute (g one has to evaluate the decoupling constant for the gluon 
propagator, the ghost propagator and the gluon-ghost vertex. We could have chosen 
also another vertex involving the strong coupling, e.g. the quark-gluon or the three- 
gluon vertex, and would have arrived at a similar expression as in Eq. ( |3.21| ). This is in 
complete analogy to the renormalization constants where due to Slavnov- Taylor identities 
the various renormalization constants are related to each other (see, e.g., Ref. P5[). 

We should mention that r^g^(0, 0) has no one-loop contribution and there are five 
diagrams at two loops, which, however, add up to zero. A non-zero contribution is 
obtained at order where 228 diagrams contribute. A typical representative is shown in 



Fig. 3.3 



At order the three contributions to C° are still separately independent of the gauge 
parameter ^, so that the ^ independence of their combination does not provide a mean- 
ingful check for our calculation. The situation changes at C(a^), where all three parts 
separately depend on ^ and only their proper combination is ^ independent as is required 
for a physical quantity. In the calculation this has been used as a check. 

In Eqs. (|3T0| ), ( CT) , (|315|) , (|3A6D and (^2^) the bare decoupling constants C3, Cl 
(2, Cm Cg expressed in terms of vacuum diagrams. Thus, if the former are to 
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Figure 3.3: Typical three-loop diagrams pertinent to 11^ (0), n°''^(0), and rgg^(0,0). 
Bold-faced, loopy, and dashed lines represent heavy quarks h, gluons G, and Faddeev- 
Popov ghosts c, respectively. 



be computed at order a\ the latter need to be known at three-loop order. In the recent 
years three-loop vacuum diagrams have been studied extensively and a FORM package, 
MATAD has been written which allows for an automated computation. 



At this point we would like to make a comment on the different kind of poles which 
appear in the calculation. If we choose zero external momentum there are in general both 
UV and IR poles in the individual diagrams. However, the renormalization constants in 
Eq. ( |3.3| ) only contain UV divergences. Thus the two different kind of poles have to be 
identified in order to arrive at a finite expression for the decoupling constants. This is 
a special feature of dimensional regularization which allows to treat simultaneously UV 
and IR divergences and to set to zero scaleless integrals. More details to this context can 



be found in Ref. ||56[ and references cited therein. 



3.3 Results 

In the following we list the analytical results for the renormalized decoupling constants 
and Cm, which relate the physical parameters in the full theory to their counterparts 
in the effective theory. In Appendix ^ we also provide the (gauge parameter dependent) 
results for ^2 and Cs- 

As already mentioned above, the computation can be reduced to one-scale vacuum 
integrals (see also Appendix |A.2| ). However, the large number of diagrams requires the 
automation of the computation. Actually the computation of was one of the first 
application of the package GEFICOM (cf. Appendix A. 3 ), which combines several 
stand-alone program packages in order to automate the computation from the generation 
of the diagrams to the summation of the individual results. For the evaluation of the 
decoupling constants all diagrams have been generated automatically with the Fortran 
program QGRAF ||5^ and the integrations have been performed using the FORM package 
MATAD HI]. 

Note, that in contrast to the renormalization constants Zi in Eq. (|3.3|) , the decoupling 



constants Q also receive contributions from the finite parts of the loop integrals. Thus, at 
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0{al), we are led to evaluate three-loop tadpole integrals also retaining their finite parts. 

For illustration of the formalism derived in the previous Subsection we want to compute 
the lowest-order contribution to C2 which is of and comes from the first diagram 

in Fig. |3.2| . According to Eq. ( p.l5| ) the vector part has to be evaluated for zero external 
momentum. This is conveniently done with the help of 



p=0 



(3.24) 



where a simple Taylor expansion up to linear order has to be performed for S°''*(p). As 
only vacuum integrals are involved in the computation one can easily perform the tensor 
reduction in the numerator and then use Eq. (|A.7| ) in order to arrive at the result for the 
bare decouphng constant 



1 + ( ^ ) CpT 



\ In-^ + 0(6) 

16£ 96 8 ml ^ ^ 



(3.25) 



Now Eq. ( p.l7| ) can be used to get a finite result. Therefore the ratio Z2/Z2 is needed to 
order where the parameters a'^ and C,' in Z2 have to be expressed in terms of and 
^. with the help of Eqs. ( |3.4|) . Actually, to our order the contributions from (g and 
exactly cancel and we are left with 



1 + 1^ 

TT 



2 CfT 



16e 



(3.26) 



Inserting Eqs. and (|]2|) in ( CT) finally leads to 



/-MS 
S2 



1 + 





" 5 


- In — 5- 








96 ~ 


8 mi_ 



(3.27) 



which agrees with Eq. (|B.5[) . 

In the same way also the scalar part of the quark self energy can be treated in order 
to obain with the help of Eq. ( |3.16 ) a result for (m- Taking also the three- loop diagrams 
into account we obtain 



-OS 



^432 36 Ml 12 M| 



407 5 1 D / 121 
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where Ca ~ 1.202 057 and 



16Li4 
-1.762 800 



2^-fC4-4C2ln^2 + ^ln^2 



(3.29) 



with Li4 being the quadrilogarithm, is a constant typical for three-loop vacuum dia- 
grams |5^. rii = Uf — 1 is the number of light-quark flavours, and Mh is the on-shell mass 
of the heavy quark h. For the numerical evaluation in the last line of Eq. (|3.28| ), we have 
chosen = M^. The 0{al) term of Eq. (^TIEj ) is computed in Ref. and the C(a^) 
term can be found in IEQI. 



The proper combination of the three ingredients entering the calculation of (g, namely 
the hard parts of the gluon and ghost propagators and the gluon-ghost vertex correction, 
lead to 




24 "^Ml + Se ^ Ml 



521 , u 
m 



131 , 9 u 
In"' 




, 3 /i- / 2479 C2 
216 Ml V31104 9 1728 



TT 



+ (-5.3239 + 0.2625 72;) 




Ml 



(3.30) 



The 0{al) term in Eq. ( 3.30 ) is published |49|. Leaving aside this term, the results in 
Eq. ( piuD can be found in Refs. 0, g. 

Notice that the 0{al) terms of Cm and (g depend on the number ni of light (massless) 
quark flavours. However, this dependence is feeble. 

The generalization of Eqs. ( p.28| ) and ( p.30| ), appropriate for the gauge group SU(Ai'c), 
is listed in Appendix ^ There, also the results expressed in terms of the MS quark mass 
can be found. In Ref. |^ the leading Yukawa corrections proportional to GpM^ have 



been computed for Cm and C. In Appendix also these results are listed. 



3.4 Applications 

Cross section cr(e+e~ bb) 

At this point we would like to pick up the example mentioned at the beginning of this 
section, namely the 0{a1) corrections to the cross section e~^e~ bb. If we consider 
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center-of-mass energies where the first five quarks can be neglected one obtains^ 



R{s) 




176 



675 135 



In 



(3.31) 



where Nc = 3 has been chosen and at order only the contribution from the diagram in 
Fig. p.l| is displayed explicitly. All other diagrams are summed in the constant C2 [0 . In 
Eq. ( ^.31| ) the definition of the coupling constant still includes the top quark as indicated 
by the superscript "(6)". Otherwise the computation of the diagram in Fig. p.l| would 
hardly be possible. One recognizes that for Mf oo the contribution from the top quark 
raises logarithmically. Note that the choice /i = Mt does not help as it introduces In Mt 
terms in C2- At this point Eqs. (|3.4|) and ( p.30|) can be used to replace ) in favour of 



af'^ which leads to 



R{s) 



a 



(5) 



TT 



TT 



C2 + 



AM? 



135 Mi 




0{al 



(3.32) 



Now the top quark is decoupled, i.e. its contribution goes to zero for ^ oo as R{s) is 
expressed in terms of af^ which is the parameter of the effective theory. 



Determination of af\Mz) and irS^\Mz) from measurements at the r mass scale 

In the previous example the decoupling relation was only needed to one-loop order. How- 
ever, the three-loop terms will be indispensable in order to relate the QCD predictions 
for different observables at next-to-next-to-next-to-leading order. Meaningful estimates 
of such corrections already exist 0, |62| . 

Once the corrections of order for R{s) are known they can be used to determine 
^^\Mz) from the knowledge of a'f\Mr) and compare it with other measurements. One 



would use R{s) to extract^ af\M^) with an accuracy of order from the data. Then 
one would use the four-loop (3 function |[T^ in order to perform the running to the bottom- 
quark threshold. There the three-loop matching relations would be necessary for a consis- 
tent decoupling. Using again four-loop running finally leads to af\Mz)- In the following 
we will illustrate this procedure. However, instead of determining a^^^Mr) via R{s) we 
directly assume a value for a^'^)(Mr) and evaluate af\Mz) for different choices of the 
matching scale /i*^^-*. For the three-loop evolution in connection with two- loop matching 



^Note that this result can immediately be obtained from the example considered in Appendix A.l 
Eq. (O). 



""The described procedure can also be applied to o's^Mt). Only for simplicity and transparency we 



have chosen a"^\Mr) 
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this has been done in Ref. |Q. We are in a position to explore the situation at the next 
order. It is instructive to include in the analysis also the tree-level and one-loop matching. 

Going to higher orders, one expects, on general grounds, that the relation between 
ai"^ and where /i' <C fj.^^^^ -C /i, becomes insensitive to the choice of 

the matching scale, fj.^'^f\ as long as fx^^'f^ = 0{mh). In the above-mentioned situation 
we consider the crossing of the bottom-quark threshold. In particular, we study how 
the n^^^ dependence of the relation between a^J^\Mr) and a^^\Mz) is reduced as we 
implement four-loop evolution with three-loop matching. Our procedure is as follows. 
We first calculate a^^\ii^^^) by exactly integrating Eq. ( p.35| ) with the initial condition 
a(4)(M^) = 0.36, then obtain ^^^^(/i^^)) from Eqs. (|3j) and (CT ) with = 4.7 GeV, 
and finally compute af\Mz) with Eq. ( |2.35| ). For consistency, A^-loop evolution must be 
accompanied by (A^ — l)-loop matching, i.e. if we omit terms of 0{a^~^'^) on the right- 
hand side of Eq. (|2.35| ), we need to discard those of 0{a^) in Eq. ( |3.30| ) at the same 
time. In Fig. |3.4| , the variation of af'^Mz) with jjL^^'> /Mh is displayed for the various levels 
of accuracy, ranging from one-loop to four-loop evolution. For illustration, /i*^^) is varied 
rather extremely, by almost two orders of magnitude. While the leading-order result 
exhibits a strong logarithmic behaviour, the analysis is gradually getting more stable as 
we go to higher orders. The four-loop curve is almost flat for /i*^^) ^ 1 GeV. Besides the 
/i*^^) dependence of af\Mz)i also its absolute normalization is significantly affected by the 
higher orders. At the central matching scale /i^^-* = M^, we encounter a rapid, monotonic 
convergence behaviour. 

Similar analyses may be performed for the light-quark masses as well. For illustration, 
let us investigate how the /i*^^-' dependence of the relation between = 'm^^HfJ'c) and 
■m^^\Mz) changes under the inclusion of higher orders in evolution and matching. As 
typical input parameters, we choose /ic = 1.2 GeV, Mf, = 4.7 GeV, and a^^\Mz) = 0.118. 
We first evolve m^^\n) from jj, = jjLc to jj, = fi^^^ with the help of Eq. (|2.12| ), then obtain 
m^^^(yU*^^^) from Eqs. (|3.4|) and ( 3.28| ), and finally evolve m^^^(yu) from /i = /i'-^^ to /i = Mz 
using Eq. ( p.l2| ). In all steps, a;i"'^^(/i) is evaluated with the same values oi nj and /i as 

mc"^^(/i). In Fig. |3]^, we show the resulting values of m'^\Mz) corresponding to A^-loop 
evolution with (A^ — l)-loop matching for A^ = 1, . . . , 4. Similarly to Fig. p.4| , we observe a 
rapid, monotonic convergence behaviour at the central matching scale fx^^^ = Mi,. Again, 
the prediction for A^ = 4 is remarkably stable under the variation of fi^^^ as long as 
^ 1 GeV. 

The various formulae describing the running and the decoupling of as and the quark 



masses are implemented in the program package RunDec |6J]. It is realized in Mathematica 
and provides a convenient possibility to check, e.g., the figures of this Subsection and 
eventually update the numerical input values. In particular. Fig. can easily be re- 
produced with the help of the procedure A1L2A1H [] . After loading RunDec the command 
A1L2A1H [0 . 36 , 1 . 777 , {{5 , 4 . 7 , mu5}} , 91 . 187 , 1] provides the result for a^p (Mz) where 
the matching has been performed at the scale mu5. / = 1, 2, 3, 4 corresponds to the num- 
ber of loops used for the evolution and the values a^'^^M^-) = 0.36, M^. = 1.777 GeV, 
Mh = 4.7 GeV and Mz = 91.187 GeV have been adopted. 
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af(M,) 




10 



Figure 3.4: /i^^^ dependence of af'^Mz) calculated from a^J^^M^-) = 0.36 and = 
4.7 GeV using Eq. ( p.35| ) at one (dotted), two (dashed), three (dot-dashed), and four 
(solid) loops in connection with Eqs. ( |3^ ) and ( p.30| ) at the respective orders. 



4 Scalar dimension-four operators in QCD 

In the last section we constructed an effective QCD Lagrangian which results from in- 
tegrating out a heavy quark. From the knowledge of its structure we determined the 
coefficient functions (i.e. the decoupling constants) by computing bare Green functions 
in the full and effective theory. 

A different point of view is adopted for the construction of the so-called non-relativistic 
QCD (NRQCD) ||6^. It has been developed in the context of heavy quarkonium physics 



in order to separate the relativistic scales associated with the mass of the heavy quark, 
M, from the non-relativistic ones which are of the order Mv, where v is the velocity of 
the quark. As a result one obtains an effective Lagrangian which is ordered in inverse 
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mf^ (M J 




(5) 



Figure 3.5: fx^^^ dependence of ml^\Mz) calculated from = ml^\nc) = 1-2 GeV, 
Mfe = 4.7 GeV, and a';^\Mz) = 0.118 using Eq. (|2T^) at one (dotted), two (dashed), 
three (dot-dashed), and four (solid) loops in connection with Eqs. ( p.4| ) and ( p.28| ) at the 
respective orders. 



powers of M. More recently a further step has been undertaken and potential NRQCD 
(pNRQCD) has been introduced in order to account for the separation of the scales Mv 
and Mf^. For recent overviews we refer to [l66i |67[] . 

In this section a different approach is considered. It is based on Wilson's operator 
product expansion (OPE) |68| which is a powerful method for the construction of effective 
theories. The basic idea is to introduce local operators Oi of appropriate dimension. They 
are formed by the light degrees of freedom and parameterize the long-distance behaviour. 
The operators are accompanied by coefficient functions which contain the remnant of the 
large parameters of the theory. 

As a simple example of the OPE one can consider the decay amplitude of the muon 
in the SM. Due to the fact that the momentum transfer is much smaller than the mass 
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of the W boson, Mw, the former can be neglected with respect to the latter. This leads 
to the famous four-fermion interaction which is generated from dimension-six operators. 
Thus, for dimensional reasons the coefficient functions contain a factor 

In QCD the scalar operators of dimension four have been studied in great detail [^, 



[70| , [7l[] . A comprehensive survey concerning their renormalization properties and mixings 
is performed in |71[]. In particular, all renormalization constants of the operators have 
been expressed in terms of renormalization constants which already appear in the QCD 
Lagrangian (cf. Eq. (|3.3| )). 

In this Section we want to discuss two applications of the scalar dimension-four oper- 
ators. In the first one we consider the decay of an intermediate-mass Higgs boson with 
M]y <, Mfj <, 2Mw into quarks or gluons. Here the top quark is considered to be heavy 
and thus it will only contribute to the coefficient functions. On the other hand, the scale 
in the operators is given by the mass of the Higgs boson. A somehow complementary 
situation is considered in the second example: the quartic quark mass corrections to the 
vector current correlators. Here the mass of the quark sets the scale in the operators. It 
is considered to be small as compared to the external momentum which manifests itself 
in the coefficient functions. 

For definiteness we want to list the operators of dimension four in this introductory part 



0,/ 

q ' 



and briefly discuss their renormalization. In |7T| the operators are classifled into gauge- 
invariant and non-gauge-invariant ones. Furthermore a distinction is made whether the 
operators vanish or not by virtue of the equation of motion. The gauge-invariant operators 
read|3 

= (Gjr 

O', = K')\ (4.1) 

where Gj^'" and are deflned after Eq. Note that O^ vanishes after the application 

of the equation of motion. 

In order to obtain a closed system two more operators have to be considered 

Ol = (V;:''9'^c°'''^)c°'''" , (4.2) 

where V^^ = S"'^d^ — gsf°'^'^G'^^ is the covariant derivative acting on the gluon and ghost 
flelds. The operators in Eq. ( |4.2| ) are not gauge-invariant and thus do not contribute to 
physical observables. 



^^For consistency, the operators should also have a prime as a superscript as they are built by quantities 
of the effective theory. However, we refrain from introducing this additional index. 
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From the practical point of view the operators of Eqs. 
Feynman rules which can be read off from Eqs. ( [4.1|) and (^.2\ 



yi]) and ( [4.2| ) define new 
E.g., in Section [4.1.4| we 



have to compute the correlator {OiOi) up to three- loop order. This makes it necessary to 
extract the Feynman rules for the coupling of Oi to one-, two-, three- and four gluons from 
Eq. ( |4.1| ). We refrain from listing them explicitly but consider as illustrative examples 
the coupling of (9° ^'^d to two gluons. This will be useful further below in Section [4.1| 
in order to demonstrate the evaluation of the corresponding coefficient functions. 

The coupling to two gluons is obtained from those terms of Eqs. and ( [4.2| ) con- 
taining two gluon fields. In the case of Oi and they read 



two gluons 



two eluons 



2 (9^G°'''"(9^G'^'°'''^ - df,Gl'''^d''Gf''^'''' 



0/,b 



(4.3) 



If we adopt the momenta fiow as shown in Fig. ^Tl] this results in the following Feynman 
rules for the vertices 



iS^b [-^/^'^ (p^ . + P2 ■ P2) + P2P2 + PlPl] 



(4.4) 



^(r(r(r(r(r(y9i3J)j)j)j)j)^ 

P2, V, b pi, |x, a 

Figure 4.1: Coupling to two gluons. The solid circle either represents or and a and 
b are the colour indices of the gluons. The momenta pi and p2 are incoming. 



In the applications we discuss below the operators O3, O4 and and the corre- 
sponding coefficient functions only appear in bare form. Thus we will not specify their 
renormalization, which can be found in |]7T|, and concentrate in the following on Oi, O2 
and Oq. The relation between the bare and the renormalized operators reads 

On = T.ZnmOl, (4.5) 

m 

where the indices n and m adopt the values 1, 2 and 6. Due to the equality 

E^M = E^^A, (4.6) 

n n 

we obtain the renormalization prescription for the coefficient functions as 

= (4,7) 
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The anomalous dimension matrix pertaining to Znm is defined through 



d/i- 



It is connected to the renormahzation matrix through 



In 



k \ 



and reads in exphcit form |]69|, |70|, ^ 



7 



/ 



V 



a 13 

dots as 



d/i2 



4a. 








■'mk 








kn 



(4.9) 



-470 
47m 



(4.10) 



Thus 7„m can be expressed in terms of the functions I3{as) (see Eq. (|2.35|) ), 7m (os 
(see Eq. ( |2.12| )) and the anomalous dimension of the vacuum energy 70 which is given 
by P, 



37| 
7o 



TT 



A as fas 
1 + + — 



3 TT 



TT 



313 2 5 

~ 3^^ ^ 12^^. 



+ 



14251 77 19 

1296 ~ y ^ "e"^^ 



1975 



54 



4109 _ 35 _ 16 
'1944 ~ 54^^ ~ 



341 2^ 



f \ 1458 9 



(4.11) 



where n/ is the number of active quark flavours. 

In the formalism presented above it is assumed that there are only currents coupling 
to quarks of the same flavour which is sufficient for the purpose discussed below. The 



more general case involving also non-diagonal terms can be found in Ref. [74 



4.1 Hadronic Higgs decay 

The coupling of a scalar CP-even Higgs boson to quarks has a very simple structure as 
no 7 matrix is involved. It essentially consists of a factor containing the mass of the 
quarks and the coupling constant. This fact makes it very simple to construct an effective 
Lagrangian and to derive powerful low-energy theorems as we will see in this Subsection. 

We want to consider QCD corrections to the hadronic decay of a Higgs boson in the 
so-called intermediate mass range, that is we compute the partial decay widths into quarks 
and gluons assuming that the top quark is much heavier than all other scales involved in 
the process. Although we have in mind the top quark as the heavy particle we will in the 
following consider a generic heavy quark with mass irih. 

To lowest order in the inverse heavy quark mass the effective Lagrangian is constructed 
by the operators of dimension four which are discussed above. As we consider light quark 
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mass effects at most in leading order (which corresponds to an overall quadratic factor in 
the case H qq and to quark mass zero in the gluonic case) there is no contribution from 
the operator Oq and thus the anomalous dimension matrix becomes two dimensional. 

We consider a theory where we have in addition to the QCD Lagrangian of Eq. (|3.1|) 
a scalar particle, if, which couples to fermions via 



1 





9 ' 



(4.12) 



where the sum runs over all quark flavours. In the limit oo Eq. ( 4.12 ) can be written 

as a sum over the operators given in Eqs. ( [4.1|) and ( [4.2|) accompanied by coefficient 
functions containing the residual dependence on the top quark: 



y,eff 



rrO 5 



(4.13) 



The relation of the bare coefficient functions and the operators to their renormalized 
counterparts can be extracted form the anomalous dimension matrix given in Eq. ( [4.10|) . 
One obtains for the renormalized operators 



O2 



1 + 2 



a' 



,9 



da' 



InZ; 



(9° 



2 5 



0°-4 



da' 



\nZL 



2 5 



(4.14) 



and accordingly for the coefficient functions 



C2 



1 



l + 2{a'^d/da'^)\nZ, 
A{a',d/da',,)\nZ'^ 



7^1 



l + 2{a'^d/da'^)\nZ, 



7C1 



Consequently, the physical part of CY,eS takes the form 



^phys 
'~'Y,eS 



-^(CiOi + CsOs) 



(4.15) 



(4.16) 



a'g, Z'g and are defined in the effective theory which is indicated by the prime. Ci 
and Oi {i = 1,2) are individually finite, but, with the exception of [Oj]; they are not 
separately renormalizat ion-group (RG) invariant. In Ref. |7^, a RG-improved version 
of Eq. ( 4.16|) has been constructed by exploiting the RG-invariance of the trace of the 
energy-momentum tensor. The ratio H^/v^ receives a finite renormalization factor, which 
is of OiCpMl). Its two- and three-loop QCD corrections have been found in Refs. 
and |]7^, respectively. 

A closer look to the operators in Eqs. ( |4.1| ) and ( |4.2| ) and the effective Lagrangian in 
Eq. ( |4.13|) suggests that there should be a connection between the coefficient functions 
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on one side and the decoupling constants as introduced in Eq. ( p.4| ) on the other side. 
Actually it turns out that Ci and C2 are obtained from derivatives of ln(^g and \n.C,mi 
respectively. This nice feature allows for the computation of Ci and C2 to order a" from 
the knowledge of and Cm to order as we will see below. 



The coefficient functions contain the remnants of the heavy quark. It is thus quite 
plausible that their computation can be reduced to Feynman diagrams where the only 
scale is given by m/^. Actually, the philosophy for the derivation of formulae for C°, . . . , C5 
is very similar to the case of the decoupling relations. One again considers Green functions 
in the full and effective theories and exploits the fact that the coefficient functions do not 
depend on the momentum configuration. We will see that for a certain Green function 
in general several coefficient functions are involved. Thus we will obtain five equations 
which can be solved for C°, . . . , C5 . 

From Eqs. (4J.) and (^4.2|) one learns that the Green function involving one gluon, two 



ghosts and a zero-momentum insertion of the operator Oh = m'lhPh^, which mediates the 
coupling to the Higgs boson, only gets contributions from O^. We define the bare IPI 
Green function in analogy to Eq. ( |3.7| ) 

P'9's{-^r'' [rL-cO.(P,^,0)]+...} 
= i2y'dxd?/e^(^'-^+'=-^)(Tc°''^(x)c°'''(O)G°'"'^(2/)0;,(O))^^\ (4.17) 

where the ellipses again represent other colour structures and p and k are the outgoing 
momenta of the c and G, respectively. The third argument of Tq^^q indicates the zero 
momentum of the operator Oh- 

In a next step we express Oh in terms of the operators given in Eqs. ( [4.1| ) and ( [4.2| ). 
Only has to be taken into account as only this operator contains the coupling of a 
gluon to two ghost fields. In the transitition to the effective theory also Eq. ( ^.41 ) has to 

be considered which leads to a factor Cg/iCs \/C^)- Note that in Eq. ( ^.17| ) we are dealing 
with amputated Green functions. Thus from the external propagators a factor (C3)^C3 
arise. Finally we arrive at 

'^"''i-ir^'yc^ci + ... , (4.18) 



p'^g. 



s " 



where the ellipses represent other colour structures and terms suppressed by the inverse 



heavy quark mass. On the right-hand side of Eq. (|4.18|) we also avoided to write down 
explicitly the contributions beyond tree-level as they vanish for p, A; ^ within dimen- 
sional regularization. In this limit only the diagrams involving the heavy quark survive 
on the left-hand side and we are left with the formula 

Clct = r^^,^^ (0,0,0). (4.19) 

In order to reduce the number of contributing diagrams and also to simplify their 
complexity we exploit that the coupling of the Higgs boson is proportional to nih. Thus 
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it can simply be generated by differentiation with respect to rrih. With the definition 

~ d{mlf' ^ ^ 

we finally obtainQ 

ClCl = 29°r^^,(0, 0) . (4.21) 

As another example, let us consider the derivation of a formula involving and 
C4. The starting point is the IPI Green function of two gluons which contains a zero- 
momentum insertion of the composite operator Oh- In momentum space, it reads in bare 
form 



IPI 



5 



ab 



g^^ p^V^QQQ ijp') + terms proportional to p'^p" , (4.22 



where p denotes the four- momentum flowing along the gluon line. In the limit 00, 
Oh may be written as a linear combination of the effective operators given in Eqs. ( |4.1| ) 
and so that 



= ^ y dxdy -(^-2^) (tG°''^''^(x)G°''^'^(i/) + C°0°) )^^^ + . . . 



IPI 

+ 



= -g^^P Qi-^^C^ + 2C4') (1 + higher orders) + . . . , (4.23) 

where the ellipses indicate terms of 0{l/mh) and terms proportional to p'^p". The factor 
1/8 results for the summation over the colour indices. In the second step, we have used 
Eq. (|3.4|) together with the fact that Tqq'q^ {p, —p) represents an amputated Green func- 
tion. If we consider the coefficients of the transversal part in the limit p — > 0, we observe 
that the contributions due to the higher-order QCD corrections on the right-hand side 
of Eq. ( [4. 231 ) vanish, as massless tadpoles are set to zero in dimensional regularization. 
In principle also the other operators contribute via loop diagrams. However, also these 
contributions lead to massless tadpoles and are thus zero. The relative weight between 
and C4 and the prefactor in the last line of Eq. ( [4.23| ) follow immediately from the Feyn- 



man rules given in Eq. (4^). On the left-hand side, only those diagrams survive which 
contain at least one heavy-quark line. Consequently, the hard part of the amputated 
Green function is given by 

r^U(O'O) = C3° (-< + 2Q°) . (4.24) 



^^We want to mention that in Ref. p9| there are misprints in the corresponding formulae: erroneously 
they contain a factor "1/2" instead of "2" . However, the initial equations and the final results are correct 

in Ggli. 
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Using Eq. ( [4.20| ) we finally arrive at 



Cl{-ACl + 2Cl) = -2dlu'^{Q). (4.25) 

In a similar way, we obtain three more relationships, namely 

C°C3° = -2a°S°/(0), 
ClCliC'.-Cl) = 1 - (0) - 29°S°/(0) , 

C.l{Cl + Cl) = 29M'(0). (4.26) 

The Eqs. (|4.21|) , ( [4.25| ) and ( [4. 261) may now be solved for the coefficient functions C° 
(z = 1, . . . , 5). They are expressed in terms of vacuum integrals with only one mass scale, 
namely the heavy quark mass. This is also the case for the decoupling constants occuring 
in Eqs. (^^ ) and ( ^:26D as discussed in Section |. 



As a simple example let us consider the computation of Ci at lowest order. Here only 
one diagram contributes to ng'^(O), namely the one where the gluon splits into two virtual 
heavy quarks. It can be evaluated with the help of Eq. ( [A.6|) . Expanded up to finite order 
in e it reads 

nS» - ^^^(i.lln^). (4.27) 

Differentiating with respect to m° leads to 

a T 

C, = (4.28) 

TT 6 

Note, that at this order no renormalization is needed and thus the bare and renormalized 
quantities coincide. Furthermore, there is no contribution from in Eq. ( [4.25|) and 
contributes for the first time at order as can be seen from Eqs. (^4.261) . 

In the next three Subsections we will describe different methods to compute the co- 



efficient functions. Finally in Section [4.1.4| we will review the state-of-the-art corrections 



to the hadronic decay width of the Higgs boson. 



4.1.1 Direct calculation of the coefficient functions 

There is the possibility to avoid the occurrence of the non-physical operators and their 
coefficient functions in the computation of Cj* and For demonstration let us consider 
the case of Cf. 

From the definition of the operators in Eqs. ( [4.1|) and ( ^4.2|) and the Feynman rules 
given in ( [4. 4] ) one can read off that the IPI Green function Tqq'q iPi,P2), where pi and 
P2 denote the incoming momenta of the gluons, only receives contributions from C^Oj* if 
Pi 7^ —p2- In analogy to Eq.( [4.22|) we can write 

^GGoSPuP2) = {9'"' Pi ■ P2 - P\P2) rGGO, (Pl.P2) 

= {^Pi ■ P2 - Pm a (-4C°) (1 + ■ ■ ■) , (4.29) 
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where the elhpses denote terms of order l/ruh, possible other Lorentz structures and 
higher order loop corrections. The latter vanish in the limit Pi,P2 and we are left 
with 




Pl,uP2,f, - Pl,fiP2,v 



iD-2){p,-p2) 



GGOh 



iPl,P2] 



(4.30) 



Pl=P2=0 



where the on-shell conditions pi = P2 = ^ have been used. This equation relates the 
coefficient function Ci directly to a physical amplitude and no Green functions involving 
ghost fields have to be considered. The price one has to pay is that the momenta pi 
and P2 can only be set to zero after the projection in Eq. ( |4.3CI| ) has been applied. This 
complicates the expressions of the individual diagrams. Furthermore it is not possible to 
use derivatives with respect to m° in order to generate the coupling to the Higgs boson 
as initially the momentum of the Higgs boson is not zero. Thus altogether 657 vertex 



diagrams like the first one in Fig. [4.2| have to be considered. 





I 



,6 





Figure 4.2: Typical Feynman diagrams contributing to the coefficients C° and Cg. 
Looped, bold-faced, and dashed lines represent gluons, heavy quarks, and Higgs bosons, 
respectively. 



In analogy to the gluon-gluon-Higgs boson three-point function one could also choose 
the part of the operator Oi involving three or even four gluons. This would lead to 
7362 four- and 95004 five-point functions at three-loop level, where one and no gluon 
momentum, respectively, has to be kept different from zero until the projection is finished. 
Sample diagrams are pictured in Fig. |]2[ In Ref. ||78|, where the decay of a pseudo-scalar 



Higgs boson has been considered, the 7362 three-loop four-point functions have been 
evaluated in order to check the calculation of the three-point function. For completeness 



we want to mention that the last diagram in Fig. [4.2| is the lowest-order graph contributing 
to C2. 

4.1.2 Low-energy theorem 

It is tempting to re-express the Green functions on the right-hand side of Eqs. ( |4.21| ), ( |4.25| ) 
and ([4.26|) in terms of the decoupling constants which were considered in Section 0. This 
leads to a straight relation between the coefficient functions and the decoupling constants. 
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For the equations involving and (7° one successively obtains 



2^ 

\ a c? ) ' 
1 did ^ d^ci die. 



2 a 

-9° In ■ 



S3 V S3 





9 ' 



(4.31) 



Next, we express ^ through renormalized quantities. Using = dh, we find 



2C° = d,ln{C 



0\2 



a 



0/ 



a), 



dhHZ'Y + dhlna', 



1 ^ "^^1 r^''^ 



9/,, In al . 



(4.32) 



Identifying the renormalization factor of Eq. ([4.15|) , we obtain the amazingly simple rela- 
tion 



-2Ci 



dh In a'^ 



(4.33) 



This relation opens the possibility to compute Ci through C(af), since one only needs to 
know the logarithmic contributions of (g in this order, which may be reconstructed from 
its lower-order terms in combination with the four-loop (3 [|1^] and 7^ [0, functions. 

It is furthermore possible to directly relate Ci to the (3 and 7^ functions of the full 
and effective theories. Exploiting the relation 



da. ^ dnih 



(4.34) 



where a'^ = a'^ifi, ag, rrih), we find 
Cl - 



TT 



2a', [1 - 27™(«,)] 



(3'{a',)-(3{a, 



I doi's 

das 



(4.35) 
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In the case of C2, we may proceed along the same hnes to obtain 

C2 = 1 
= 1 



29/jlnCm 
2 



1 - 27^ (a. 



1 Qth! 



m' da. 

q t, 



(4.36) 



where m'^ = m^(/i, a^, m^j). It should be stressed that Eqs. (|4.35| ) and ( [4.36D are valid to 
all orders in a,. 



Fully exploiting present knowledge of the (3 and 7^ fT^, |T6[ functions, we may 
evaluate Ci and C2 through O^a^) via Eqs. ( |4.35| ) and ( [4.361 ). In the pure MS scheme, 



we find p 
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where, for simplicity, we have chosen ^ = fi^ = mh{^h) in the approximate expressions. 
The corresponding results expressed in terms of the pole mass read 
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where we have put ^ = in the numerical evaluations. In the leading Yukawa 
corrections of 0{a'^Gpm'^) {n = 0,1,2) to the coefficient functions have been evaluated 
in the SM. The analytical results are listed in Appendix ^ 

With the knowledge of Ci and C2 the construction of the effective Lagrangian is 



completed. In Section [4.1.4| it will be used in order to compute the hadronic decay rate 
of the Higgs boson. 

Recently the effective Lagrangian has been used in order to consider the Higgs boson 

m 



a ffist step to the next-to- 



production process via gluon fusion. In Refs. 
next-to-leading order QCD corrections has been done. At this accuracy it is necessary to 
compute two-loop virtual corrections to the process gg ^ H using the effective ggH 
vertex of Eq. ( |4.16| ). Thus the coefficient function Ci enters as a multiplicative constant. 

For completeness we want to mention that the low-energy theorem derived in this 
section has also been specified to the 'j'jH coupling. The analytical expressions can be 



found in Ref. 49 . 



4.1.3 H ^ gg in the background field method 

An interesting alternative to the considerations of the Sections [4.1.1| and [4.1.2| is based 



on the Background Field Method (BFM) [^3|, ^ |8^, In this framework the gluon 
field is decomposed into a quantum and a background part where the former only appears 
as a virtual particle inside the loops. On the contrary, the background field serves as an 
external field in the Green functions. In the BFM the gauge invariance is maintained 
while quantizing the theory. This was one of the main motivations for its development. 

A comprehensive discussion for the case of QCD and the computation of the two-loop 
f3 function as a practical application of the BFM can be found in p3. In particular, it is 
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shown, that the relation between the charge renormahzation constant, Zg, (cf. Eq. ( |3.3| )) 
and the wave function renormahzation constant of the background gluon, , reads 

= ^- (4.41) 

Thus, only the background gluon self energy has to be computed in order to obtain the (3 
function whereas in the conventional approach also vertex functions have to be considered. 
E.g., next to the gluon propagator also the ghost two-point function and the gluon-ghost 
vertex have to be computed. The price one has to pay for these kind of simplifications are 
more complicated Feynman rules for those vertices involving the background gluon . 



In particular, there are three- and four-gluon vertices which contain additional terms 
proportional to 1/(1 — where ^ is the gauge parameter defined in Eq. p.2| ). As a 
consequence, it is not possible to choose Landau gauge, which corresponds to ^ = 1 and 
which would avoid the renormahzation of ^, from the very beginning of the calculation. 
Rather one has to perform the calculation using a general gauge parameter and either 
adopt Landau gauge at the end of the calculation or renormalize the gauge parameter. 

Concerning the decoupling properties of the background field we can essentially take 
over the discussion of Section ||. In analogy to Eq. (|3.10| ) one obtains for the decoupling 
constant of the background gluon field 

'BO _ 1 I -nO,h 



C^' = l + n^^(0), (4.42) 



where 11^^ is the hard contribution of the transversal part of the bare background gluon 
polarization function. coincides with at one-loop order as only the diagram with 
a heavy quark loop contributes and the coupling of fermions to background gluons is 
identical to the one of quantum gluons. Starting from two loops, however, virtual quantum 
gluons appear inside the Feynman diagrams and the analytical expressions are different. 
In Fig. ^]3| some typical one-, two- and three-loop diagrams contributing to n^s(O) are 



shown. One arrives at a finite expression for the decoupling constant with the help of 



Cs" = (4.43) 

It is not surprising that one has the following connection between and (g 

a = (4.44) 

which is in analogy to Eq. (|4.41|) . We explicitly checked this relation and computed 
11^^(0) at three-loop order using a general gauge parameter. After renormahzation one 
obtains 



a 




24 24 36 
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Figure 4.3: Typical Feynman diagrams contributing to the two-point function with exter- 
nal background gluons. Looped, solid, and dashed lines represent gluons, heavy quarks, 
and ghosts, respectively. 
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1728 mJ^ 



(4.45) 



which is in agreement with Eq. ( p.30| ). The 219 diagrams contributing to 11^^(0) have 
to be compared with the 189 -1-25-1- 228 = 442 diagrams which are necessary in order to 
obtain the result of Eq. ( |3.30| ). 

Eq. ( [4.44|) also significantly simplifies the computation of the coefficient function Ci. 
Due to Eq. ( [4. 331 ) it is simply obtained via 



(4.46) 



This coincides with the naive expectation that the effective coupling of the Higgs bo- 
son is generated by taking derivatives of the gluon polarization function with respect to 
the heavy quark mass. Furthermore Eq. ( [4.46|) agrees with the low-energy theorem for 
the photon-Higgs interaction as derived in In this sense we could claim that the 

background gluon is "more physical" than the gluon in the conventional approach. 

In this context we want to refer to |8^ where low-energy theorems in Higgs physics have 
been considered at one- and two-loop order. In particular it was realized that the decay 
rate of a scalar Higgs boson into two photons can be obtained by naive differentiation of 
the photon self energy if the latter is computed in the framework of the pinch technique. 

As a further check on the result for Ci we also perform the direct calculation by 
considering the quantity V^bqbq^{pit'P2) which is defined in analogy to Eq. ( [4.22| ): 



GBGBOh 
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The only difference is the presence of external background (instead of quantum) fields. 
We evaluated the 732 vertex diagrams in the spirit of Eq. (|4.30| ) and could confirm the 
order terms in Eqs. ( |4.37| ) and (|4.39| ). 

Due to the fact that there are no diagrams involving external ghost fields contributing 
to Eq. ( |4.46| ) there is no admixture from 0° while matching T'^b^bq^{pi,P2) with the 
effective theory even for the choice pi = —p2- Thus we can immediately set pi = —p2 = p 
in Eq. (|4.471 ) and consider the limit p 0. No complicated projector like in Eq. (|4.30|) 
is necessary which simplifies the evaluation of the diagrams with external background 
gluons. 



4.1.4 Hadronic decay rate of the SM Higgs boson 



In the following we present the current status of the total decay rate of the SM Higgs bo- 
son]^ in the intermediate mass range. This is done on the basis of the effective Lagrangian 
Eq. ( [4. 161 ), where the top quark takes over the role of the heavy quark h. 
The total decay rate into hadrons can be cast in the form|3 



r (i^ ^ hadrons) = (1 + | ^ A,,- [(1 + AI2) (^2^)' + A?2 W 



2q 



Ar + A, 



weak 



xt=0 



(4.48) 



with Aqq = 3GFMHml/{47cV2) and Agg = AG pMjj/ {71^/2). The terms in the first hue 



and the first term in the second line proportional to Agg have their origin in Eq. ( [4.16|) . 
In particular, the universal corrections 5u arise from the renormalization of the factor 
/v^ and are known to order a^GpfTif [f7]. The factors A^- contain the QED and QCD 
corrections from the light degrees of freedom only, while the terms A'^'^° summarize the 
corrections coming from higher dimensional operators. They are at least suppressed by 
a factor a'^M'^ / . In Eq. ([4.48|) we separately display the weak contribution where the 
leading term of order Gprn^ is stripped off. It is denoted by|^ ^weak 

Typical diagrams contributing to the QCD corrections of A22 are pictured in Fig. |474 . 
At one-loop order they have been evaluated in 



1[ keeping the full dependence of the 
quark mass. The dominant corrections at order 0{a1), i.e. those obtained keeping only 

The calculation 



9|. 



the factor from the Yukawa coupling, have been evaluated in 
has later on been improved and the correction terms proportional to m^/M^ became 
available []93| , A naive expansion in rriq is sufficient for their computation. Beyond 
the quadratic term, however, one either has to adopt the method discussed in Section ^]2] 
or apply the large-momentum expansion. The latter has been performed in |p5[ and 



iSpQj. j-gyjg^ articles we refer to |8^ . 

^*The additional index "q" for the coefficient function C2 indicates that there might be an explicit 
dependence on the flavour through elektroweak corrections. 



^^For the definition of xt see Eq. (B.7) 
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Figure 4.4: Typical Feynman diagrams contributing to A22. The solid circles represent 
the operator 



mass correction terms up to order (m^/M^)* have been evaluated. For a Higgs boson 
in the intermediate-mass range considered in this Subsection the corrections beyond the 
quadratic term are quite small and can safely be neglected. In the higher order terms 
have been considered in the context of a heavy Higgs boson which can also decay into top 
quarks. In this case it turned out that even the quartic terms are important and only an 
expansion up to order {mf/Mfj)^ gives satisfactory results. In Section ^]2| we review the 



calculation and numerical results are presented in Tab. [5.1|. For completeness we want 



to mention that the imaginary part of the correlator {O2O2) has even been evaluated at 
four-loop order using the technique described in Section p.2| . This leads to corrections of 
order a? 



|0 



Next to the pure QCD corrections there are also the contributions of order a and the 
mixed QED/QCD terms of order aa^ p6| , p7|] which can easily be extracted form the 



QCD results. In summary, the numerical result for the discussed terms read 
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(4.49) 



where rii is the number of light quarks. 

The imaginary part of the mixed correlator {O1O2) is denoted by Af2- Next to contri- 
butions to the partial width into quarks it also involves purely gluonic final states as can 
be seen from the diagrams shown in Fig. |4.5| . The contributions at two- and three-loop 

respectively. In numerical form AI2 is 



order have been computed in 
given by 



B1I981, 1751 and M 
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(4.50) 



In case one is only interested in final states involving quarks the purely gluonic cuts have 
to be subtracted. Currently they are only known at order a^. To this order the subtracted 
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Figure 4.5: Two- and some of the three-loop diagrams contributing to A^2- The sohd 
circles represent the operators 0\ and C2, respectively. 



result reads 75 



TT 



76 .4, 



Ml, 



+ 



(4.51) 



where the ellipses indicate terms of order (x\. Note the logarithmic singularity in the 
light-quark mass which arises from the fact that only parts of the final state are contained 
in Eq. (|43l| ). 




Figure 4.6: Typical Feynman diagrams contributing to the correlator {(D\0\). Looped, 
solid, and dashed hues represent gluons, light quarks, and Higgs bosons, respectively. 
Solid circles represent insertions of 0\. 

The correlator formed by the operator 0\ mainly contains cuts arising from gluons. 
However, starting from two loops there are also contributions from light quarks and at 
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order there are even cuts involving no gluons at all (cf. second diagram in the lower row 
of Fig. |4.6| ). In Fig. ^]6|some typical diagrams are pictured. In particular, the combination 



(Ci)^Aii contains the contribution from all the diagrams pictured in Figure O. The two- 



loop contribution has been evaluated in ||100| , |101|| and the order al terms can be found 
if ||79[. If we set /^^ = and evaluate the correlator for = Mfj we obtain 
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af\MH) 



vr 



;i8.250 - 1.167ni) 



"^ll^^MllD^ (242.973 - 39.374 rn + 0.902 n'f) . (4.52) 




Figure 4.7: Typical diagrams generating 0{a'^) corrections to r{H gg). Bold-faced 
(dashed) lines represent the top quark (Higgs boson). 



The contributions to A^*^° are not covered by Eq. ( ^.16| ). In the language of the effective 
Lagrangian it would require to deal with operators of dimension six and higher. However, 
up to now they have not been studied in detail. The approach adopted in Refs. [Bl, PB| is 



based on asymptotic expansion which is applied to the propagator-type diagrams involving 
a top-quark loop. At order there are two classes of such diagrams contributing to A^*^", 
namely the double-triangle (or singlet) and the double bubble diagrams. The exact result 
for the imaginary part of the latter with massless external quark lines and heavy virtual 
top quark can be found in 



102| 



The leading term contribution to A^*^" reads 



A 



hdo 



a 



(5)- 



TT 



a 



(5)- 



Ml 



Mi 



5863 
24300 



113 , Mjr' 

In- ^ 



1620 M^ 
^ ' 0.241 - 0.070 In " 



Mi 



(4.53) 



53 



The higher order contributions in Mfj/M^ are very small and can safely be neglected 
even for = Mf. For Ajj*^" one obtains 
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where again the higher order terms in Mfj/M^ are much smaller. The contribution to 
A^'^" is obtained from the application of the hard-mass procedure to the propagator- 
type diagrams where the external Higgs bosons are attached to top quarks and they 
subsequently couple to gluons. To leading order in the expansion for large rrit one obtains 
the one- and two-loop diagrams of Fig. 4^. The comparison of Eq. ( [4.54 ) with Eq. (|4.52|) 
shows that also in this case the contributions from the higher dimensional operators are 
small. We want to mention that for H ^ gg a numerical calculation of order exists 
which takes the complete mass dependence into account [|101[ . 



At one-loop the complete weak corrections have been computed in analytical from ||103 
T04| . If we put m„ = and consider the limit Mh < 2Mw, A"*'''''|^,=o takes the form [|104 
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(4.55) 



The leading m\ term is stripped off as it is already contained in the universal factor 



Expressed in terms of the MS top quark mass the latter reads in numerical form [|105| , [77 
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We want to mention that also non-universal radiative corrections to Ci and C25, which 
are enhanced by a factor GpiTif, are available up to the three-loop order [77, 75, 59 1. They 
will be listed below in comparison with the terms of 0{al). 

At this point we want to compare the relative size of the individual terms. In particular 
we have in mind terms of order 0{al), a^ml/Mfj, aa^, a'^Xf and a^Mjj/Mi. For this 
purpose we will consider all quarks with mass lighter than mi, as massless. This means 
that the sum in Eq. ( 4.48 ) reduces to q = b which is conveniently written in the form 



T{H hadrons) = Af,i (l + A^^ + A^) + 
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where contains only corrections from light degrees of freedom. All top-induced terms 
proportional to ^4^5 from Eq. ( 4.48| ) are contained in At, which we express in terms of 
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af'\n). Ag contains the corrections from the gluonic final state. Choosing jj^ = Mfj and 
n/ = 5 we find 
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(4.59) 
(4.60) 



with a^lp 



a(5)(MH)/7r, Lt = InM^/M^, and Xt = GpM^ /{Stt^V^). In Eqs. (|438|) 



and ( [4. 591 ) also the quadratic mass correction terms are listed. In A; they are obtained 
from the naive expansion of the diagrams. The (m[^^)^/M^ corrections in At arise from 
the singlet diagram with one top and one bottom quark triangle. In this case a naive 
expansion fails as can be seen by the logarithmic term in Eq. ( |4.59| ). Instead the asymp- 
totic expansion has to be applied [Q. Both for A/ and Aj one observes that the 0{a^) 
term proves to be numerically more important than the power suppressed contribution of 
0{a1ml/Ml). Note, that Eq. ( [4.59|) contains contributions with pure gluonic final states 



which is due to diagrams of the type in Fig. |4~5 . 

In the approximation considered here we have —2 ^ Lf < 0. This means that the 
logarithm needs not necessarily to be re-summed as in addition the coefficients in front 
of Lf are much smaller than the constant term. 

A comparison of Eqs. (|4.58| ) and ([4.59|) shows that the top-induced corrections in 
At of 0{al) are numerically of the same size as the ones arising from "pure" QCD. 
Furthermore one should mention that the coefficient of the Mt-suppressed terms are tiny 
and, as OLsjXt ~ 30, also the oP^Xt enhanced terms are less important than the cubic 
QCD corrections. This is also the case for Eq. ( |4.60| ). For comparison in Eq. ( |4.58|) 
also the two-loop corrections of order aag are listed. In principle also higher order mass 
corrections are available [Q. However, in the case of bottom quarks it turns out that 
they are tiny. 

In summary, we have shown that for an intermediate-mass Higgs boson the application 
of the effective Lagrangian (cf. Eq. ([4.13|) ) is quite successful and enables the computation 
of the hadronic Higgs decay up to orders and a\ for the quark and gluon final states, 
respectively. The smallness of the higher dimensional operators (cf. Eq. ( [4.53|) and ( [4.54| )) 
justifies this approach. In conclusion we can state that the perturbative expansion of the 
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hadronic width of the Higgs boson is well under control. 



4.2 Quartic mass corrections to cr(e^e hadrons) 

The total cross section for hadron production in electron-positron annihilation, o"(e'^e^ 
hadrons), is one of the most fundamental observables in particle physics (for a review 
see ||106| |). For energies sufficiently far above threshold it can be predicted by perturbative 
QCD, and it is well accessible experimentally from threshold up to the high energies of 
LEP and a future linear collider. It allows for a precise determination of the strong 
coupling as and, once precision measurements at different energies are available, for a 
test of its evolution dictated by the renormalization group equation. 

Often the center-of-mass energy is much larger than the quark masses which then can 
safely be neglected. However, there are also many situations where it is important to 
take into account the effect of finite quark masses |p,07|| . E.g., one can think of charm or 



bottom quark production not far above their production thresholds ||108| , |109|| , or of top 
quark production at a future linear collider | |110| |. 



The complete mass-dependence at order to R{s) = a{e~^e~ hadrons)/cr(e"'"e~ 
has been evaluated quite some time ago in analytical form [|1 1 1|| in the context 
of QED. At order a1 this task already becomes much less trivial. The massless approx- 
imation became available quite some time ago [^]. However, only for a certain class 



of diagrams — the ones containing a second massless quark pair — the full quark mass 
dependence could be obtained in analytical form using conventional methods ||112|| . For 
all other contributions different methods have to be applied which are discussed in detail 
in Section |^. Here, we only want to mention that a crucial ingredient is the application 
of the large-momentum procedure which provides an expansion in m?/s. At order 
also this method fails as it would be necessary to evaluate massless four-loop propagator- 
type diagrams. At the moment this is not possible. Thus a different strategy has to be 
employed which we will describe below. 

In addition to the massless result, which has been obtained in [^, the m? / s 
terms of O(a^) have been calculated in ||113|] . They were obtained by reconstructing 



the logarithmic a\w? j s terms for the polarization function H((3'^) from the knowledge of 
the full three loop 0[oPi.w? j s) result of ||1 14|| with the help of the renormalization group 
equations. These are sufficient to calculate the j s terms of the imaginary part in the 
time-like region. A generalization of this approach has been formulated for the quartic 
mass terms in [^2], ^ and was originally adopted for the calculation of a^m^/ terms ||7^ . 



The basic ingredients are the OPE and the renormalization group equations 



(RGE). The idea is to apply the OPE to the correlator of two currents and compute 
its imaginary part which immediately leads to corrections for The current corre- 

lator is expressed as a sum over local operators multiplied by coefficient functions which 
represent the short distance part of the process. Afterwards one exploits the RGEs in 
order to relate different pieces in the sum and to construct the logarithmic terms of the 
polarization function. In addition to the anomalous mass dimension and the /3-function, 
the anomalous dimensions of the operators of dimension four are required in appropriate 
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Figure 4.8: Large-momentum procedure for the one-loop photon polarization function. 
The quark lines carry the mass m. 



order. 

Let us be more specific and consider the time-ordered product of two currents. The 
application of the OPE leads to 

T^{q) = t fdVTj{x)j{0) Y.Cniq)On, (4.61) 

where the dependence of the coefficient functions Cn{q) on the large scale q is made 
explicit. In Eq. (|4.61|) we have to consider all operators On of dimension four as exactly 
those contribute to the quartic mass corrections. They can be found in Eqs. ( [4.1|) and ( [4.2|) . 

Actually the contributing operators essentially coincide with the ones of Section |4.1| . 
There, however, the top quark mass took over the role of the large scale and the external 
momentum was supposed to be much smaller. As a consequence the coefficient functions 
depend on Mt and are computed with the help of vacuum integrals. On the contrary, we 
will see below that the coefficient functions of Eq. (|4.61|) are expressed through massless 
propagator-type integrals. 



The virtue of Eq. ( [4.61| ) becomes obvious if one takes the vacuum expectation value. 



Then the left-hand side turns into the polarization function which in the following we 
generically call n(g'^). On the right-hand side we obtain a sum of vacuum expectation 
values of the local operators multiplied by the coefficient functions. As already mentioned, 
the latter only depend on q whereas the scale of the vacuum expectation values is given 
by the quark mass. This strongly resembles the large-momentum procedure. Also there 
a factorization of the scales is achieved. In fact, it can be shown that for a certain choice 
of the operator basis an identification of the Cn{q) with the hard subgraphs and of the 
vacuum expectation values with the co-subgraphs is possible. 

For illustration we consider the correlator of two vector currents at one-loop order. The 



application of the large-momentum procedure is visualized in Fig. |4.8| . The prescription 
tells us that the first term on the right-hand side of the equation has to be expanded in 
the quark mass, m, leading to massless integrals. Also the subdiagram on the very right 
in the second term has to be expanded in m and then has to be inserted in the blob of 
the loop-diagram. Thus one ends up with a vacuum integral. Taking into account terms 
up to quartic order one gets^^ (see, e.g., Ref. [ |115| , |116|| ) 




^^For a precise definition of the polarization function see Eqs. ( 5.15 ) and ( 5.20 ) in Section || 
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(4.62) 



with /g^ = ln(— g^//i^) and Igm = ln(— g^/m^). The elhpses in Eq. ([4.62| ) denote higher 



order mass correction terms. On dimensional reasons it is clear that the m° and 



term only comes from the first diagram on the right-hand side of Fig. |4.8| . Also the 
m^ln(— g^) term can only arise from this diagram. On the other hand, the massive 
logarithm, m^lnm^, originates from the vacuum integral of Fig. ^]8. 

Let us now look at the result provided by the OPE. At lowest order in only the two 
operators O2 and Oq of Eq. (|4.1|) contribute. The product of the corresponding vacuum 



expectation values and the coefficient functions reads (see, e.g., Ref. ||115| , Q) 




Thus, Cg(Og) exactly repoduces the m^ln(— g^) of the massless diagram of Fig. |47^ , 
whereas {O2) provides the massive logarithm. The sum of the two contributions in 
Eq.( [4.63D reproduces the m*^ terms of Eq. ( [4.62|) . 



Note that in the above consideration no normal-ordering prescription has been used. 
Otherwise the vacuum expectation value of the operator O2 would be zero. Furthermore, 
the coefficient function Cq would necessarily contain Inm^ terms in order to reproduce 
the quartic terms of Eq. (|4.62| ). Thus, in case the normal-ordering prescription is applied, 
there is no separation of the two scales q and m. 

This example, in particular Eq. ( [4.63| ) and Fig. shows that in principle one 
could still use the large-momentum procedure for the practical computation of Cn{q) 
and {On{0)). However, in practice it turns out that this is quite tedious. 

The aim of the calculation is to obtain R{s) up to order m^a^, which means that due 
to the equation 

R{s)\ ^ = 127rlmn(g2 = s + ze) 

" ~ Im [C^{0^) + ^2(02) + ^^6(^6)] , (4.64) 

one has to evaluate the coefficient functions and vacuum expectation values up to suffi- 
ciently high order. The imaginary parts can only arise from the coefficient functions as by 
construction only they can develop logarithms of the form ln(— s — ie). Note further that 
the information about the considered current only enters into the coefficient functions; the 
matrix elements of the operators are universal. For the vector current correlator it turns 
out that Ci develops an imaginary part starting from 0{a^). As (Oi) is proportional to 
as there is no contribution to order from this term. The lowest order of {O2) is and 
the imaginary part for C2 starts at order a^. This implies that C2 is needed up to order 
(three loops) and {O2) up to order as (two loops). So far the occuring integrals are all 
available in the literature. However, in the case of Cg the logarithmic terms up to 0{as) 
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are needed, which would require a four-loop calculation. This can be avoided as we will 
see in the following . 

We consider the renormalization group properties of the polarization function. In 
general n(g^) is not renormalization group invariant. Considering, however, only the 
quartic mass terms we have 







a^terms 



^i''jl^{C,{0,) + C,{02) + C,{0,)) 



(4.65) 



as non-zero contributions may at most have mass dimension two. Using Eq. ([4.10| ) and 



fx 



_d_ 
dj7^ 



+ f3TT 



_d_ 

da. 



(4.66) 



one obtains 
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(4.67) 



with L = ln(/i^/(— g^)). With the help of this equation the logarithmic terms of Cq at 
order a^ can be obtained through two- and three-loop calculations. In particular Cq itself 
appears on the right-hand side of Eq. ( [4.67| ), however, only at order a^ which corresponds 
to massless three-loop integrals. 

As a simple example let us evaluate the order a° term of Eq. ( |4.67| ). In this limit there 
is only a contribution from the last term. With the help of Eq. ( [4.11| ) and C2 = we 
obtain 



167r2 




(4.68) 



which after integration reproduces the logarithmic terms of the renormalized version of 
Eq. ( ra) : 




-r (-4 + 8/, 



(4.69) 



At this point it is instructive to make again a comparison with the large-momentum 
procedure. Applied to the polarization function there is always one term where the 
hard subgraph constitutes the complete diagram Taylor expanded in the masses (cf. Ap- 
pendix A.l). At ra-loop order, i.e. considering QCD corrections to order a""^, this means 



that n-loop massless propagrator-type integrals have to be solved. Actually, as we are only 
interested in the imaginary part of the polarization function only the logarithmic parts 
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of the integrals is needed. This part exactly constitutes the left-hand side of Eq. ( [4.67| ). 



The right-hand side of Eq. (|4.67| ) contains lower-order terms of C^fJ] and contributions of 



(n — l)-loop diagrams. Thus the price one has to pay in order to avoid the computation of 
the (imaginary part) of the n-loop diagrams is the construction of appropriate operators, 
the computation of their anomalous dimension and their coefficient functions. 

For the practical computation of the coefficient functions the so-called "method of 
projectors" [|1171 , ^ is used. For the projectors, 7r„, an appropriate combination of initial 
and final states, \i) and |/), and derivatives with respect to masses and momenta is chosen 
in such a way that one has 



m 



5nm. (4.70) 



Here (9° is one of the operators defined in Eqs. ( [4.1|) and (|4.2| ) and vr^ has the form 




(4.71) 

p=m=0 



It is understood that the nullification of p and m happens before the loop integrals are 
performed. Thus, in Eq. ( [4.70| ) there is only the tree-level contribution; all loop correc- 
tions become massless tadpoles which are set to zero in dimensional regularization. The 
application of the projectors to Eq. ( ^4.61| ) immediately leads to 

C'M = [tKq)] , (4.72) 

which relates the bare coefficient function C° to massless two-point functions. For Cg, e.g., 
the projector is quite simple. It essentially consists of four derivatives of the polarization 
function with respect to m. From the structure of the operators it is clear that the 
projectors for and are more complicated, as the corresponding diagrams also involve 
external gluons and quarks, respectively [[118| , . 

There is quite some similarity between the "method of projectors" and the procedure 
we have used for the computation of the decoupling constants in Section ^ In fact, if 
one tries to construct a projector for one arrives at a similar system of equations as 
in Eq. ^ 0]. 



As already mentioned, the computation of the vacuum expectation values reduces to 



the evaluation of vacuum integrals which have been calculated up to three-loop order [|11£ 
72|, g 0,|7|]. 



We refrain from listing the individual results for the coefficient functions and the 
vacuum expectation values of the operators but provide directly the results for R{s). 
Thereby we want to list the results of those terms which contribute to the production of 
the heavy quark pair QQ via the exchange of a photon which will be denoted by Rq{s). In 
this Subsection we want to list Rq{s) up to order a^m'^/s'^ and postpone the discussion of 



the remaining terms (in particular the full mass dependence at order a^) to Section |5.2.1 



''Note that 7^ and /3 start at 0{as) and 0{al), respectively. 
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It is convenient to decompose the contributions to Rq{s) into three parts 



\ 9 



smg 



(4.73) 



where the sum runs over all massless quark flavours q and Qq denotes the charge of quark q. 
In Eq. (|4.73| ) we distinguish the contributions where the massive quark Q directly couples 
to photon (rg) from the ones where in a first step a massless quark is produced which 
subsequently splits into the massive quark Q. Furthermore, the singlet contributions are 
displayed separately. They arise from diagrams where the external current couples to a 
closed quark line which is different from the one involving the final-state quarks. In the 
following rg^sing will not be considered. Its contribution is numerically small and can, e.g.. 



be found in [121 



Both rq and r^g are expanded in m^/s and can be written as 

tq = ro + rg,2 + rQ,4 + . . . , 
rqQ = ro + rgg,2 + r^Q^^ + • • • , 



(4.74) 



where tq belongs to the massless approximation, while rg^„ and r^g^^ contain the mass 
terms of order mg. A look to the contributing diagrams shows that the contributions to 
r^g 2 and r^g 4 arise for the first time at order a^. However, it can be inferred from general 
renormalization group considerations that the corresponding coefficient of r^g 2 has to be 
zero 



1131 , |122|| , which means that it starts out only at order a^. 



The numerical result for the massless approximation reads |]60|, ^ 

2 



Cts 

1 + — + 

vr 



(1.98571 



0.115295n 



7r 



+ — - 6.63694 - 1.20013 Uf - 0.00517836 nj) + 



(4.75) 



where the ellipses indicate higher orders in as- Uf is the number of active quark fiavours. 
The quadratic mass corrections are given by ||114| , |113| 



ml as 



rqQ,2 



s 
+ 

s 



12 + — f 126.5 -4.33333 rif 



1032.14 - 104.167nj + 1.21819nJ 



TT 



- 7.87659 + 0.35007 



+ ... , 



(4.76) 



and, finally, for the quartic terms we have [[73| , [7^ , |121|| 

/ .2 \ 2 



6 - 22 



TT 
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+ 



+ — (148.218 - 6.5 Ims + (-1.84078 + 0.333333 Ims) rif) 



(4800.95 - 244.612 + 13/: 



2 

ms 



+ (-275.898 + 18.1861 - 0.666667 /^J Uf 
+ (4.97396 -0.185185 /, 



m 



2 \ 2 
Q \ 



a. 



ms) "^f) 



0.474894 - L 



TT 



+ — 4.59784 - 22.8619 l^s + 2 It 



+ (0.196497 + 0.88052 Ims) rif) + 



(4.77) 



with Ims = In / s. 

In Eqs. ( 4.76 ) and ( 4.77 ) the MS quark mass has been chosen as a parameter. This is 
inherent to the method used for the computation. Actually also the mass which is present 
in the renormalized operators O2 and Oq is defined in the modified minimal subtraction 
scheme. In order to transform the expressions into the on-shell scheme the two- ||123|| and 
three-loop ||124| , |125| , |126|| (see also Section relation between the masses is necessary^^. 



We obtain for the quadratic terms (see also Ref. ||127 



9Q,2 



M2 



Cts 



TT 



a.. 



12 + — (94.5000 + 24L^, - 4.33333^/) + 



TT 



(347.168 + 378L„ 



QLms + rif 



-67.6190 - 17.3333L„, + 2Li, ) + 1.21819 nj) + . . . 



M2 



as_ 

TT 



7.87659 + 0.35007 n/ 



+ ... , 



(4.78) 



where Lms = ^"^{Mq/s). Note the presence of mass logarithms which are introduced via 
the transition to the on-shell scheme. 

The quartic terms read in the on-shell scheme 



,os 



'M, 



2 \ 2 



Ots 



TT 



a. 



6 + — (10 - 24L^,) + — 570.519 - 155.5L 



-26.8336 + <^Lms - 2L^ 



+ 



77 

vr 



15L 



(9157.82 - 444.899L„ 



- 147.750L^, + 7.50000L;^, + Uf (-936.140 + 243.009L, 

- 0.666667L^J + n) (20.6385 - 7.86797L^, + 1.44444L 

- 0.222222LL 



26.1667L: 



2 

ms 



^^Due to the absence of a Born term in Eq. (|4.76| ) the two-loop relation between the MS and on-shell 
quark mass is sufScient in this case. 
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^os 



- 0.474894 -Lms + — (9.79728 - 21.4282L„, - 2L^, 



\ S I \ TT 

+ Uf (0.196497 + 0.880520L™,)) . . . 



(4.79) 



Note that in this case the transition to the on-shell scheme even introduces cubic mass- 
logarithms. 

In Fig. [4.9| rc,4, r^^i and rt^i are shown as a function of the center-of-mass energy, ^/i, 
where successively higher orders in are taken into account. It can be seen that the 
major part of the result is given by the Born approximation. The correction terms of 
order a^, a1 and are significantly smaller than the leading terms. However, it can be 
observed that with increasing order they remain roughly comparable in magnitude which 
could indicate a bad behaviour of the perturbative expansion. Nevertheless, the higher 
orders are small compared to the m^/s^ Born terms. It was shown in | |121[ | that the 
overall prediction for R{s) is stable and a variation of /i between ^/s/2 and for rc,4 
at 6 GeV varies by ±0.0005 and for rfe 4 at 14 GeV by ±0.0016. In the case of 4 the 
variation is negligible. 

Thus a prediction for Rq{s) up to order is available. It includes mass terms in an 
expansion up to the quartic order. 



5 Asymptotic expansion and Pade approximation 

Using the techniques currently available, diagrams beyond two loops can only be computed 
for special cases. In particular, only diagrams depending on one scale, like an external 
momentum or an internal mass, have been studied systematically. This section is devoted 
to a method which allows for the numerical reconstruction of a function depending on two 
dimensionful parameters like an external momentum q and a mass M. It constitutes a 
powerful combination of asymptotic expansion and the analytic structure of the function 
to be approximated. 

The basic idea is as follows. Let us consider a Feynman diagram depending on an 
external momentum q and one mass parameter, M, which can occur in some of the 
internal lines. Then the final result is a function of M^/g^, which in general is quite 
involved and at three-loop order — at least with the current techniques — not computable 
in an analytical form. On the other hand, it is straightforward^ to evaluate the diagram 
in the limits <C and 3> M^. The information from the different kinematical 
regions is combined and a semi-numerical function of M^/g^ is constructed. Below we 
will demonstrate on typical examples that it provides a very good approximation to the 
exact result. 



In Subsection ^rT| we start with a detailed description of the method and present 



explicit results for a two-loop example. The physical processes discussed afterwards in 

^^With "straightforward" we mean that program packages exist which allow the computation of the 
corr espo nding expressions with the help of computers. In this context see also the Appendices A. 2 



and A. 3 
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Figure 4.9: Quartic mass corrections (oc m^) to the non-singlet contribution of Tc (a), r^, 
(b), and (c) arising from diagrams where the external current couples directly to the 
massive quark. 



Sections and ^]4| point out different aspects and fields of application. In the 

first application we consider current correlators at three loops both for a diagonal and a 
non-diagonal coupling to quarks. In the latter case we assume that one of the quarks is 
massless. The correlators depend on the external momentum q and the quark mass M. 
The main interest is in the imaginary part which represents a physical observable. 

Also the two-point function considered in Subsection — proper combinations of the 
quark selfenergy — depends on an external momentum and one mass parameter. However, 
the main interest is not on the functional behaviour but on the value at threshold, i.e. 

= M^. The successful application of our method in that case is not obvious, especially 
as only the information for <^ and -C are incorporated. 

In Subsection four-loop integrals are computed in order to obtain the order QED 
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corrections to the muon decay. As we are only interested in the imaginary part they can 
be reduced to three-loop integrals using asymptotic expansions. From the technical point 
of view we want to obtain the value of a function for = using only expansion terms 
for <C M^. Also here the Fade approximation turns out to be quite successful. We 
want to mention already here that the example of the muon decay does not fit completely 



into the philosophy which is developed in Section |5.1j as no high-energy information can 
be incorporated. This leads to worse approximations. Furthermore, in a first step the 
imaginary part is taken and afterwards the Fade approximation is applied. 

5.1 The Method 

The basic ingredients and tools of our method are moments of the function to be con- 
sidered, conformal mapping and Fade approximation. In this Section we will explain the 
role of each of them. 

The aim is to obtain an approximation to the function f{z) which can not be com- 
puted directly. The analytical properties of f{z) are exploited and expansions of f{z) for 
small and large argument are used in combination with a conformal mapping and Fade 
approximation. In what follows we have in mind the computation of approximations to 
a physical function, like a vector boson self energy, at a given loop-order[^. In particular, 
we assume that f{z) is analytical for 2; 0. However, in the limit l/z ^ we allow for 
a non-analytical behaviour. 

In general a Fade approximation of a function f{z) is defined through 

K-]w = • • . (5-1) 

l + biz + ... bmZ"' 

where the coefficients and bj are determined from the requirement that the Taylor 
expansion of Eq. coincides with the first n + 171 + 1 terms of the Taylor expansion of 
f{z) around^ z = 0. Thus, in case the Taylor expansion is known up to terms of order 
z^ Fade approximations [n/m] fulfilling the condition 

k-1 > n + m, (5.2) 

can be computed. 

In the approach discussed in the remaining part of this Section the considerations of 
the previous paragraph are improved with respect to two points. First, we perform the 
Fade approximation not only in z but also in a new variable which is confined to the 
interior of the unit circle and thus provides better convergence properties. Second, we 
want to include in the approximation for f{z) not only Taylor coefficients around z = 
but also information from other kinematical regions, in particular from z — > 00. 



^°This is in contrast to the considerations of Ref. [128| where the Pade approximation has been per- 
iled in the coupUng cc 
'^In general also Tayk 
choice zq = is sufficient 



formed in the coupling constant in order to estimate higher order terms in as- 

^^In general also Taylor expansion around zq ^ can be considered. However, for our purpose the 
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Our procedure is restricted to the approximation of functions which depend on one 
dimensionless variable, z. Concerning the physical applications we have in mind two- 
point functions where z is given by the ratio of the squared external momentum and the 
square of an internal mass. This motivates the following form of /(z) for small and large 
argumentQ 



^ cuz^ for 2: 

nial . ) (5-3) 

XI Xl'^fc,* for 2 ^ -00 



L fe=0 i>0 



where the normalization /(O) = has been chosen. As we will see below, in the compu- 
tation of the coefficients and d^^^i — in the following also refered to as moments — one 
of the scales drops out and the integrals to be evaluated are much simpler. The moments 
Cfc and (ifc^j will serve as input for our procedure. In addition we also admit information 
about the behaviour at the physical threshold, which we choose to be at z = 1, as input. 

In order to obtain a semi-numerical approximation of the function /(z) the following 
steps have to be performed: 

1. Compute as many moments as possible for small and large z. As we require an- 
alyticity for 2; — > one gets in this limit a simple Taylor series of the Feynman 
diagrams in the external momentum. The expansion in q can be performed before 
the momentum integrations are performed. As a result the external momentum no 
longer appears in the integrand and one ends up with vacuum diagrams. They are 
analytically known up to three-loop order in case of one internal mass parameter 



(cf. Appendix |Al2|) . 



However, for z — > —00 the rules of asymptotic expansion ||129|| have to be applied. 
As a consequence the number of individual terms to be considered in the practical 
computation is larger. However, also here the number of scales in the individual 
diagrams is reduced. One ends up with either vacuum integrals or massless two- 
point functions. The latter are responsible for the logarithmic terms in z. 

2. Incorporate the information for z ^ 1 which we denote as (2;). In the physical 
examples considered below this information is either logarithmical, i.e. of the form 
\n{l — z), or proportional to Xj \J\ — z. The latter occurs, for instance, in the abelian 
contribution to the vector current correlator (cf. Section |5.2| ) and corresponds to 
the Coulomb singularity. The further steps slightly depend on which case is present. 

In case the leading threshold behaviour is logarithmic one constructs f^^^{z) in such 
a way that the singularity is reproduced for z — >■ 1. One has to take care that f^^^{z) 
does not destroy the behaviour of f{z) for small and large z. In particular, the 
expansion of f^^^{z) has to be analytical for z 0. By construction the difference 
f{z) — f^^^{z) is regular for z 1 and has the same limiting behaviour as the one 



22 



It is advantageous to consider the space-like region of z where no imaginary part occurs. 
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required in Eq. ( [5.3| ). We should mention already here that due to the construction 
of the Pade method the resulting function has a vanishing imaginary part at z = 1. 
Thus it is crucial to implement the leading threshold behaviour in this way. 

Threshold singularities of the form 1 / VI — z are not taken into account at this 



step, i.e. formally / 



thrf 



is chosen. They are treated below. In contrast 



to the logarithmic singularities they are removed via multiplication and not by 
subtractions. 



3. Construct a function f''"^{z) in such a way that the combination 



(5.4) 



is polynomial both in z and 1/z, i.e. in the small- and high-energy region. Further- 
more no logarithmic singularities may be introduced for z ^ 1. 

In this step a large part of information (e.g., the large high-energy logarithms) 
which is known analytically is extracted and only a small remainder f{z) is left. It 
parameterizes the unknown part of f{z). 



4. Perform a conformal mapping. The change of variables ||T3i 



(5.5) 



(1+CU)2' 

maps the z plane into the interior of the unit circle of the a; plane. Thereby the cut 



[1, oo) is mapped to the perimeter. The conformal mapping is visualized in Fig. [O . 
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Figure 5.1: The conformal mapping ( p. 51) maps the z plane into the interior of the unit 
circle in the lo plane. 



5. In a next step a function is defined for which finally the Pade approximation is 
performed. Due to the discussion in the context of Eq. ( p.l|) we are interested to 
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shift the information available for cu— Itocj^O. Furthermore we have to take 
care of possible power-like threshold singularities. 

Following Ref. [[L31|l we define 

(5.6) 

where for / the moments up to order must be known. The function p^^^{uj) is 

equal to 1 in case f{z) has logarithmic divergences at threshold and p*^'"(u;) = {1—uj) 
if f{z) is proportional to 1/Vl — z for 2 — 1. The available information from the 
moments transforms into P„;^^(— 1) and P^'^^ (0), {k = 0,1, . . . , riiar+ngrna — ^), where 
Usma is the number of moments for z ^ 0. Whereas for a logarithmic threshold 
behaviour the corresponding information is already taken into account in step 0, 
the — z behaviour is treated with the factor p^^^{uj) and in addition Prn^^{uj) 

is known for u = 1. 

6. In the last step a Fade approximation is performed for the function Pni^^{uj). This 
means that Pni^^{uj) is identified with a function [n/m]{uj) as defined in Eq. ( p. ID , 
where the number of coefficients on the right-hand side depends on the amount of 
information available for Pni^^{uj). In particular, one has n + m = niar + ^sma + 1 if 
P„;^^(l) is available and otherwise n + m = niar + Usma- This leads to a system of 
(non-linear) equations which can be solved for the coefficients a, and bj in Eq. ( |5.1|) . 
For large values of n + m the analytical solution becomes quite lengthy and time 
consuming. Thus it is preferable to solve the equations numerically (using high 
precision) . 



7. Finally, Eq. ( ^.6|) has to be solved for f(z) and from Eq. ( |5.4|) an approximation for 
the function f{z) is obtained. 



Due to the structure of Eq. ( |5.1| ) some Fade approximants develop poles inside the unit 
circle (|a;| < 1). In general we will discard such results as they would induce unphysical 
poles in the z-plane. In some cases, however, a pole coincides with a zero of the numerator 
up to several digits accuracy. These Fade approximations will be taken into account in 
constructing our results. If not stated otherwise we will, in addition to the Fade results 
without any poles inside the unit circle, also use the ones where the poles are accompanied 
by zeros within a circle of radius 0.01, and the distance between the pole and the physically 
relevant point q'^/M^ = 1 is larger than 0.1. 

There are situations where the information for z ^ —oo can not be used as this 
would lead to physically not allowed scenarios (cf. Subsection p.3| ). In this case it is not 
necessary (and even not possible) to define the function Pm^^^uj). Instead one can directly 
perform a Fade approximation either in z or in u. 

At this point we should spend some words on the estimation of the errors to be assigned 
to the final results. It is difficult to provide general rules for their determination as it very 
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much depends on the problem under consideration. Experience on the estimation of the 
error can be gained from the comparison with known results at lower order or for other 
colour structures. A reasonable choice for the systematic error due to the used method is 
to take the spread of the individual Pade results. 

Concerning the above list some comments to point |^ are in order. In principle there 
are many ways to subtract the high-energy logarithms. However, one has to keep in mind 
that the subtraction must not spoil the polynomial behaviour for small z. Furthermore, 
no divergences may be introduced. In particular, f''°^{z) has to be regular for z = 1. 

For the construction of f^°^{z) it is convenient to use the function 

GW = (5.7) 



with 



- i - 1 

u = \ 1 , (5.^ 

a - i + 1 



which naturally occurs in the result of the one-loop photon polarization function, as a 
building block. That this is possible in a systematic way can best be seen by looking at 
the expansion of G{z) in the different kinematical regions 



G{z) 



2 8 

1 + -Z+ —z^ + O (z^) z ^ , 

3 15 ^ ^ 

-ln(-l/4z) + — (l + ln(-l/4z)) + o(-) ^ ^ -oo , (5.9) 

TT 



From this equation one can see that f°^{z) can be chosen as a linear combination of 
terms {l/zy{l — z)^ z™'{G{z))^ {j,l,m,n > 0) where the corresponding coefficients are 
determined as follows: 



Consider the term in the second equality of Eq. (|5.3|) which has the lowest value of 
k and the highest value of z 7^ and fix the index n such that the powers of the 
logarithms coincide. 

Determine / in such a way that there is no singular behaviour for z ^ 1. 

As G{z) starts with order 1/z for z —00 one eventually has to correct for it with 
the help of the index m. 

In a similar way the index j is used in order to subtract the logarithms suppressed 
by higher powers in 1/z. 

Finally, terms involving / = m = n = and j > are added in order to restore the 
behaviour for z Q. 
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• Repeat the procedure with the next values of k and i, i.e. lower i by one unit until 
2 = is reached; then increase k by one unit. 

By construction this algorithm terminates once the linear logarithm of the largest high- 
energy moment is treated. 

To our knowledge a simple version of the method was first applied in [ [I32| ] for the 
evaluation of certain four-loop contributions to the anomalous magnetic moment of the 
muon. It was obtained by a convolution over the photon polarization function. For the 
latter an approximation formulae was obtained with the help of the Pade-method. One 
should stress that for this application only the integral over the space-like momenta of 
the approximation was used. In the applications which are discussed in this review the 
approximated function itself and in particular its analyticity properties are of interest. 
A brief introduction to the Pade-method and the discussion of some results can also be 



found in ||133|| . 



5.1.1 Explicit example at two loops 

For clarity let us present an example where in all steps explicit results are given. We 
consider two-loop QCD corrections to the correlator of two vector currents (for a precise 
definition see Eq. ( |5.15D below). 

1. The expansion of the diagrams for small external momentum leads to two- loop 



vacuum integrals which, e.g., can be computed with the help of MAT AD |^. After 
renormalizing the quark mass in the on-shell scheme and subtracting the constant, 
the result for the first three expansion terms reads 

TTm^,/ 3 /328 17962 999664 3 \ , 

' 167r2 V 81 675 496125 / ^ ^ 

In the high-energy region we restrict ourselves to the first two terms which can be 
obtained by solving massless integrals with the help of MINCER Using again the 



on-shell quark mass definition and taking into account the condition 11(0) = gives 

n(i)''^(z) = ^ - 4C3 - ln(-4z) - - ln(-4z) + ...). (5.11) 
167r^ \6 z / 

At threshold 11(2;) has a logarithmic singularity which can be cast in the form 

n«''''*'^^(^) = Aln(^) . (5.12) 

Thus, the combination Il{z) — Il^'^^{z) is constant for z = 1, has a polynomial 
behaviour for 2; ^ and at most logarithmic singularities for z ^ —00. 
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3. The high-energy logarithms are subtracted with the help of the function 



- 21 + 2 + 7r2(3 + 5z) 

(5.13) 



167r2 3z ^ 

+ 3(-l + 2) (-7 - 2^ + 7r2(l + 2^)) G{z 



which is constructed using the algorithm outlined above. The resulting function 
flW,v(^z) = n(^)''^(z) - n(^)''^'*'^''(z) - n(^)'^''°s(^) is analytical for z ^ and free of 
logarithms in the first two high-energy terms. 



After the conformal mapping (cf. Eq. ( |5.5| )) and the [2/2] Pade approximation in u 
are performed one obtains for 11^^^''" (z) 

, (0.874397 + 0.905702cj + 0.165184a;2w;L^ )2 

n(^)'''^ = -0.874397+^ 5.14 

^' 1 + 0.860764 CJ + 0.069525^2 ' 



which finally leads to li^amA^) = n(i)'^(2) + n(i)'^'*'^''(z) + Ii^^'^^''^^"9{z). In Eq. ( pTHj ) the 



numbers are truncated. Usually high numerical precision is needed in order not to loose 
significant digits in the final result. 

By construction, Yi^^^^^{z) has the same analyticity properties as the exact function. 
As in this case the latter is known one can also check the quality of the approximation. 
It turns out that even for the relatively small amount of input used here there is a 
perfect agreement between n^p^5!(z) and Ii^^^''"{z). E.g., it is not possible to distinguish 
the imaginary parts plotted in the range < 2m/y^ < 1 ||134|, |135||. 



5.2 Current correlators in QCD 

A variety of important observables can be described by the correlators of two currents. If 
the coupling of the currents to quarks is diagonal quantities like e^e~ annihilation into 
hadrons and the decay of the Z boson are covered by the vector and axial- vector current 
correlators. Total decay rates of CP even or CP odd Higgs bosons can be obtained by 
the scalar and pseudo-scalar current densities, respectively. For these cases the full mass 
dependence at order a1 has been computed in | 134 | for the non-singlet and in | 131 | for 
the singlet correlators. In Subsection |5.2.1| these results will be briefly reviewed. 

On the other hand the correlators involving different quarks describe, e.g., properties 
of a charged gauge or Higgs boson. In particular a certain (gauge invariant) class of 
corrections to the single-top-quark production via the process qq — > tb becomes available. 
As an application of the (pseudo-) scalar current correlator we want to mention the decay 
of a charged Higgs boson, which occurs in extensions of the SM, into a massive and a 
massless quark. Another important application of the non-diagonal current correlator 
is connected to the meson decay constant. Within the heavy quark effective QCD it is 
related to the spectral density, evaluated near threshold. The latter can be obtained from 
the correlator of the full theory which is considered below. 
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Let us in a first step introduce the polarization functions for the four cases of interest. 
The vector and axial-vector {6 = v, a) correlators are defined as 



{-q'9,u + q,qu)n\q') + q,qMiq') = z | dxe^^^(0|T4(x)jf (0)|0) , (5.15) 
and the scalar and pseudo-scalar ones (5 = s,p) read 

g2n^(g2) = iJdxe'''^{0\Tf{x)f^iO)\0). (5.16) 



The currents are given by 

j; = V^i7m^2, = ^iTmTsV's, / = ^^i^2, / = ^^V^i75^2. (5.17) 



Here m is the MS and M the on-shell quark mass. In Eqs. ( 5.15 ) and (|5.16|) two powers 
of q are factored out in order to end up with dimensionless quantities 11'^ (g^). As we are 
mainly interested in the imaginary part, the overall renormalization can be performed in 
such a way that this is possible. Furthermore it is advantageous to adopt the QED-like 
renormalization n'^(O) = 0. 

The physical observable R{s) is related to Il{q^) by 

R\s) = 127rlmU\q^ = s + ie) foTS = v,a, (5.18) 

R\s) = Svr ImU^q"^ = s + ie) ioTS = s,p. (5.19) 

It is convenient to define 



U\q') = U^<^)^\q') + ^^CFU^'^''{q')+(^^]\^'^^\q') + 

IT y TT y 

n(2).^ = n?)'^ + CaCfU^^Y + CpTmu'i^^'^ + CfTuP'^ + CfTU^'^ , (5.20) 

and similarly for R^{s). The abelian contribution II^'*''^ is already present in (quenched) 
QED and T^na originates from the non-abelian structure specific for QCD. The polariza- 
tion functions containing a second massless or massive quark loop are denoted by np"*'^ 
and n^-'''', respectively, n^^^'"^ represents the double-triangle contribution. 

Actually, we are mainly interested in the imaginary part R^{s) which in principle 
could be obtained from tree diagrams with five external legs, from one-loop four-point 
integrals and from two-loop three-point integrals. However, in particular the latter can 
not be evaluated analytically using current methods. Also numerically the treatment of 
these integrals is inconvenient. On the other hand, if one has to rely on approximations 
like small or large external momenta it is much more advantageous to stick to two-point 
functions simply because the resulting integrals are easier to solve and the corresponding 
techniques are much more advanced. 
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5.2.1 Diagonal current correlators 



The diagonal correlators with all their applications have extensively been discussed in the 
literature. Therefore we will be brief in this subsection and mainly refer to the original 
literature where in most cases also the analytical results can be found. 
It is useful to define dimensionless variables 



Aw? ' 



X 



2m 

7^' 



(5.21) 



where q is the external momentum of the polarization function and s corresponds to the 
center-of-mass energy in the process e~^e~ hadrons or the mass of the boson in case 
decay processes are considered. Then the velocity, v, of one of the produced quarks reads 



(5.22) 



Every time the generic index 5 appears without further explanation it is understood that 
8 represents one of the letters a, v, s or p. 

Vector and axial-vector correlators. The vector correlator certainly plays the 
most important role, mostly because it covers the processes induced by the photon. Al- 
ready in 1979 the massless a] corrections have been evaluated [30| and roughly ten years 
later even the order corrections became available BH]. However, due to the impres- 



sive experimental precision the massless approximations are not sufficient for a reliable 
comparison as we will see below in the case of the hadronic contribution of Aa (cf. Sec- 
tion |5.2.2| ). However, a complete analytical computation of W{q^) at three- loop order 
is currently not feasible. In Ref. ||134|, 135 | the Fade method described above has been 



applied and semi-numerical results have been obtained. At that time only the mass cor- 
rections of order im? /q^ [[114|| , for which no asymptotic expansion has to be applied, have 
been available in the high-energy region. They have been combined with the terms up to 
order (g^/m^)^ ||135|| in the small-g^ expansion in order to get semi-analytical results for 
the individual colour structures. 

As the main interest is in the imaginary part, one could also adopt the attitude 
to compute as many terms in the high-energy expansion as possible. Going, however, 
beyond the m? terms a naive expansion fails and the large-momentum procedure has to 
be applied. The calculation becomes very cumbersome if it has to be applied by hand. 
For this reason the large-momentum procedure has been automated and the program 
Imp l?^ has been developed. As a first application correction terms up to order {m? /q^Y 
have been evaluated [ 136 , 137|| for the vector correlator. In Fig. the comparison of the 



individual expansion terms with the Fade result from Ref. [134 is shown. 



For all three functions R^^'"" ^ ^na" ^'^d rI'^^''" and values between a; = and x 



0.6 

{x = 2m/ Y^) the expansions including terms of order {m? / sY (or more) are in perfect 
agreement with the semi- analytical Fade result. Conversely this provides a completely 
independent test of the method of [|134| [ which did rely mainly on low energy information. 
Including more terms in the expansion, one obtains an improved approximation even in 
the low energy region. However, the quality of the "convergence" is significantly better for 
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Figure 5.2: The abelian contribution R^'^, the non-abehan piece R\jA ■, the contribution 
from hght internal quark loops uf''^'^ and the contribution R^p^'"" from the double-bubble 
diagram with the heavy fermion in both the inner and outer loop as functions of a; = 
Imj Wide dots: no mass terms; dashed lines: including mass terms [m? / s)^ up to 
n = 5; solid line: including mass terms up to (m^/s)^; narrow dots: semi-analytical result 
(except for Rp^'^). The scale //^ = has been adopted. 



i^P"*'^ and R^^a" than for R^"^'" . Two reasons may be responsible for this difference: (i) In 
a high energy expansion it is presumably more difficult to approximate the 1/v Coulomb 
singularity in R^"^'^ than the mild Inf singularity in R^^^ Rf^'^ . (ii) The function 
Rf"*'" can be approximated in the whole energy region 2m < ^/s < oo by an increasing 
number of terms with arbitrary accuracy. This is evident from the known analytical 
form of |m| , a consequence of the absence of thresholds above 2m in this piece. In 

contrast the functions R^^'"" , R^na" Rp^'^ exhibit a four particle threshold at ^/s = 4m. 
The high energy expansion is, therefore, not expected to converge to the correct answer 
in the interval between 2m and 4m. In particular, in the case of 4')''' it can be seen that 
for X > 0.5 no convergence is observed. 

At this point it is tempting to combine both approaches — asymptotic expansion in the 
high-energy region and the Fade method — and evaluate high-order Fade approximants. 
A detailed study in the case of the vector correlator can be found in |^4| from which 
Fig. 15.31 is taken. In Fig. 15.31 the influence of the number of low- and high-energy input 
data is studied. In all four plots the quantity 6R)^ = R\ '^ — 3(7r''/(8f) — 37r^) is shown. 
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(c) (d) 

Figure 5.3: Abelian part of the vector correlator as a function of v. The leading threshold 
term is subtracted as described in the text. In (a) and (b) the number of low-energy 
moments are varied and the high-energy terms are fixed to two and seven, respectively. 
In (c) and (d) the low-energy moments are fixed to four and eight, respectively, and the 
high-energy terms are varied. The notation means that terms up to order 

and \j z^^^ are taken into account in order to construct the Fade approximation [i/j]. 



i.e. the leading threshold term is subtracted. Otherwise it would not be possible to detect 
any difference between the individual Fade results. Moreover, the abscissa only extends 
to w = 0.2 as for v ^ 0.5 all curves coincide. In Figs. |5.3| (a) and p.3| (b) the higher-energy 
terms are fixed to two and seven, respectively, whereas an increasing number of low-energy 
moments are considered. On the other hand, on the plots in the lower row the number of 
moments is fixed to four and eight, respectively, with varying high-energy input. 

From Figs. ^]B|(a) and p.3| (c) a clear stablilization of the results can be observed with 
increasing degree of the Fade approximation. The same is true for Figs. ^]3|(b) and p73|(d): 
although the degree of the Fade approximants is higher from the very beginning a further 
stabilization is visible. 

From these considerations one can conclude that both the small- and high-energy 
expansion terms are crucial as input for the Fade procedure. A significant stabilization of 
in the threshold region is observed if more terms are taken into account. However, 
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for practical purposes it is probably more than enough to consider, e.g., only the quadratic 
terms in the high-energy region. The reason for this is that also the leading threshold 
behaviour is incorporated into the analysis. The situation is different in those cases where 
no information about the threshold is available or one even wants to determine the value 
of the (real part of the) considered function for z = 1. We will come back to this point in 
Section |]|. 

At this point we refrain from listing the results of the individual terms contributing 
to R{s) as all of them are available in the literature. In particular we want to refer to the 
Appendix of Ref. [p.07|| where detailed results up to order a^m^/s^ and a^rn^/s are listed. 
Concerning the full mass dependence at order «^ a complete discussion and a detailed 
compilation of the individual terms in the on-shell scheme can be found in Ref. |p.27|| (see 
also ||138|| ). Together with the quartic on-shell terms at order given in Eq. ( ^4.79| ) this 
constitutes the current state-of-the-art radiative corrections for R{s). 

In ||135| next to the vector case also the axial- vector, scalar and pseudo-scalar correla- 
tors have been considered. Moments up to order (g^)* have been combined with quadratic 
mass terms and the leading threshold behaviour in order to obtain semi-numerical ap- 
proximations for the correlators. 

QCD corrections to the axial- vector correlator have also been evaluated in [|110|| . In 
the limit 3> the first seven terms could be evaluated using automated asymptotic 
expansion. In Fig. |5.4| the expansion terms are compared with the Fade results [|135|] for 
the individual colour factors. Again perfect agreement is found in the region where the 
asymptotic expansion is expected to converge to the exact result. 



In [ [I10| | the results have been used in order to obtain in combination with the vector 
correlator, order corrections to the top quark production in e~^e~ annihilation above 
the threshold. The results for the cross section are shown in Fig. p.b\ where also the 
electro- weak corrections have been included ||139||. 



Scalar and Pseudo-scalar correlators. For completeness we want to mention 
that the high-energy expansions of the scalar and pseudo-scalar correlators have been 
considered in The results have been applied to the decay of a scalar and pseudo-scalar 
Higgs boson into top quarks where higher order mass effects are important. This can be 
seen in Tab. |5.1| where a (pseudo-) scalar Higgs boson mass of 450 GeV has been considerd. 
Both the results for the individual mass-correction terms {M'^/Mj^y {i = 0, ... ,4) and 
their proper sum is listed up to order a'j.. 

Singlet contribution. A special feature of the diagonal current correlator is the 
occurence of so-called singlet diagrams as pictured in Fig. 5^. They are often also referred 
to as double-triangle diagrams as the external currents are not connected through the 
same fermion line. Note that for the vector correlator there are no singlet diagrams at 
three-loop level according to Furry's theorem [p.4CI|| . In case of axial- vector couplings one 
has to take both members of a weak isospin doublet into account in order for the axial 
anomaly to cancel. It is therefore convenient to replace the current in Eq. (|5.15| ) by 
js,fi ~ '^1til5'^~XliJ.'l5Xy where ip and x are isospin partners. The diagrams contributing to 
the singlet part, 11^^^'" (g^), are depicted in Fig. |5.6| where in the fermion triangles either 
ip or X niay be present. Note that for a degenerate quark doublet Ilg'^'°'{q'^) vanishes. 
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Figure 5.4: Rr'"", i = A, NA, I, F, S , Sb as functions of x = IMtj^l at /i^ = M^. 
Successively higher order terms in {Mf / sY are included: Dotted: n = 0; dashed: n = 
1, . . . ,5; solid: n = 6. Narrow dots: exact result -R^^'") or semi-analytical results 

(i?^^'", R^^^J^). In the case of the top-bottom doublet Rsb contains only the imaginary 
parts of the singlet diagrams which arise from the gluon and bottom quark cuts (see also 
Fig. O). 



Having furthermore in mind the physical case (■?/', x) = {t, b), it is justified to set = M 
and = 0. In the case of the scalar and pseudo-scalar currents there is no anomaly. 
Furthermore, as the coupling is proportional to the quark mass, only one diagram has to 
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Figure 5.5: The electroweak (a) and QCD (b) corrections to e'^e tt. In this Figure 
the normalized cross section Rt = a{e~^e~ — > tt)/apt is shown, apt = Aira'^/Ss. 





Figure 5.6: Singlet or double-triangle diagrams. In the fermion lines either the quark 
or its isospin partner x may be present. 



be considered, namely the one where in both quark lines the same heavy quark flavour is 
present. 

As mentioned in Section it is essential that the expansion for — is analytical. 
However, in the case of the singlet diagrams this is not fulfllled as there exist massless 



cuts containing gluons and light fermions. Thus the method of Section 5.1 cannot directly 



be applied to n^^'''^(g^). Rather it is applied to 



n 



S,modiQ 
S.mod 



127r2 
1 



K 



Sb \ 
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dr- 
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(5.23) 
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(M2)0 










S 


exact 




scalar, on-shell 








1.000 


-0.907 


0.137 


0.014 


0.003 


0.247 


0.248 




-0.778 


5.646 


-3.080 


0.010 


0.004 


1.802 


1.802 




-35.803 


77.056 


-17.792 


-5.347 


-0.680 


17.435 


- 


S^(as/7r)' 


0.943 


-0.663 


0.026 


0.009 


0.003 


0.318 






scalar, MS 








0.784 


-0.558 


0.066 


0.005 


0.001 


0.298 


0.299 




4.445 


-3.723 


-0.238 


0.143 


0.032 


0.659 


0.673 




21.799 


-7.606 


-15.334 


0.680 


0.546 


0.086 


— 


S.(as/7r)' 


0.941 


-0.679 


0.045 


0.010 


0.002 


0.319 






pseudo-scalar, on-shell 






i?(0)'P/3 


1.000 


-0.302 


-0.046 


-0.014 


-0.005 


0.633 


0.629 




-0.778 


3.495 


0.417 


0.047 


0.027 


3.208 


3.238 


i?(2).P/3 


-35.803 


25.024 


12.173 


3.780 


1.293 


6.467 


- 


S.(as/7r)' 


0.943 


-0.172 


-0.022 


-0.009 


-0.003 


0.737 






pseudo-scalar, MS 






^i?(°)'V3 


0.784 


-0.186 


-0.022 


-0.005 


-0.002 


0.569 


0.569 


^CFi?(i)^P/3 


4.445 


-0.248 


-0.215 


-0.119 


-0.043 


3.821 


3.791 


W n(2),p /o 


21.799 


9.854 


2.455 


-0.781 


-0.520 


32.807 




S.(a./7r)' 


0.941 


-0.185 


-0.026 


-0.010 


-0.003 


0.717 





Table 5.1: Numerical results for i?* and W both in the on-shell and MS schemes. The 
contributions from the mass terms {MfY, their sum (E) and, where available, the exact 
results are shown. I]i{as/7ry is the sum of the 1-, 2- and 3-loop terms. The numbers 
correspond to Mt = 175 GeV and M^jj^ = 450 GeV. The renormalization scale fi^ is set 
to s = Mjjf^. 



where -R^^'" and R^-?}''^ denote the contribution of these massless cuts to and 



gg 



respectively. Thus by definition, n^^|^od(5'^) has the same analytical properties as the 
non-singlet polarization functions. The notation already suggests that in the scalar and 
pseudo-scalar case there is only the cut through the two gluons. On the contrary, this 
cut is zero in the axial-vector case according to the Landau- Yang-Theorem ||141|| and only 



cuts involving the massless quark contribute. The analytical expressions for R^gl'"" and 
Rfg''^ can be found in ||142|| and ||131|| , respectively. Expansions of the former are also 
listed in 11131 



In 1131 1 the asymptotic expansion has been applied to the singlet diagrams and terms 
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up to {M"^ /q^Y the high-energy region have been evaluated. For small external momen- 
tum moments up to order (axial-vector) and (scalar and pseudo-scalar) are available. 
This input was used to compute roughly 30 Fade approximants for each correlator. In 

(sohd lines), 



r(2)v 



and lif^" 



Fig. |5.7| (a)-(c) the results for the imaginary part of Ilg , 11 
together with the first seven terms of the high energy expansion (dashed and dotted lines) 
are shown as functions of x = Imj Note, that in the displayed region, < a; < 1, 



imn 



(2) 



Imn 



(2) 

S.mod' 



Therefore, if one is interested, e.g., in (inclusive) production of the 
heavy quarks only, the corresponding massless cuts (depicted as dash-dotted lines) have 
to be subtracted. The resulting curves are shown in Fig. |5]^ (d). 
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Figure 5.7: 



The imaginary parts i?!^"*'" 



and of 



(a) the axial-vector, (b) 



d(2)„ 
S ' -"-5 

scalar and (c) pseudo-scalar singlet diagrams, respectively. For the abscissa the variable 
X = 2m/y/s is chosen. Solid line: Fade result; wide dots, dashes and narrow dots: 
(m^/g^) "-expansion for n = 0, n = 1, . . . , 5 and n = 6, respectively; dash-dotted line: 
purely massless cuts R^st"^, RfJ'^ and RfJ'^- Fig. (d) shows the difference between the 
solid and the dash-dotted line (i.e., the contribution of the massive quarks) of Figs, (a), 
(b) and (c) as solid, dashed and dotted line, respectively. The curves in (a) can also be 
found in the last row of Fig. E^. 



The requirement that no logarithmic terms may appear in the expansion for z — > 
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(cf. Eq. ( p.3| )) seems to be quite restrictive. However, the example of the singlet diagrams 
shows that also problems which at first sight do not match the definition of f{z) in 
Eq. ( |5.3| ) can be treated. 

From the systematic evaluation of the polarization function at low and high energies 
and the information about the leading threshold behaviour it is possible to construct 
— with the help of Pade approximation and conformal mapping — an approximation 
to n(-2) taking into account the complete quark mass dependence. The imaginary part 
immediately leads to important physical quantities like R{s), top quark pair production 
above threshold and the total decay rate of scalar and pseudo-scalar Higgs bosons. 



5.2.2 An application of R{s): Aahad(^z) 



An important application of the vector current correlator, as discussed in Sections ^ 
and |5.2.1| , is the evaluation of the fine-structure constant at the scale of the Z boson, 
a(M|). It plays a crucial role in the indirect determination of the Higgs boson mass which 
constitutes one of the most important goals of precision experiments. In particular, the 
error on a(M|) induces one of the largest uncertainties. 

In Ref. |p.43|| a conservative analysis has been performed which exclusively relies on 
data below a center-of-mass energy, ^/s, of 40 GeV. Consequently it suffers from sizable 
experimental errors. Only for i/s > 40 GeV perturbative QCD has been used. Recently 
several suggestions have been made, which significantly reduce the uncertainties in a(M|). 
Most of them were actually triggered from the knowledge of the complete mass dependence 
at order (see Section ^.2|) . There is a number of so-called theory-driven analyses ||144| , 



145| , |146| , |147] , |148| , p.49| , |150|| which replace to smaller and larger extent unprecise data by 



«(^) = : — : . . / (5) . . . (5-24) 



the results of perturbative QCD (pQCD). They will be discussed below. 

Let us start with some basic definitions. Re-summation of the leading logarithms leads 

to 

a(0) 

X - Aaiep(s) - A«[Jd(s) - A«topi 

with a = a{0) = 1/137.0359895. Aaiep denotes the leptonic contribution and is known 
up to the three-loop order |p.51| ] 

Aaiep(M|) = 314.97686 X 10"^ (5.25) 

Perturbation theory is also applicable to treat the contribution arising from the top quark. 
Including three-loop QCD corrections one gets ||145|| 



Aatop(M|) = (-0.70 ± 0.05) X 10"^ (5.26) 

The contribution from the remaining five quarks has to be taken into account using the 
dispersion integral 



67r J Ami s[s — ik/| — le) 
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a-i 










273.2(4.2) 


128 


985(58) 


280(7) 


128 


892(95) 


280(7) 


128 


892(95) 


275.2(4.6) 


128 


958(63) 


281.7(6.2) 


128 


869(85) 


278.4(2.6)* 


128 


914(35) 


277.5(1.7) 


128 


927(23) 


277.6(4.1) 


128 


925(56) 


277.3(2.0)** 


128 


929(27) 


277.0(1.6)* 


128 


933(22) 


277.8(2.5) 


128 


922(34) 


274.3(1.9) 


128 


970(26) 


276.1(3.6) 


128 


946(49) 


277.3(2.1) 


128 


930(29) 



Reference 



Pi 



PI 



HI 



[148] 



|4g 



PI 



Martin et al. '95 
Eidelman et al. '95 
Burkhardt et al. '95 
Swartz '96 
Alemany et al. '97 
Davier et al. '97 
Kiihn et al. '98 
Groote et al. '98 
Erler '98 
Davier et al. '98 
Jegerlehner '99 
Martin et al. '00 
Burkhardt et al. '01 
Jegerlehner '01 



comment 



(low-order) pQCD 

data 

data 

r data 
+ pQCD 
+ CT 

SR (pQCD) 

r data + UDR 

r data + pQCD + SR 

pQCD (ER) 

new data + pQCD (ER) 
new data 

pQCD (ER) + new data 



Table 5.2: Comparison of the different evaluations of Aq;[^'|j(M|). The column "com- 
ment" reminds on the different methods used in the analysis as described in the text. 
(*Aatop(M|) subtracted; ** value corresponding to as[M\) = 0.118 adopted.) 



with R{s) = a{e^e —>■ hadrons)/cr(e"'"e —>■ fi~^fi ) (cf. Eq. (|5.19| )). It is not possible 



to use perturbation theory for R{s) in the whole energy region. Thus one has to rely to 
some extent on experimental results. The results of the recent evaluations can be found 
in Tab. ^.2| where also the resulting values for a~^(M|) are listed. The latter has been 



obtained with the help of Eqs. ( |5.24| ), ( |5.25| ) and ( |5.26| ). In the fourth column of Tab. ^ 



some keywords are given which shall indicate the method used for the analysis. In the 
following we briefly describe the abbreviations. For more details we refer to the original 
papers. 

T data | |155|| . r data from ALEPH have been used in order to get more information 
about R{s) for energies below roughly 1.8 GeV. The hypothesis of conserved vector current 
(CVC) in combination with isospin invariance relates, e.g., the vector part of the two-pion 
T spectral function to the corresponding part of the isovector e~^e~ cross section through 
the following relation 



a^=^ (e+e- ^ TT+TT-) = ^^^^^vj=, {r ^ nn%^) . (5.28) 

s 

A similar equation holds for the four-pion final state. Their incorporation into the analysis 
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has been performed in ||155|| leading to a slight reduction of the error on Aq;|^^^^[^. 

Perturbative QCD (pQCD). The first attempt to replace unprecise data by pQCD 
can be found in | 152 |. At that time, however, mass effects were barely known. Thus 
pQCD could only be applied far above the particle thresholds. Meanwhile R{s) can be 
calculated in the framework of pQCD up to order if quark masses are neglected |3^, Q 
and up to 0{a'i) with full quark mass dependence |T3|, |13|, |123. In |144|, |14| pQCD 
has been used down to an energy scale of ^/s = 1.8 GeV and it has been shown that the 
non perturbative contributions are small. This leads to a further reduction of the error 
of about a factor two. For convenience we list in Tab. the perturbative hadronic 



contributions for a variety of energy intervals. As our default values we adopt /i 



a(5)(M|) = 0.118 [|T5|], M, = 1.6 GeV and = 4.7 GeV. In separate columns we list 
the variations with a change in the renormalization scale, the strong coupling constant 
and the quark masses: 



6a. 



±0.003, 5Mc = ±0.2 GeV, 6Mb = ±0.3 GeV. 



(5.29) 



The typical contributions which have to be taken into account look as follows [|145|| . In 
the perturbative regions one receives contributions from light [u, d and s) quarks whose 
masses are neglected throughout, and from massive quarks which demand a more refined 
treatment. Below the charm threshold the light quark contributions are evaluated in order 

plus terms of order a1s/{AMl) from virtual massive quark loops. Above 5 GeV the 
full Mc dependence is taken into account up to order a^, and in addition the dominant 
cubic terms in the strong coupling are incorporated, as well as the corrections from virtual 
bottom quark loops of order a1 s/(4M^). Above 11.2 GeV the same formalism is applied to 
the massive bottom quarks and charm quark mass effects are taken into account through 
their leading contributions in an / s expansion. All formulae are available for arbitrary 
renormalization scale /i which allows to test the scale dependence of the final answer. 
This was used to estimate the theoretical uncertainties from uncalculated higher orders. 
Matching of as between the treatment with n/ = 3, 4 and 5 flavours is performed at the 
respective threshold values (cf. Section^). The influence of the small 0{a1) singlet piece 
which prevents a clear separation of contributions from different quark species can safely 
be ignored for the present purpose. 

Charm threshold region (CT) ||145|| . Perturbative QCD is clearly inapplicable in 
the charm threshold region between 3.7 and 5 GeV where rapid variations of the cross 
section are observed. Data have been taken more than 15 years ago by the PLUTO ||159|| , 
DASP | |160| |, and MARK I collaborations [ |161[ ]. The systematic errors of 10 to 20 % exceed 



the statistical errors significantly and are reflected in a sizeable spread of the experimental 
results. In | |145[| the experimental data are normalized to match the predictions of per- 
turbative QCD both below 3.7 and above 5.0 GeV. Two models have been constructed 
which describe the differences of the normalization factors below and above the considered 
energy interval. 



^■^On the contrary, the inclusion of the r data leads to a significant reduction of the error of the 
anomalous magnetic moment of the muon as it is more sensitive to the low-energy region. 
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Energy range (GeV) 


central value 


6fi 


5as 


SM, 


dMh 


1.800 - 2.125 


5.67 


0.22 


0.04 


0.00 


0.00 


2.125 - 3.000 


11.66 


0.21 


0.06 


0.01 


0.00 


3.000 - 3.700 


7.03 


0.06 


0.03 


0.00 


0.00 


1.800 - 3.700 


24.36 


0.48 


0.13 


0.01 


0.01 


5.000 - 5.500 


5.44 


0.03 


0.03 


0.06 


0.00 


5.500 - 6.000 


4.93 


0.03 


0.02 


0.04 


0.00 


6.000 - 9.460 


25.45 


0.11 


0.08 


0.10 


0.00 


9.460 - 10.520 


5.90 


0.02 


0.01 


0.01 


0.00 


10.520 - 11.200 


3.48 


0.01 


0.01 


0.00 


0.00 


5.000 - 11.200 


45.20 


0.19 


0.15 


0.21 


0.01 


(without bb) 












11.200 - 11.500 


1.63 


0.00 


0.01 


0.00 


0.00 


11.500 - 12.000 


2.62 


0.00 


0.01 


0.00 


0.00 


12.000 - 13.000 


4.93 


0.01 


0.01 


0.00 


0.00 


13.000 - 40.000 


72.92 


0.08 


0.12 


0.02 


0.02 


12.000 - 40.000 


77.85 


0.09 


0.14 


0.02 


0.02 


40.000 - CO 


42.67 


0.03 


0.06 


0.00 


0.00 


11.200 -00 


124.77 


0.12 


0.21 


0.03 


0.02 


1.8 - Qo (pQCD) 


194.33 


0.79 


0.49 


0.24 


0.03 


QED 


0.11 











Table 5.3: Contributions to Aa|^2i(^I) (i^ units of 10^^) from the energy regions where 
pQCD is used (adopted from ||127|| ). For the QED corrections the same intervals have 
been chosen. For the variation of as(M|), Mc and Mb (cf. Eqs. (|5.29|) ) have been used. 
H has been varied between and 2^/s. 



Recently, the BES collaboration has measured R{s) in the energy range between 2 and 

162|| . We applied the method of ||145|| and 



5 GeV with substantially improved precision 
obtained perfect agreement with the results of PLUTO, DASP and MARK I. 



QCD sum rules (SR) |4|, |48 



Global parton-hadron duality is used in order 
to reduce the influence of the data in the different intervals. This is achieved by choos- 
ing a proper polynomial, Qn{s), which is supposed to approximate the weight function 
Mz/s{s — M|) in Eq. ( ^.271) as good as possible. Adding and subtracting Q^rl 



m 



Eq. ( p.27| ) and exploiting the analycity of the subtracted term leads to 



^1 , R{s)M} 



so 



sis -Mi: 



dsR{s) 



Ml 



s is 



Ml] 



Qn{s) 



+Qm<f dsU^^^{s)QN{s). 

^151=51 



(5.30) 



Thus the influence of the experimental data is significantly reduced in the first term of 
the right-hand side and pQCD only has to be used for |s| = si which is indicated by the 
superscript QCD. 
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Unsubtracted dispersion relations (UDR) [ |1471| . They are used in order to eval- 
uate the electromagnetic coupling in the MS scheme. Four-loop running is accompanied 
by three-loop matching in order to arrive at a(M|), which subsequently has to be trans- 
formed to the on-shell quantity a(M|). Via this method no complications in connection 
with the J/\E' or T resonances occur. However, one encounters a much stronger depen- 
dence on the quark masses. 

Perturbative QCD in Euclidian region (ER) jliOi [l50i [15^ . The authors of 



Ref. ||143|| also re-evaluated /S.a^l^{Mz) using pQCD. In a first step Aa^^j was calculated 
at the large negative scale s = — M| and then analytically continued to s = M|. Thus 
pQCD has only been applied in the Eulidian region where it is supposed to work best as 
one is far from resonances and thresholds. Furthermore pQCD has been used down to 
-(2.5 GeV)2. 

New data in the low-energy region. In the meantime new experimental data for 
R{s) in the low-energy region have become available. Besides improvements in the energy 
interval below 1.4 GeV by the CMD-2 detector at the VEPP-2M collider in Novosibirsk 
a measurement of R{s) in the range 2 GeV < ^/s < 5 GeV has been performed by 
the experiment BES 11 at Beijing. In Ref. [rST] these data have been incorporated and 



accompanied with pQCD in the regions 3 GeV < < 3.74 GeV and ^/s > 5 GeV in 
order to evaluate Aa^^^{Mz) ■ The result shown in Tab. 5^ has been obtained using only 



inclusive measurements of R{s) for ^/s <, 1.9 GeV. Based on the comparision of time-like 
and space-like (i.e. in the Euclidian region) evaluations of Aa^^^{Mz) it has been argued 
in 



150|| that this is preferred to the exclusive measurements of R{s). 



More recently the data-based analysis of ||153|| has been updated ||156|| using pQCD 



only above ^/s = 12 GeV. The main improvements are due to the new BES measurements. 



Tab. |5.2| shows that the inclusion of pQCD leads to a significant reduction of the error 
in Aa[|^^^(M|). The new analyses (with the exception of ||15CI|| p^) agree well both in their 
central values and in their quoted errors. These promising developments suggest to use 
the new values in the interpretation of the electroweak measurements. Once more precise 
experimental results on R are available it can replace the corresponding parts in the 
theory-motivated analyses. Certainly these measurements would be extremely valuable 
for a cross check of the theory-driven results. 

5.2.3 Heavy-light current correlators 



The method of Section |5.1| has also been applied to the non-diagonal correlators where one 
of the quark fields in Eq. ( |5.17| ) has mass M and the other one is massless. In this limit 



the vector (scalar) and axial-vector (pseudo-scalar) correlators coincide. Furthermore it 
is convenient to work with the variables 

^ = 1^' a; = ^, V = ^— ^. (5.31) 



24 Actually, in [|T50| also a result based on exclusive data is given (Aal^'^^dl*^!) = (276.49± 2.14) x 10"'*) 



which is in better agreement with the other values. 
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For the application of the Fade method the polarization function 11^(2;) has to be 



considered in the different kinematical regions. In [|163| , |164|| seven terms for small and 
eight terms for large external momentum have been computed both for the vector and 
scalar correlators of Eqs. ( |5.15|) and ( ^.16| ), respectively. This means expansion coefficients 
up to order , respectively, 1/2;^ are available for the Fade procedure. 

The threshold behaviour constitutes the main difference as compared to the diagonal 
correlator discussed in the Section |5.2.1| . The Born result for of the diagonal vector 
and pseudo-scalar correlators are proportional to v for small velocities. At order 
i?*^^)'*^ approaches a constant and at order one either has a 1/u or a Inw behaviour — 
depending on the colour structures [[134|, [135|| . The axial- vetor and scalar correlators show 



the same pattern with a additional factor v at each order. 

On the contrary the imaginary part of the non-diagonal correlators are proportional 
to (possibly accompanied with Inf terms) independent of the order in and of the 
colour structure. This is valid in every order in as follows from Heavy Quark Effective 
Theory (HQET) [|165|| . Actually the latter can be used to obtain the leading threshold 
behaviour of R"{s) and R^{s) at 0{a'^) from the corresponding correlators in HQET. 
In particular, the renormalization group equation in the effective theory is used to get 
the leading logarithmic behaviour at order ||164|| . Afterwards the decoupling relation 
between the currents in the full and effective theory [ |166| , |167| , |168|| is exploited to get 
the information about R'"{s) and R'^{s) \\163[ |164|] . It turns out that linear and quadratic 
logarithms occur. This translates into quadratic and cubic logarithms of the corresponding 
polarization function which are incorporated into the Fade method as descibed in step |] 
of Section 

Before discussing the results in the full theory we want to spend time on the spectral 
function in HQET. In numerical form it reads jl64|| 



R'{cu) 



1 + (8.667 + L J + f ' (g^^ + 35.54L. 

TT \ TT J 



+ 1.875L2 + m (-1.583L^ - 0.08333L^ 



(5.32) 



where u = ^/s — M is the only dimensionful quantity in the effective theory and = 
ln(/iVcu2). The tilde and the prime remind that the quantity on the left-hand side of 
Eq. ( |5.32|) is defined in the effective theory where the heavy quark is decoupled ||167| , |168| . 
In Eq. ( 5.32| ) the corrections of order ag are known since quite some time [ 169 | whereas 



the constant c„, , which is not accessible using renormalization group techniques, has been 



determined in [|164|| with the help of the Fade results in the full theory. Its dependence 
on the number of massless quarks is given by 



46(15) - 1.2(4) ni. 



(5.33) 



In the meantime the coefficient c^^ has been computed analytically | |1 70| | with the result 
Cm = 58 — 1.7 n;. Agreement with Eq. ( |5.33|) can be observed for all physically reasonable 
values of rii. 
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We want to mention that R' enters as a building block into the sum rules which are 



used to determine the meson decay constants like, e.g., fs (see, e.g., Ref. ||169|| ). The 
typical scale where Eq. (|5.32| ) has to be evaluated is 1 GeV. Thus the first order QCD 
corrections amount to about 100% and the terms of order contribute with additional 
60(20)% where the sign is the same as for the LO correction. 



In Figs. |5.8| and ^.9| the imaginary parts of the individual three-loop colour structures 
for the vector and scalar correlators are plotted as a function of the variable x. Each solid 
line contains of the order of 15 Fade approximants which show perfect agreement among 
each other. For comparison also the high-energy expansion terms including order l/z' 
are shown as dashed curves. Excellent agreement with the semi-numerical Fade results is 
observed up to x ~ 0.5 which corresponds to f ~ 0.60. Some colour structures {Rf'^'"{s), 
R^2^'^{s), R^^^/{s), -Rp''''*(s)) show even an agreement up to a; ~ 0.7 {v ~ 0.34) which is 
already fairly close to threshold. We want to remind that the functions which exhibit next 
to the cut at ^/s = m also a cut at ^/s = 3m are not expected to converge to the correct 
answer in the interval between m and 3m. This explains the somewhat crazy behaviour 
of Rp^'^{s) for larger values of x. 

The non-diagonal correlators describe properties connected to the W boson. In par- 
ticular, a certain (gauge invariant) class of corrections to the Drell-Yan process, i.e. to 
the production of a quark pair through the decay of a virtual W boson generated in pp 
collisions, are covered by the vector and axial-vector correlator. The absorptive part is 
directly related to the decay width of the (highly virtual) W bosons into quark pairs and 
gluons. Of particular interest in this connection is the single-top-quark production via the 
process qq —>■ tb. The imaginary part of the transversal W boson polarization function 
constitutes a finite and gauge invariant contribution of 0{a^). 

The corrections of order as to the (total) single-top-quark production rate are quite 
large. They amount to about 54% and 50% for Tevatron and LHC energies, respec- 



tively [171 , where 18%, respectively, 17% arise from the final state corrections. This calls 
for a complete C(a^) calculation. 

If one considers the leading term of the large- A'^c limit it is possible to use the results 
for R" to perform a theoretical analysis of order to the single-top-quark production. 
The production cross-section of the virtual W* boson is identical to that of the Drell-Yan 
process qq eVf.- The latter is known to 0{a1) from Ref. | 172 |. Thus we can take proper 



ratios to make predictions in the large-iVc limit at NNLO free from any dependence on 
parton distribution functions. As an example, we consider 

jli{PP^W*^th) _ Imp^g^)] 



^{pp^W*^eu,) Im[ne.(g2)] 

= N,\Vtb\''R''is). (5.34) 



The numerical significance of the order corrections is shown in Fig. |5.10| where the 
LO, NLO and NNLO result of R''{s) is plotted in the range = 200 . . .400 GeV. For 
the numerical values Mt = 175 GeV and as{Mz) = 0.118 has been chosen. Whereas the 
0{as) corrections are significant there is only a moderate contribution from the order 
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Figure 5.8: R^2^''"{s), R^^^^{s), and Rp'^'^{s) as a function of x. The dashed curves 

correspond to the analytical expressions obtained via asymptotic expansion containing the 
terms up to order l/z"^. 



terms. In the range in shown in Fig. |5.1CI| they are below 1% of the Born result. Note 



that (at least for /i^ = M|) the NNLO corrections to the Drell-Yan process are also small 
and amount to at most a few percent (see e.g. ||173|| ). 

As an application of the scalar and pseudo-scalar current correlator we want to mention 
the decay of a charged Higgs boson which occurs in extensions of the Standard Model. 
The corrections to describe the total decay rate into a massive and a massless quark. 
To be more precise, the hadronic decay rate of the charged Higgs boson takes the form 

T{H+^UD) = ^^^MH+{a^ + b^)R'{Mfj+), (5.35) 

where a and b parametrize the coupling of the Higgs boson to the massive quark U and 
the massless quark D 

Jh+ = [a (1-75) + 6 (1 + 75)]^. (5.36) 



88 




Figure 5.9: R^^^'^{s), R^^'^/{s), and Rp^'^{s) as a function of x. The dashed curves 

correspond to the analytical expressions obtained via asymptotic expansion containing the 
terms up to order l/z"^. 



In Fig. |]TT| R^{s) is plotted at LO, NLO and NNLO |T6|, |T6J. Again it turns out 
that the radiative corrections are well under control as order terms contribute at most 
of the order of 1%. 



5.3 QED corrections to muon decay 

The Fermi coupling constant, Gp, constitutes together with the electromagnetic coupling 
constant and the mass of the Z boson the most precise input parameters of the SM of 
elementary particle physics. Gp is defined through the muon lifetime which in turn is 
obtained from the decay rate. The one-loop corrections of order a have been evaluated 
more than 40 years ago in Only recently the two- loop corrections of order have 
been computed by two independent groups fl^, |19] . The inclusion of the new terms leads 



to the following value for the Fermi coupling constant ||174|| 



Gf = 1.16639(1) X 10"^ GeV~^ (5.37) 
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Figure 5.10: LO (dotted), NLO (dashed) and NNLO (solid) results of i?^(s). 



where the error is reduced by a factor two. It is now entirely experimental. In this 
Subsection we want to concentrate on the method used in [1^ which is based on Fade 
approximation and conformal mapping. We also want to mention that in Ref. | |175| | 
optimization methods have been used in order to estimate the coefficient of order a^. 

It is common to both calculations of the order corrections that the muon propagator 
is considered in the framework of the effective theory where the W boson is integrated out. 
The QED corrections to the resulting Fermi contact interaction were shown to be finite 
to all orders [|176|| . It is quite advantageous to perform a Fierz transformation which for a 
pure V ~ A theory has the consequence that afterwards the two neutrino lines appear in 
the same fermion trace. Thus the QED corrections only affect the fermion trace involving 
the muon and the electron. This also provides some simplifications in the treatment of 
75 since in the case of vanishing electron mass a fully anticommuting prescription can be 
used. 
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Figure 5.11: LO (dotted), NLO (dashed) and NNLO (sohd) results of R'{s), M = Mt = 
175 GeV. 



The decay rate can be written in the formQ 



where 



r = 2Mlm[zS: 



OS 
V 



qOS 



(5.3^ 



S' 



OS 



OS ryOS 



are functions of the variable 



qOS 



z. 



OS 



(1 + ^'v) 



(5.39) 



(5.40) 



and Ey represent the scalar and vector part of the muon propagator. They are 
functions of the external momentum q and the bare mass m°. In our case they further 
depend on the bare electromagnetic coupling a° and are proportional to the square of the 
Fermi coupling constant, G'j^. Typical diagrams contributing to and Sy are depicted 

can be found 



^ ^Som e discussion of Eq. (|5.38|) 
in tml . 



in particular the additional factor z in front of Sy^ 
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V^^j^ V^^y^ 



Figure 5.12: Sample diagrams for the muon (thick line) self energy. Two of the thin 
lines represent the electron and the corresponding neutrino and the third one the muon 
neutrino. All (one-particle-irreducible) diagrams involving the coupling of the photon 
(wavy lines) to the muon and electron have to be taken into account. 



in Fig. 5.12 . and represent the wave function and mass renormalization in the 
on- shell scheme. 

It is convenient to parameterize the QED corrections for the muon decay in the fol- 
lowing form 

- ro L(o) , ^ 4(1) , i^^V f 4(2) , 4(2) , 4(2) A , 



r(/i 



, (5.41) 



with 



G|.M^/(1927r^). a = a{fi) represents the electromagnetic coupling in the 
MS scheme and A^^]. corresponds to the contribution involving a virtual r loop. It is 



suppressed by M^/M^ and almost four orders of magnitudes smaller than the other 

1(2) 



terms ||178|| . The hadronic contribution is denoted by |p.7S 



To lowest order the result is known for finite electron mass, 



1 



Ml 
Ml 



m5 m5 



Ml 



Ml 



and in the limit Mp 



we obtain for the one-loop corrections 
25 _ 



4^^ 



-1.810. 



(5.42) 



(5.43) 



The approach chosen in |T8| to evaluate Aj^^^ has some similarity to the computation of 
the four-loop contribution to the (3 function (see also Section p73|) . In fact, since only 
the imaginary part has to be calculated one is only interested in the pole part like in 
the case of a MS renormalization constant. However, only the pole part arising from the 
cuts through the electron and neutrino lines have to be taken into account and not the 
ones through the muon line. This is ensured by considering the on-shell muon propagator 
which requires the evaluation of (the imaginary part of) four-loop on-shell integrals with 
external momentum = M^. Integration-by-parts relations are applied to reduce 
the number of occuring integrals to a small set of so-called master integrals. Only for 
the latter a hard calculation is necessary. Usually a large number of terms is generated 
in intermediate steps while using integration-by-parts relations. However, those four-loop 
integrals that have no imaginary part can immediately be discarded. This includes four- 
loop vacuum graphs and diagrams with a through-going on-shell line. A comprehensive 
discussion and lots of intermediate results can be found in Ref. [ITI 
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The approach of |T9[ used for the computation of the QED corrections of order a is 
based on an expansion of the full fermion propagator in the limit ^ g^, where q is the 
external momentum and is the on-shell propagator mass of the muon. The on-shell 
limit —>■ will be performed afterwards with the help of Pade approximations. This, 
of course, only provides an approximation to the exact result. However, the integrals to 
be evaluated are simplified considerably. Furthermore, the accuracy obtained with this 



method is sufficient to check the result of 181 and enables the same reduction of the 



theoretical error on Gp. Good convergence to the exact result is also expected from the 



analysis performed in ||179|1 where the hard-mass procedure has been applied to scalar 



two-loop integrals involving massless thresholds. 

Due to the Fierz transformation the loop integration connected to the two neutrino 
lines can be performed immediately as it constitutes a massless two-point function. This 
is also the case after allowing for additional photonic corrections. As a result one encoun- 
ters in the resulting diagram a propagator with one of the momenta raised to power e 
where D = 4 — 25 is the space-time dimension. This slightly increases the difficulty of 
the computation of the resulting diagrams. Especially for the order corrections, where 
the original four-loop diagrams are reduced to three-loop ones with non-integer powers 
of denominators, it is a priori not clear that these integrals can be solved analytically. 
However, it turns out that for the topologies needed in our case this is indeed possi- 
ble. For the computation of the massless two-point functions we have used the package 



MINCER [^. Only slight modifications are necessary in order to use this package also for 
the computation of the new type of integrals. 

In contrast to the current correlators considered in Section p]^ the expansion terms 
for z — ^ — oo can not be used as they describe the unphysical process fi ^ fi + j. Another 
difference is the presence of massless cuts in the limit — > 0. Thus a naive expansion 
is not possible and rather the asymptotic expansion has to be applied, which generates 
from the 44 contributing four-loop diagrams 72 sub- and cosub-diagrams that have to be 
evaluated. The analytical results are rather lengthly and cannot be listed. Instead the 
results are presented in numerical form . 



In order to get reliable results it is necessary to compute as many terms as possible in 
the expansion parameter z. Subsequently a Fade approximation is applied as described in 



Section |5.1| . We want to recall that before the Fade procedure a conformal mapping can 
be used which maps the complex z-plane into the interiour of the unit circle. Following 
Ref. [|177|| we denote those results by cj-Fades and the ones obtained without conformal 
mapping by z-Fades. 

The calculation is performed with the help of the package GEFICOM [^. It uses 
QGRAF for the generation of the diagrams and EXP ||180|| for the application of the 
asymptotic expansion procedures. 

The method has been successfully tested at Born level and at order a where a large 
number of moments can be evaluated. This gives a hint on how of many terms are 
necessary at 0(ct^) in order to obtain a reliable answer. The results for A^^^ of the 



93 



i nnnt, 


P.A. 


z 


UJ 


6 


[3/2] 


5.836 


7.249 


6 


[2/3] 


5.836 


7.057 


7 


[4/2] 


5.935 


7.040 


7 


[3/3] 


5.833 


7.076 


7 


[2/4] 


5.938 


7.080 


8 


[4/3] 


6.110 


6.873 


8 


[3/4] 


6.113 


7.060 


exact: 


6.743 



Table 5.4: Pade results for the corrections of 0{a'^) to the muon decay, A^^\ The first 
row indicates the order in z which has been used to construct the Fade approximations. 



individual Pade approximations are shown in Tab. |5 

= 6.5(7), 



It leads to the final answer |T9[ 

(5.44) 



where the deviation of the central value from the exact result of 6.743 is less than 3% 
and well covered by the extracted error of roughly 10%. Thus the sole knowledge of our 
results would also reduce the theoretical error on Gp. Finally, the total decay rate of the 
muon takes the form 



F(/i ^ u^eue) 



where A^^^^^ = -0.042(2) [TTSI and Ag[ = -0.00058 [TT^ has been used. Furthermore 



0.9998- 1.810 



a(M^ 



TT 



+ 6.700(2) 



'a(M^) 



+ . . . 



TT 



(5.45) 



fi^ = has been adopted. 

As already noted in ||1 75|| the numerical coefficient in front of the second order correc- 
tions becomes very smallQ if one uses the on-shell scheme for the definition of the coupling 
constant a. Then the MS coupling is given by q;(M^) = a{l + a/(37r) ln(M^/Mg )) and 
there is an accidental cancellation between the constant and the logarithm in the second 
order corrections. 

A similar kinematical situation as in the decay is also given for the semileptonic 
decay of a bottom quark, b — > ueug. From the technical point of view the difference is 
only due to the non-abelian structure of QCD. The two methods described above have 
been applied to the total rate 



r(6 ^ ueue) 



2.413 



71 



'a. (Mb 



+ 



TT 



(5.46) 



where Fg = G^M^^lKfel V(1927r3 



The order a? results read a9'^ = —21.296 111 8 111 and 



A^j^^ = —21.1(6) ||T9|. Again perfect agreement between the two methods is observed. 



26lnstead of "6.7" one has "0.27" in Eq. (|5^) 
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At the end of this Subsection we want to mention that similar methods have been 
used to compute the decay rate of a top quark into a W boson and a bottom quark ||177 . 



Contrary to the case of the muon the W boson is not integrated out from the Lagrangian. 
Thus, at order three-loop diagrams contributing to the top quark propagator have to 
be considered. In ||177|| a double expansion in q'^/M'^ and M^/M^ has been performed, 
where q is the external momentum of the top quark propagator, which leads to a reliable 
prediction for T{t Wb) including finite W-mass effects. 



5.4 The relation between the MS and on-shell quark mass 

In this Subsection we consider the relation between the on-shell and the MS quark mass 
at three- loop order in QCD. The result has been obtained for the first time with the 
help of the Pade method |[L24| ] and has been confirmed half a year later by a completely 



independent calculation [12(; 



Here, we want to describe the approach of [|124|| as it constitutes an other facet of 
applications of the method as described in Section 5J. In contrast to the examples 
described before one is directly interested in the real part of the considered function. 
Furthermore, in the relation between the MS and the on-shell value of the masses we 
want to know the function f{z) as defined in (|5.3|) for z = 1, whereas it is only available 
for small and large values of z. Nevertheless, the method is powerful enough to get the 
value for /(I) with an error of 2 — 3% (see below). 

The basic object entering the mass relation is the fermion propagator, S(q'). However, 
the Pade method cannot be applied directly to ^{q) as it contains (unknown) singularities 
at threshold. Thus, proper combinations have to be considered which are regular for z = 1. 
They are obtained from the requirement that the inverse fermion propagator has a zero 
at the position of the on-shell mass. 

In the following three different types of masses will appear: the bare mass, mo, the 
MS, m(/i) and the on-shell mass M. The relation between them is given by 

m(/i) = Z„mO = zMM, (5.47) 

where Zm is finite and has an explicit dependence on the renormalization scale /j,. It is 
the purpose of this Section to describe its calculation at order a^. 

As already mention above, in order to obtain the mass relation we have to consider 
the fermion propagator as shown in Eq. ( |2.15| ). The renormalized version can be cast in 
the form 

(SM)-' = i[{M-i) Sy{z) + M {zMSsiz) - Sv{z))] , (5.48) 

withQ 



Sv{z) = Z2(l + S° 



V) ) 

■ 



Ss{z) = Z2Z^{l-TPs). (5.49) 

the wave function 



^^Note that in contrast to the quantities defined in Eq. (|5.39| ) (see also Ref. 177 
renormahzation for functions S's/y is still defined in the MS scheme. 
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^2 denotes the wave function renormalization in the MS scheme which is sufficient for 
our considerations. Note that the functions Ss and Sy are MS quantities which later on 
are expressed in terms of the on-shell mass. The two-loop relation between m and M is 
enough to do this at order a^. It is convenient to write the functions Ss/v in the following 
way 



Ss/v = l + T.S^s/v(-) , (5-50) 



n>l 



where the quantities exhibit the following colour structures (the indices S and V are 
omitted in the following): 



5(1) = CpSp, 

'>FFH 
^FLH 

+ CfT^Sfhh. (5.51) 



S^"^^ = C'^Sff ~\~ CfCaSfa~\~ CfTtiiSfl ~\~ CfTSfh 1 
S^^^ = C^pSfff ~\~ C']^CaSffa~\~ CfC\Sfaa~\~ C^pTtiiSffl ~\~ CI^TSfi 
+ CfCaTuiSfal + CfCaTSfah + CfT'^u'^Sfll + CfT'^uiSi. 

I /^_T^2r< 



The same decomposition also holds for the function Zm- In (|5.51| ) rii represents the number 
of light (massless) quark flavours. Cf and Ca are the Casimir operators of the fundamental 
and adjoint representation. In the case of SU (iVc) they are given by Cf = {N^ ~ 1)/ (2iVc) 
and Ca = N^. The trace normalization of the fundamental representation is T = 1/2. 
The subscripts F, A and L in Eq. (|5.51|) shall remind us on the colour factors Cf, Ca 



and Trii, respectively. H simply stands for the colour factor T. 

A formula which allows for the computation of the MS-on-shell relation for the quark 
mass is obtained from the requirement that the inverse fermion propagator has a zero at 
the position of the on-shell mass: 



(SFiq))-' 



= 0. (5.52) 

=M2 



In the literature there are two different approaches to compute the occuring Feynman 
diagrams. In the first evaluation of the three-loop MS-on-shell relation the method of 
Section ^]l] has been applied which we will discuss below. On the contrary, a subsequent 
analysis | ]1 26| | has chosen = from the very beginning which makes it necessary to 
solve three-loop on-shell integrals. This can effectively be done using the integration-by- 
parts method within dimensional regularization [0. It enables the derivation of recur- 
rence relations which express complicated integrals in terms of simpler ones. At the end 
one arrives at a small set of integrals — so-called master integrals — which actually have 
to be evaluated. In |p.26| , p.82|| the considerations of ||183|| have been extended and the 



missing master integrals have been evaluated. The technique used for the computation is 
based on the hard-mass procedure for large M which represents the on-shell integrals in 
terms of a power series in q'^/M'^. The coefficients contain nested harmonic sums which in 
the on-shell limit, i.e. for = M^, can be reduced to known mathematical constants. For 
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explicit examples we refer to [ p.82|| . At the end of this section we will list the analytical 



result for z^ifi) obtained in ||126| . 



The starting point for the approach of |124, 125 1 is Eq. ( 5.52|) . Applying it to Eq. ( 5.48 
leads to the condition 

h{z) = Zmifi') Ss{z) — Sv(z) =0 for z = l. (5.53) 



At a given loop-order L, Eqs. ( ^.49|) are inserted and the resulting equation is solved for 
zl^\ Thus Eq. ( ^.53|) can be cast in the form 



h{z) = f{z) + zl^^i^^y . (5.54) 

Our aim is the computation of /(I). Note that the individual self energies and Ey 
develop infra-red singularities when they are evaluated on-shell. The proper combination 
which leads to the relation between the MS and on-shell mass is, however, free of infra-red 
problems. 

The explained procedure has been applied to each colour structure occuring in f{z) 
separately as outlined in [ |125|| . E.g., collecting all terms proportional to Cp in Eq. ( |5.54| ) 
leads to 

fpFFiz) = Ss,FFf{z) - Sv,FFf{z) + Z^Ss,FF + Z^^Ss,F , (5.55) 

which has to be evaluated for z = 1. The individual terms on the right-hand side develop 
a (unknown) singular behaviour which is encoded in the corresponding moments. In case 
the Pade method is (naively) perfomed with the individual pieces, the Pade approximants 
try to imitate the threshold singularity. However, due to the very construction of the 
method the analytical structure of the Pade result is polynomial for z — > 1 and the typical 
threshold logarithms can not be reproduced. Thus, the results show instabilities in the 
vicinity of z = 1 . On the contrary, the proper combination as given in Eq. ( |5.55| ) has to be 
regular for 2; ^ 1, as the on-shell mass does not contain any infra-red singularity ||184| , |185 . 
The corresponding Pade results demonstrate great stability. 

To summarize, although superficially only information about small and large momenta 
enter the Pade procedure, it is sensitive to the analytical structure at threshold as the 
information about the singularity is, to some extend, also contained in the moments of 
the analytical function f{z). 

At this place we will refrain from the discussion of the individual colour structures 



which can be found in [p.25|| but only present the results for the sum where the numerical 
values for Ca, Cp, T and rii have been inserted. Still care has to be taken because of the 
diagrams involving a closed heavy quark loop. 



It was already realized in [p.34|| that the Pade procedure shows less stability as soon 
as diagrams are involved which exhibit more than one particle threshold. In our case the 
interest is in the lowest particle cut which happens to be for = M^. The Pade method 
heavily relies on the combination of expansions in the small and large momentum region. 
The large momentum expansion, however, is essentially sensitive to the highest particle 
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n 


P.A. 


= 


ni = \ 


n/ - 2 


ni = 3 


n/ = 4 


n/ = 5 


5 


[4/5] 


-200.2787 


-173.6663 


-148.3787 


-124.4156 


-101.7771 


-80.4628 


5 


[4/6] 


-201.6419 


-174.8844 


-149.4553 


-125.3553 


-102.5862 


-81.1625 


5 


[5/4] 


-203.9394 


-176.7290 


-150.8970 


-126.4411 


-103.3591 


-81.6482 


5 


[5/5] 


-201.4721 


-174.7445 


-149.3422 


-125.2644 


-102.5104 


-81.0786 


5 


[5/6] 


-198.7799 


-172.8884 


-148.2001 




-102.2739 


-81.0336 


5 


[6/4] 


-202.8435 


-175.8651 


-150.2298 


-125.9387 


-102.9929 


-81.3939 


6 


[4/6] 


-201.0906 


-174.3880 


-149.0165 


-124.9749 


-102.2619 


-80.8758 


6 


[5/5] 


-200.9265 


-174.2458 


-148.8927 


-124.8668 


-102.1673 


-80.7929 


6 


[5/6] 


-200.4927 


-173.9600 


-148.7433 


-124.8358 


-102.2290 


-80.9131 


6 


[5/7] 


-200.3603 


-173.7018 


-148.3764 


-124.3940 


-101.7753 


-80.5533 


6 


[6/4] 


-201.6970 


-174.9293 


-149.4861 


-125.3673 


-102.5725 


-81.1016 


6 


[6/6] 


-200.3195 


-173.6857 


-148.3751 


-124.3879 


-101.7244 


-80.3848 


6 


[7/5] 


-202.1300 


-175.2569 


-149.7173 


-125.5125 


-102.6443 


-81.1143 



Table 5.5: Pade results for the sum of those contributions which don't have a closed 
heavy fermion loop, n/ has been varied from to 5. 



threshold. Thus, if this threshold numerically dominates the lower-lying ones it cannot be 
expected that the Pade approximation leads to stable results. In such cases a promising 
alternative to the above method is the one where only the expansion terms for ^ are 
taken into account in order to obtain a numerical value at = M^. This significantly 
reduces the calculational effort as the construction of the Pade approximation from low- 
energy moments alone is much simpler. In practice this approach will be applied if the 
Pade results involving also the high-energy data look ill-behaved. 

In the present analysis diagrams with other cuts than for = are already present 
at the two-loop level (see Fig. pl^ ) which allows us to test these ideas. Indeed, taking into 
account terms up to order and performing a Pade approximation there is an agreement 
of four digits with the exact result [p.25|| . Also at three-loop order either = or 

= 9M^ cuts appear. Cuts involving five or more fermion lines are first possible starting 
from four-loop order. Note that cuts involving an even number of fermions cannot occur. 

Concerning the colour structures introduced in Eq. (^.51 ) we use for the sum of the 
structures FFF, FFA, FFL, FAA, FAL and FLL both the low- and high-energy mo- 
ments whereas for the sum of the structures FFH, FAH, FLH and FHH only the 
expansion for 2; — is used. In Tabs. and |5.6| the results for different Pade ap- 
proximations are listed, n indicates the number of low-energy moments involved in the 
analysis, i.e. n = 6 implies the inclusion of terms of 0{z^). The number of high-energy 
terms can be obtained in combination with the order of the Pade approximant {[x/y\) 
and is given hy x + y + 1 — n. 



The final result of ||124| , |125|| for the mass relation can be found in Tab. where a 
comparison with the results of [|126|| is performed. Note that there is perfect agreement 
for all values of ni. At this point we want to mention that the result of |p.24| , |125|| is 
more general as the function f{z) in Eq. (|5.54|) has been computed for all values of z and 
not only for the special point z = 1. This opens the possibility to obtain the fermion 
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n 


P.A. 




71; = 1 


ni = 2 


ni = 3 


71; = 4 


71/ = 5 


4 


[1/3] 


-0.9345 


-0.9572 


-0.9798 


-1.0024 


-1.0249 


-1.0475 


4 


[2/2] 


-0.9321 


-0.9546 


-0.9770 


-0.9995 


-1.0218 


-1.0442 


4 


[3/1] 


-0.9324 


-0.9551 


-0.9777 


-1.0003 


-1.0229 


-1.0455 


4 


[4/0] 


-0.9604 


-0.9828 


-1.0053 


-1.0277 


-1.0501 


-1.0725 


5 


[1/4] 


-0.9271 


-0.9495 


-0.9720 


-0.9944 


-1.0169 


-1.0393 


5 


[2/3] 


-0.9219 


-0.9440 


-0.9661 


-0.9882 


-1.0103 


-1.0324 


5 


[3/2] 


-0.9086 


-0.9347 


-0.9591 


-0.9827 


-1.0060 


-1.0290 


5 


[4/1] 


-0.9254 


-0.9478 


-0.9703 


-0.9927 


-1.0151 


-1.0375 


5 


[5/0] 


-0.9495 


-0.9719 


-0.9942 


-1.0166 


-1.0389 


-1.0613 


6 


[1/5] 


-0.9217 


-0.9441 


-0.9665 


-0.9888 


-1.0112 


-1.0336 


6 


[2/4] 


-0.9140 


-0.9364 


-0.9589 


-0.9813 


-1.0037 


-1.0261 


6 


[3/3] 


-0.9125 


-0.9352 


-0.9578 


-0.9803 


-1.0028 


-1.0252 


6 


[4/2] 


-0.9126 


-0.9352 


-0.9578 


-0.9803 


-1.0028 


-1.0253 


6 


[5/1] 


-0.9202 


-0.9425 


-0.9649 


-0.9872 


-1.0096 


-1.0319 


6 


[6/0] 


-0.9416 


-0.9639 


-0.9862 


-1.0085 


-1.0308 


-1.0532 



Table 5.6: Pade approximations performed in the variable z. No high-energy results 
have been used. Again ni has been varied from to 5, the dependence, however, is very 
weak. 





Zm(M) - "1 M)/M 






"1 




0(af) [124 






CI(q^) |124| 


O(a-l) [L26| 




0(ai) 124] 






1 

2 
3 
4 

5 


-14.33 
-13.29 
-12.25 
-11.21 
-10.17 
-9.13 


-209T55 

-176(4) 
-150(3) 
-126(3) 
-103(2) 
-82(2) 


-i^rf- 

-172.4 
-147.5 
-123.8 
-101.5 
-80.40 


-11.67 
-10.62 
-9.58 
-8.54 
-7.50 
-6.46 


-ivlJ(!i! 

-146(4) 
-123(3) 
-101(3) 
-81(2) 
-62(2) 


-iBT? 

-142.5 
-120.0 
-98.76 
-78.86 
-60.27 


13.44 
12.40 
11.36 
10.32 
9.28 
8.24 


ijSTO 

168(4) 
143(3) 
119(3) 
96(2) 
75(2) 


im% 

164.6 
139.9 
116.5 
94.42 
73.64 



Table 5.7: Two- and three- loop coefficients of the relation between on-shell and MS mass. 
The choice /i^ = M^, respectively, /i^ = has been adopted. 



propagator in QCD at three-loop order for arbitrary external momentum extending the 



considerations of |p.86| ] by one more loop. 

For the relation between the MS and on-shell quark mass one finds up to three 

p3|, pmpRra 



loops 

m(/i) 
M 



1 + 



TT 



4 
3 



445 19 2 



+ 



71 



2 ^ 



TT 



3019 , 2, , 1^ 
288 ^ 3^ 6^ 



(144 



13 



+ i^C2 + —l^M + —llM)ni-- ^ 



12 



M 



+ 



TT 

11779^, 



475 



165635 25 , 25 , , 55 ^ , , 

C2 C2ln2 H Cs Lm 

2592 3 ^ 9 ^ 36^ ' ^ 



864 ^"^^ 432 



/ 10051 37 



7 



1 H C2 + -C2 In 2 + -Cs Lm + 

V 1296 18^ 9^ 9^^ ' ^ 



911 

432 W 



99 



11,3 

+ 54 



+ nf 



89 
'648 



gC2 j lUi 



(5.56) 



where (2 = vr^/6 and I^m = In/x^/M^. The constant z^\M) is given by [|126|] 



.(3) 



(M) 



9478333 55 4 
^ 162 ^ 



644201 587 44 
6480 ^ ^ "^^^ 27 



93312 

3475, 1439,, 1975, 220 



61 



432 ' 72 ' ' 216 ' 27 

/967 22 , 4 , 2 \ , 241 

+ + — n 2 n^ 2 C2 + 

V1O8 27 27 )^ 72 

2353 13, 7 , 
77C2 - 77C3 



246643 1,4^ 

in 2 

23328 81 



A 305 

C3 - Y5^C4 



27 



04 



23328 54 



54 



(5.57) 



where 04 = Li4(l/2) ^ 0.517479. 

The function A(r) in Eq. ( |5.56| ) takes into account the effects of secondary hght quarks. 
If < r < 1 then the function A(r) may be conveniently approximated as follows [ |123|| 



TT 



0.597 + 0.230 



(5.5^ 



which is accurate to 1%. Up to now there is no calculation available taking into ac- 
count the complete mass dependence of the light quarks at order a^. The subclass of 
diagrams containing a one-loop light quark vacuum polarization insertion has been con- 
sidered in ||187|| , where it was observed that the dominant contribution is provided by the 
linear mass corrections as at order a1 (cf. Eq. ( p.58| )). 

At this point it is worthwhile to compare the results of [|124| , |125| , |126|| with estimations 
for the 0{al) terms obtained with the help of different optimization procedures. In 
the fastest apparent convergence (FAC) ||188|| and the principle of minimal sensitivity 
(PMS) ||189|| have been used in order to predict the three-loop coefficient of M/m{m). 
For 111 = 2 one observes a discrepancy of 9%. It even reduces to below 1% for rii = 5, i.e. 
in the case of the top quark. The results obtained in the large /?o-l™it [p-90|| , where Pq is 
the first coefficient of the QCD /? function, agree to 2% for rii = 3, roughly 8% for n; = 4 
and 17% for n/ = 5. 

Among the various applications of the order term in the MS-on-shell relation we 
only want to mention the improvement in the determination of the top quark mass to be 
measured at a future e~^e~ linear collider. To be specific, let us consider the production of 
top quarks in e~^e~ collisions. The corresponding physical observables expressed in terms 
of Mt show in general a bad convergence behaviour. In the case of the total cross section, 
e.g., the next-to- next-to-leading order corrections partly exceed the next-to-leading ones. 
Furthermore the peak position which is the most striking feature of the total cross section 
and from which finally the mass value can be extracted depends very much on the number 
of terms one includes into the analysis. The commonly accepted explanation for this is 
that the pole mass is sensitive to long-distance effects which result in intrinsic uncertainties 
of order Aqcd [P-91| , |192|| . In other words, it is not possible to determine the pole mass 



from the analysis of the cross section at threshold with an accuracy better than Aqcd- 
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Several strategies have been proposed to circumvent this problem ||193| , [L94| , |195|| . They 



are based on the observation that the same kind of ambiguities also appear in the static 
quark potential, V{r). In the combination 2Mt + V{r), however, the infra-red sensitivity 
drops out. Thus a definition of a short-distance mass extracted from threshold quantities 
should be possible. The relation of the new mass parameter to the pole mass is used 
in order to re-parameterize the threshold phenomena. On the other hand a relation of 
the new quark mass to the MS mass must be established as it is commonly used for the 
parameterization of those quantities which are not related to the threshold. In order to do 
this consistently the three-loop relation between the MS and the on-shell mass is needed. 

In ||193|| the concept of the so-called potential mass, rrit^ps, has been introduced. Its 
relation to the MS mass, mtlmt), reads 

mt,P5 (20 GeV) = (165.0 + 6.7 + 1.2 + 0.28) GeV , (5.59) 

where the different terms represent the contributions from order q;° to af. For the nu- 
merical values mt{mt) = 165.0 GeV and a^^\mt{mt)) = 0.1085 have been used. The 
comparison of Eq. ( |5.59| ) with the analogous expansion for Mt, 

Mt = (165.0 + 7.6 + 1.6 + 0.51) GeV, (5.60) 

shows that the potential mass can be more accurately related to the MS mass than Mf. 
Further details and more examples can be found in ||125| . 
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A Technical remarks 
A.l Asymptotic expansion 

A promising approach to compute — at least in certain kinematical limits — multi-loop 
diagrams is based on asymptotic expansion. An asymptotic expansion can be considered 
as a generalization of a Taylor expansion. In both cases one obtains an expansion in powers 
of a small quantity. However, in case of the asympotitic expansion the corresponding 
coefficients are not constant but contain non-analytical functions of the small parameter. 
As a simple example let us consider the function f{x) = Li2(l — x) for which the Taylor 
expansion for x — does not exist beyond the leading order. Nevertheless there is an 
asymptotic expansion which reads 

f{x) = — + x(lnx-l) + C(x2lnx) , (A.l) 
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with the non-analytic (for x ^ 0) function Inx in the coefficient of the hnear term in x. 

In the case of Feynman diagrams the situation is similar. Systematic procedures have 
been developed in the case of a large external mometum ( "large-momentum procedure" ) 
and a large internal mass ("hard- mass procedure") ||129|| . Both procedures apply to prob- 
lems which can be formulated in Euclidean space. This is the case for all calculations 
presented in this review. In contrast to that there are also phenomena which are tightly 
connected to the Minkowskian space-time. Also in such cases rules for an asymptotic 



expansion have been developed. In particular, two- loop on-shell two-point diagrams [19C: 



and two-loop vertex diagrams in the Sudakov limit [|197|| have been considered. Also the 



threshold expansion ||198|| belongs to this class of phenomena. 

One may treat the large-momentum and hard-mass procedures on the same footing. 
Thus, in what follows we only present the general formulae in the case of large external 
momenta — the transition to the hard-mass procedure is straightforward. The prescrip- 
tion for the large-momentum procedure is summarized by the following formula^: 



T{Q,m,q) ^r/7(g,m) -k T^q^,ra^}'jiQ,m^,q^) . (A.2) 

7 

Here, F is the Feynman diagram under consideration, {Q} ({m, q}) is the collection of 
the large (small) parameters, and the sum goes over all subgraphs 7 of F with masses 
rriy and external momenta q^, subject to certain conditions to be described below. T^q^m} 
is an operator performing a Taylor expansion in {q, m} before any integration is carried 
out. The notation F/7 ^ T{g indicates that the subgraph 7 of F is replaced by its 
Taylor expansion which should be performed in all masses and external momenta of 7 
that do not belong to the set {Q}- In particular, also those external momenta of 7 that 
appear to be integration momenta in F have to be considered as small. Only after the 
Taylor expansions have been carried out, the loop integrations are performed. In the 
following we will refer to the set {7} as hard subgraphs or simply subgraphs and to {F/7} 
as co-subgraphs. 

The conditions for the subgraphs 7 are different for the hard-mass and large- 
momentum procedures]^. For the large-momentum procedure, 7 must 

• contain all vertices where a large momentum enters or leaves the graph 

• be one-particle irreducible after identifying these vertices. 

From these requirements it is clear that the hard subgraphs become massless integrals 
where the scales are given by the large momenta. In the simplest case of one large 
momentum one ends up with propagator-type integrals. The co-subgraph, on the other 
hand, may still contain small external momenta and masses. However, the resulting 
integrals are typically much simpler than the original one. 
In the case of hard-mass procedure, 7 has to 

^^In the case of the hard-mass procedure one essentially has to replace Q by the large mass M. 

Actually they are very similar and it is certainly possible to merge them into one condition using a 
more abstract language. For our purpose, however, it is more convenient to distinguish the two procedures. 
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(a) (b) 
Figure A.l: Fermionic double-bubble diagrams with generic masses mi and 



+ 2x 




Figure A. 2: Hard-mass procedure for the double-bubble diagram. The hierarchy <^ 
is considered where m is the mass of the inner line. The hard subdiagrams (right of "t^t" ) 
are to be expanded in all external momenta including q and reinserted into the fat vertex 
dots of the co-subgraphs (left of "t*:"). 



• contain all the propagators carrying a large mass 

• be one-particle irreducible in its connected parts after contracting the heavy lines. 

Here, the hard subgraphs reduce to tadpole integrals with the large masses setting the 
scales. The co-subgraphs are again simpler to evaluate than the initial diagram. 

An example demonstrating the practical application of the large-momentum expansion 



was already presented in Section^ (cf. Fig. see also ||116|| ). As an application of the 
hard-mass procedure let us consider the double-bubble diagram of Fig. |A.1| with the 
hierarchy ml <^ -C m|. The imaginary part leads to contributions for the total cross 
section cr(e'^e~ —>■ hadrons). One may think of charm quark production (mi = Mc) in the 
presence of a virtual bottom quark (m2 = Mb). It turns out that already the first term 
provides a very good approximation almost up to the threshold ^/s = 2M5 ||199| , |20CI| , |201 
For simplicity we set mi = and m2 = m in the following. 



The corresponding diagrammatic representation is shown in Fig. [A.2|. There are three 
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sub diagrams, where one again corresponds to the naive Taylor expansion of the integrand 
in the external momentum q. After Taylor expansion, the subdiagrams are reduced to tad- 
pole integrals with mass scale m. The scale of the co-subgraphs is given by g, thus leading 



to massless propagator-type integrals. The result for the first three terms reads |202, 116 
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with Iq^ = ln(— g^//i^) and /qm = ln(— g^/m^). 

We want to stress that the main simplification, which is common to all kinds of 
asymtotic expansions, comes from the fact that the expansions in the small parameters 
are done before any momentum integration is performed. The proof that this leads to 
correct results is based on the so-called strategy of regions p03| ] • There different regions of 
each loop momentum are selected and in each of them Taylor expansions with respect to 
the small parameters are performed. In the limits of the hard-mass and large-momentum 
procedures an interpretation of the different regions in terms subgraphs and cosub-graphs 
is possible (see above). This is different in the case of the threshold expansion ||198|| where 
a graphical representation becomes much less transparent. However, the application of 



the strategy of regions [203] leads to correct results. 



A. 2 Single-scale Feynman diagrams up to three loops 

The problem of evaluating one-loop Feynman diagrams is — at least in principle — solved 
(see, e.g., Refs. [0, |^, |204| , ^). However, in case many legs and lots of different masses 



appear also one-loop computations can become very tedious, in particular if degenerate 
momentum configurations are involved. 

At two-loop order the class of Feynman diagrams which have been studied in detail 
is much more restricted. There is a good understanding of two-point functions (see, 
e.g., |]T0[) which also has found important physical applications ||205|| . Concerning three- 



and four-point functions one is essentially restricted to the massless case. In this context 
we want to draw the attention to the recent activity in the computation of the two-loop 
box diagrams (for a brief overview see |jll[ ) . Within the last two years the basis has been 
estabhshed to compute two-loop virtual corrections to the four-point Feynman amplitudes 
where all internal lines are massless and at most one external leg is off-shell. This opens 
the door to investigate next-to-next-to-leading order processes like the two-jet production 
at hadron colliders or three-jet production in e'^e~ annihilation. 

It is obvious that the complexity of the computation of a Feynman diagram strongly 
depends on the number of different scales involved. There is one class of diagrams which 
is studied in great detail up to three loops, namely diagrams which only depend on 
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one dimensionful scale. Next to massless propagator-type diagrams with one external 
momentum, g, we have in mind vacuum integrals with one non-zero mass, M and so- 
called on-shell integrals for which the condition = is fulfilled. 

The basic idea for the computation of the integrals is common to all three types: after 
the numerator is decomposed in terms of the denominator recurrence relations are applied 
which express the diagram as a linear combination of so-called master integrals. Only for 
the latter a hard computation is necessary. However, since in the case of single-scale 
diagrams the master integrals are essentially pure numbers it is also possible to use high- 
precision numerical methods in case an analytical calculation is not possible. We want 
to mention that for the massless propagator-type and the massive vacuum integrals two 
and nine master integrals are needed, respectively. In the case of the three-loop on-shell 
integrals a list of all master integrals can be found in the Appendix of Ref . ||182|| . Counting 
also those integrals which are composed of products of lower-order diagrams they amount 
to 18. 

For convenience we want to provide the analytical results for the massless one-loop 
two-point functions, Pab{Q), the on-shell two-point functions, Oab{Q), and the one- (K) 
and two-loop (Vabc) vacuum integrals in Euclidian space. 



PabiQ) 



d^p 1 



Oab{Q) 
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The results of a general three-loop diagram can not be expressed in terms of F func- 
tions. Moreover, due to the large number of contributing diagrams and the complexity 
of intermediate expressions it is absolutely necessary to use computer algebra programs 
for the computation of multi-loop diagrams. It is thus also hardly possible to provide 
intermediate results which eventually could be used in other calculations. Therefore, on 
one side one is left with the description of the method used for the evaluation of the dia- 
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grams. On the other hand it is possible to provide the program code which was used for 
the computation. Thus everybody can repeat the calculation or apply it to own problems. 

Both for the massless propagator-type integrals and the massive vacuum integrals 
FORM ||5^ packages have been published: massless integrals up to three loops can be 
computed using MINCER |^T[; the package for the massive integrals is called MAT AD In 
the following we will present an example which demonstrates the use of MATAD. The use 
of MINCER is very similar. Actually, in the package GEFICOM ^ MATAD and MINCER are 
used in parallel using the same notation for the input. 



Let us consider the triangle diagram as pictured in Fig. |A.3| . It is one of the 657 
diagrams which contribute to the coefficient function C° appearing in the effective La- 
grangian of Eq. ( [4.13|) . According to the Lorentz structure the result can be written 

as 



where qi and q2 are the momenta of the gluons with polarization vectors e'^(gi) and e^{q2). 
Thus the vertex diagrams have to be expanded up to linear order both in qi and and 
an appropriate projector has to be applied in order to get K{Mt). 




Figure A. 3: Sample diagram contributing to the decay of the Higgs boson, 
looped lines represent quarks and gluons, respectively. 



Solid and 



MATAD requires one file containing the diagrams and the projectors which have to be 
applied. In our case it could look as follows 

*— #[ TREATO: 
multiply, ( 

a*deno(2 , -2) * (ql .q2*d_ (mu,nu) -q2 (nu) *ql (mu) -q2(mu) *ql (nu) ) 
+b*deno (2,-2)* (-ql . q2*d_ (mu , nu) + (3-2*ep) *q2 (nu) *ql (mu) +q2 (mu) *ql (nu) ) 

); 

. sort 

*— #] TREATO: 

*— #[ TREATl: 
*— #] TREATl: 

*— #[ TREAT2: 
*— #] TREAT2: 
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*— #[ TREATMAIN: 
*— #] TREATMAIN: 



*— #[ d31335: 
((-1) 
*M 

*Dg(iiul,nu2,pl) 
*Dg(iiu7,iiu8,-p4) 
*Dg (nu3 , nu4 , ql , -pi ) 
*Dg (nu5 , nu6 , q2 , pi ) 

* S ( -ql , -p3m , iiu7 , -ql , p5m , nu4 , -p2m , nu6 , q2 , p5m , nu8 , q2 , -p3m) 
*V3g(mu,ql ,iiul ,-pl ,nu3,pl-ql) 
*V3g(nu,q2,nu2,pl ,nu5,-pl-q2) 
*1) ; 



ttdefine TOPOLOGY "04" 
*— #] d31335: 



where the diagram d31335 corresponds to the one shown in Fig. |A.3| (for details concerning 
the nomenclature of the momenta see |^5[). The fold TREATO contains (up to an overall 
factor {qi ■ ^2)"^) the projector on the coefficients in front of the structures g^'^ and q^q!^ of 
Eq. (|A.8|) . They are marked by the symbols a and b, respectively. Thus the transversality 
of Eq. (|A.8|) can be explicitly checked in the sum of all contributing diagrams (the result 
of a single diagram does in general not have a transverse structure). 
A second file which is required, the so-called main-file, reads 



#define PRB "hgg" 
#define PRDBLEMO "1" 
#define DALA12 "1" 
#define GAUGE "0" 
#define POWER "2" 
#define CUT "0" 
#define FOLDER "hgg" 
#define DIAGRAM "d31335" 
#- 

Sinclude main. gen 



The fifth line ensures that an expansion of the integrand up to the second order in the 
external momenta is performed and the third line sets qf and g| to zero and factors out 
the scalar product qi ■ q2- #def ine CUT "0" sets e to zero in the final result. In this 
example we choose Feynman gauge which is achieved with #def ine GAUGE "0". 
After calling MAT AD it takes of the order of a minute to obtain the result: 

d31335 = 

+ ep'-2 * ( 40*qi.q2*M'2*a + 344/9*Ql . Q2~2*a - 232/9*Ql . Q2"2*b ) 

+ ep--l * ( - 308/3*Ql.Q2*M-2*a - 3530/27*Ql . Q2-2*a + 1786/27*Q1 . 02^2* 
b ) 
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+ 60*qi.q2*M'2*z2*a + 734/3*qi . q2*M'2*a - 1936/9*qi . q2'2*z3*a + 1136/9* 
qi.q2-2*z3*b + 172/3*qi.q2-2*z2*a - 116/3*Q1 . Q2-2*z2*b + 46817/81* 
qi.q2-2*a - 26239/81*Ql.Q2-2*b; 

Note that the terms proportional to Q1.Q2*M~2 cancel after adding all contributing dia- 
grams. 



A. 3 Automation of Feynman diagram computation 

In this section we briefly want to mention the program packages which have been used to 
obtain most of the results discussed in this review. For a general overview concerning the 



automation of Feynman diagram computation we refer to ||116 . 

The large number of diagrams which occurs in particular if one considers higher loop 
orders makes it necessary to generate the diagrams automatically. The Fortran program 
QGRAF provides the possiblility to implement own models in a simple way. Furthermore 
it is quite fast and generates several thousand diagrams in a few seconds. One of the 
disadvantages of QGRAF is that the user has to put the Feynman rules himself. On the 
other hand, this provides quite some flexibility in the choices of the vertices. E.g., it is 
straightforward to implement the vertex involving the coupling of the operator Oi to four 
gluons (cf. Eq. (^)). 

In general, the application of asymptotic expansions, in particular if serveral of them 
are applied successively, generates many subdiagram which have to be expanded in several 
small quantities. Even for a single multi-loop diagram this becomes very tedious if it has 



to be performed by hand. For this reason the programs LMP [^] and EXP | 180 | have 
been developed. LMP was especially developed in order to apply the large-momentum 
procedure to the diagonal current correlators [|136| , |137| , |95| , |110|1 . In some sense EXP can 



be considered as the successor of LMP. Next to the hard-mass procedure also the succesive 
application of large-mometum and/or hard-mass procedure is possible. This broadens the 
area of applications. Here we just want to mention as examples the correction of 0{aas) 
to the Z boson decay ||206|| , O^a"^) corrections to the top quark decay [ p.771| or two-loop 



QED corrections to the muon decay [0 (cf. Section 

The very computation of the integrals is performed with the program packages 
MINCER and MATAD [gSf (see Appendix The former deals with massless 



propagator-type integrals up to three loops and MATAD was written to deal with vac- 
uum diagrams at one-, two- and three-loop order where several of the internal lines may 
have a common mass. The area of application for each of the individual packages seems 
to be quite restrictive. However, in particular the combined application offers a quite 
flexible use. 

In order to handle problems where a large number of diagrams are involved and where 
eventually an asymptotic expansion has to be applied in a convenient way the program 
package GEFICOM has been written. A very limited number of small input files allows the 
user to rule the flow of the computation. Qgraf is called to generate the diagrams. A 
Mathematica |^2[ script determines the toplogy of each individual diagram and provides 



input which either can be directly read from MINCER and/or MATAD or can be passed to EXP 
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or LMP. At the end the results of the individual diagrams are summed and the bare result 
is stored. Moreover a convenient environment is provided which, e.g., makes sure that all 
result files are up-to-date. Thus, processes involving a large number of (sub-) diagrams 
can be treated without taking care of each individual result. 



B Decoupling constants and coefficient functions 

Transforming the decoupling constants of Eqs. (|3.28|) and (|3.30|) to the MS scheme one 
obtains 
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In the following, we list the decoupling constants Cm and Q appropriate for the general 
^auge group SU(A^c)- The results read 
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For Ai'c = 3, we recover Eqs. ( |B.1| ) and ( |B.2[ ). 

The renormalized decoupling constants C2 and ^3 for the quark and gluon fields, respec- 
tively, arise from Eqs. ( |3.15| ) and ( |3.10D . Of course, C2 and ^3 are both gauge dependent. 
Restricting ourselves to the case iVc = 3, we find in the covariant gauge ( p.2|) 
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Figure B.l: Feynman diagrams contribution to Z2, Zm and . The dashed hne either 
represents the Higgs boson (h) or the neutral (x) or charged (0^) Goldstone boson. 



-MS 




121 /i2 13 , 2 1 3 /i2 ^ 

In - — In^ At + — In^ At 



(B.5) 



91 



1152 96 



29 3,2/^ 
+ t::: In At + — In^ At 



32 



86183, 99 , 1 
27648^^ 128^^ 32 ^ 



+ 



V 27648 64 



- — Cs ) In — + 



383 
2304 



, 2 119, 3 /i' 

In^ + — In'^ At 



/ 3307 
+ nii 



15552 
677 
" 1536 



12^^" 
32^^ 



293 , /i 



768 



mt 



1152 
233 



In + 737 In At - — In'' At 



m 



h 



36 



96 



mt 



In4-Ain< 
1024 ml 32 ml 



256 m|^ 



m 



(B.6) 



In Ref. [g9l the Yukawa corrections to Cm and (g enhanced by the top quark mass have 
been evaluated. They are conveniently expressed in terms of the variable 
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where mt is the top quark mass defined in the MS scheme. Corrections proportional to Xt 
arise if in addition to the pure QCD Lagrangian also the couplings of the Higgs boson (h) 
and the neutral (x) and charged (0^) Goldstone boson to the top quark are considered. 
Sample diagrams contributing to (m and Q are shown in Figs. p.l| and p.2| , respectively. 
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Figure B.2: Feynman diagrams contribution to Z3 and (g. The dashed hne either 
represents the Higgs boson {h) or the neutral (x) or charged (0^) Goldstone boson. 



The decouphng constant for the u, d, s and c quark mass read^ 
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where I31 is the third component of the weak isospin, i.e. I31 = +1/2 for up-type quarks 
and I31 = —1/2 for down- type quark flavours. For the bottom quark one receives 
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where we have used Cp = 4/3, Ca = 3 and T = 1/2. The constants 
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For convenience, the corrections of order are repeated from Eq. ( B.l 
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where ^4 = CI2 is Clausen's function and Li4 is the quadrilogarithm, occur in the 

evaluation of the three-loop master diagrams |207| , p08| , p09|| . 
Finally, for Q we obtain the following result: 
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where the contribution of the diagrams in Fig. |B.2| (c) corresponds to the last entry in the 
last line of Eq. (|RTTD . For convenience also the pure QCD result of 0{al) is listed. The 
corresponding three-loop terms can be found in |^ . 



In the remainder of this Appendix we want to provide the analytical results for the 
Xj-enhanced corrections of order Gpm^ to the coefficient functions Ci and C2 |^9|]. For 
Ci we have 
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where mt is the MS top quark. The 0{asXt) terms can be found in 2101 , 75| and the 0{al 
results were computed in [|100| , |101|] The corrections of O(asX^) are taken from fS^ 
For the light quarks we get for C2 
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and in the case of the bottom quark the coefficient function reads: 

C2b = C2d + Xf 
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In |£5l C2i and C2b are hsted for general gauge group SU{Nc). 
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C Analytical results for R{s) 



As in the literature the quadratic and quartic correction terms to R{s) are not yet avail- 
able in analytic form we want to provide the corresponding results using the notation 
introduced in Eqs. ( [4.73|) and ([4.74|) . For completeness also the massless approximation 
is given. Of course, it is the same in the MS and on-shell schemes 
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with a4 = Li4(l/2) ^ 0.517479. 
Note added: 

In the meantime a third independent evaluation of the order QED corrections to the 
muon decay became available |pll|| . 
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